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Introduction 



Nowadays, it is a well known fact that most of the matter in our Solar System, 
in our Galaxy and, probably, within the whole Universe, exists in the form of 
ionized particles. On the other hand, in our natural environment, that is to say in 
'normal' conditions as on our little planet, gaseous matter generally consists of 
neutral atoms and neutral molecules: only under specific conditions, in fact, like 
along lightnings' paths, a large amount of both atoms and molecules can be found 
as ionized. For decades, from this point of view, astronomers and astrophysicists 
believed that only positively charged ions were worthy of relevance in drawing 
the networks for possible chemical reactions in the interstellar medium, as well as 
in modeling the physical conditions in most of astrophysical environments. Thus, 
negative ions (and especially molecular negative ions) received minor attention 
until their possible existence was observationally confirmed (discovery of the first 
interstellar anion, CeH"), about thirty years after the first physically reasonable 
proposal on their actual detection was theoretically surmised by E.Herbst. 

From a purely theoretical point of view, negatively charged ions play a peculiar 
role as they can be formed in large quantities in the gas phase by attachment of 
low (1.5 eV < < 10 eV) and very low-energy (E < 1.5 eV) free electrons, 
while sometimes such a formation process can occur even close to zero eV. In an 
astrophysical context, roughly speaking, their role should be then found in their 
involvement in the charge balance as well as in the chemical evolution of the con- 
sidered environment: depending on their amount and on the global gas density, in 
fact, the possible evolutive scenario could be susceptible of marked variations on 
the estimated time needed for reaching the steady state, their presence having thus 
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also important repercussions on the final chemical composition of a given environ- 
ment. In contrast, the formation of positive ions requires an energy amount equal 
to or greater than the first ionization potential of the neutral molecular species 
considered, i.e. an energy around 8-10 eV for most of organic molecules, a value 
which is considerably high for several astronomical contexts, like cold dense in- 
terstellar molecular clouds and, with less importance, within diffuse molecular 
gaseous regions. 

Low-energy electrons, usually considered in the 0-10 eV range, can be very reactive 
in the sense that they are effectively captured by many molecules, which can then 
undergo either rapid stabilization toward the undissociated species or fast uni- 
molecular decomposition toward specific fragments, where, in the latter case, the 
extra electron initially captured by the neutral molecular target, can in turn remain 
stuck to one of the fragments as well as be re-emitted again in the surroundings 
with a different (usually lower) kinetic energy. Consequently, it can play a role in 
heating the gaseous medium at a molecular level. Furthermore, low-energy elec- 
trons can also be 'simply' deflected and then causing vibro-rotational excitation of 
the target molecules, therefore, playing once more a role in the heating processes 
at a molecular level. On a much larger scale, they could even produce conduc- 
tivity inhomogeneity, which in turn disturbs the radiowave propagation in such a 
medium, like on the edge or on the tail of an interstellar shock, or in principle 
also moving across a dense cloud: from this point of view, electrons can be also 
seen as an important and versatile source for 'decoupling', at least in the low fre- 
quency range, the energy transfer from the interstellar radiative field to molecules 
and nanoparticles. Hence, in a qualitative sense, all of these nanoscopic mecha- 
nisms could play a crucial role in understanding more deeply the local interactions 
within different 'phases' of the ISM, ranging from the cold dense molecular gas 
up to the diffuse molecular component. In fact, if we look at the stars, with high- 
energy thermonuclear burning reactions in their interiors, as the primary sources 
of almost all the energy that is released in the seemingly empty interstellar space. 



then we can look at the free electrons (Q) as a flexible 'means' which, in compe- 
tition with photons, is able to participate into the complex processes responsible 
for the coupling between the ionized, the dust and the neutral components of the 
ISM. 

Generally speaking, the main reasons that originally motivated us to undertake the 
present work, were at least two. 

First of all, we intended to demonstrate the importance of resonances in forming 
molecular anions in different astrophysical environments. 

Secondly, we were attracted by the possibility of investigating with a reliable 
and suitably approximate approach, the occurrence of radiationless paths hke in- 
tramolecular vibrational redistributions (IVRs) to account for the dissipation of 
the extra energy initially carried by the impinging electron. 

The former aim, which is of course a common feature to each of the molecules in- 
vestigated here, can be easily justified on the basis of the most important attempts 
recently dedicated to the understanding of anion formation processes in several 
astrophysical environments, where indeed only s-wave electron attachment pro- 
cesses were expressly considered. In the framework of the electron-molecule 
collisions, once the physical conditions like the mean free electron density and 
the mean kinetic electron temperature are defined according to the specific astro- 
nomical context under investigation, taking in consideration only s-wave electron 
attachment processes to account for the formation of molecular anions provides, 
in fact, a reliable but narrow point of view. In the qualitative sense, resonant 
electron attachment can be viewed as the doorway for the ensuing possible for- 
mation of the stable molecular anion, provided the neutral molecule has a positive 
electron affinity. It therefore follows that the resonant contributions associated to 
partial waves referred to non-spherical angular momentum states, can account for 
the formation of molecular anions also for non- vanishing collision energies which 
can be still relevant in the given environment. In this framework, once one has de- 
termined which are the astrophysically relevant resonances for a given molecule, 

4n the ISM essentially produced either by high-energy UV photons or by stellar winds 
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it becomes of interest to also investigate, for each of them, which is the possible 
evolution of that resonant species: the radiative stabilization for the extra energy 
content that characterizes the resonant species, in fact, is usually a slow process, 
that can even exceeds the resonance's lifetime, so that in such a case the autode- 
tachment becomes strongly competitive. On the other hand, the radiationless IVR 
constitutes in general another reliable evolutive path that, conversely to the sponta- 
neous emission of a high-energy(|^ photon, can efficiently lead to the stabilization 
of the metastable anion; at the same time, according to the non-ergodic vibra- 
tional redistribution, which for a pletora of molecules is an actual possibility, such 
a process can also account for specific fragmentation channels, as will be shown 
in details in the last chapter, where all the findings shall be analised and discussed. 

The present PhD thesis is focusing on electron-molecule interactions with astro- 
physically relevant molecules, and therefore represents a theoretical/computational 
work which deals with an area placed at the boundary between (molecular) astro- 
physics, quantum collision thery, and of course theoretical chemistry. The three 
molecular species that will be computationally investigated for their behaviour un- 
der low-energy electron collisions are the ortho-benzyne (0-C6H4), the coronene 
(C24H12), and the carbon nitride (NC2N), respectively. Due to their sizes, their pe- 
culiar structures, their chemical reactivity and their physical and chemical proper- 
ties as well as according to astronomical observations, when available, the above 
three molecules provide representative examples for different astrophysical con- 
texts. However, each of them is linked to a specific astrophysical aspect that 
currently constitutes an intriguing scientific challenge, hence the present study. 

The 0-C6H4 is surmised to play a role in larger PAHs astrohysical synthesis; de- 
spite its peculiar physical and chemical properties like its large dipole moment, 
the neutral parent molecule has not been yet observed in spatial regions where it 
could appear as a reasonable molecular product. By admitting the presence of the 
neutral in some astrophysical environments, like the proto-planetary atmosphere 
CRL-618, the ensuing reaction of the negative undissociated molecular species 



^of several eV: in the 0-5 eV according to most of astrophysical environmental conditions 
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with cations present in the same context (toward more complex structures), could 
be considered as one of the possible paths according which it is possible to justify 
its disappearence while, at the same time, also providing one further step within 
the mechanism responsible for the interstellar PAHs formation. Therefore, on the 
basis of the above considerations, we have investigated its feasibility in forming 
negative ions under low-energy electron collisions. 

The C24H12 constitutes another aromatic species, whose size is markedly larger 
than the previous aromatic molecule. Its astrophysical importance rests first of 
all on its (alleged) almost ubiquitous presence in several astronomical contexts, 
ranging from the cold dense molecular clouds up to the diffuse interstellar re- 
gions. Therefore, it seems reasonable to assume this molecule as one of the most 
representative interstellar medium-size PAHs([^, at least in the sense of one pos- 
sible 'founder' of a family of larger PAHs. Furthermore, according to most of 
astronomers and astrophysicists, it is assumed that 'large' Q PAHs should almost 
soak up free electrons, if present in the same spatial region, having thus deep 
repercussion on the chemical evolution of such an astrophysical environment. In 
this framework, due to its positive electron affinity, the coronene molecule can be 
viewed as one of the key interstellar polycyclic aromatic species so that, by com- 
puting the integral, the differential and the momentum transfer cross sections for 
the electron scattering process within a suitably large and astrophysically mean- 
ingful energy range, and with an ab initio approach based on a rigorous quan- 
tum dynamical treatment, it was possible to investigate its capability of forming 
(meta)stable anions under single collision conditions as well as its role in deflect- 
ing the colliding electrons. Therefore, on the basis of the above findings, it was 
also possible to characterize its global efficiency in producing the transient nega- 
tive species, as reflected by the estimate for the metastable anion formation's rate 
coefficient. As will be fully discussed in the last chapter of the present thesis, it 
is not surprising that this molecule interacts with free impinging electrons leading 

^it is one of the largest and most symmetric between the medium-size members of interstellar 

PAHs family 

"^with about more than 25/30 carbon atoms 
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to a relatively large number of resonances; for the moment it is important to point 
out that its large static polarizabiUty coupled with the large number of tt electrons 
characterizing the neutral enables one to qualitatively expect such a behaviour, 
although the response of this molecule to colliding electrons for vanishing ener- 
gies as well as its relative 'transparency' to incident electronic projectiles in the 
sense of a markedly small deflection of the colliding particle, are not immedi- 
ately obvious. Moreover, from the 'pure' theoretical point of view, the C24H12 
behaviour under low-energy single electron collisions was never computationally 
investigated before; in this connection, such a theoretical work enables to shed 
light on some of the possible mechanisms responsible for the anion formation. 
Due to its molecular size, it also represents a demanding challenge to relate the 
present theoretical findings with the available experiments, especially when one 
takes in consideration the lowest collision energy range (£■ ~ eV) under the light 
of the astrophysical role played by the C24H^2 anion. 

Last but not least, the NC2N. This molecule was choosen as an interesting candi- 
date that, in principle, can participate to the formation of one of the most abundant 
anions observed in the Titan's upper atmosphere, the CN~ anionic species, as well 
as for the production of the undissociated anion, NC2N~. The former would be 
formed if an efficient dissociative electron attachment does take place, while the 
latter formation involves the stabilization of the whole molecular structure with- 
out fragmentation paths. It therefore follows that, once the metastable species is 
formed at the equilibrium under the observed conditions of Titan's atmosphere, 
the investigation of the nuclear deformation paths is a foundamental step to es- 
tablish which anion is mainly produced as well as its formation probability. More 
in general, such an investigation could allow to find suitable environmental con- 
ditions that favour the formation of a specific product, thus providing useful in- 
formation that can be used to draw different evolutive chemical scenarios. In 
this connection, I find important to stress that the available developed chemical 
models of Titan's atmosphere, in fact, consider the dissociative path as dominant, 
according which the NC2N should mainly provide the CN~ species. 



Chapter 1 

Electron driven processes in 
astrophysical molecular gases 

1.1 Molecular negative ions, temporary molecular 
ions (TNI) and their stability: the electron affin- 
ity (EA) 

It is well known and experimentally demonstrated that many atoms and molecules 
can also exist as negatively charged species in the gas phase. Such ions can in fact 
be produced by electron transfer processes from neutral species (indirect ioniza- 
tion), like in the following examples 

N + M — >N+ + M- (1.1) 

+ M — >N + M- (1.2) 
or simply by a direct collision with an incident electron (electron attachment) 

e- + M — >M- (1.3) 

In the present section, I'm going to focus on the last case, the so called free 
electron attachment process and the eventually subsequent decomposition of the 
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molecular anion once it is formed as in equation 1.3. Of course, the previous 
two reactions refer to other possible channels by means of which it is possible 
to produce the gaseous anionc species M" ; it is however important to emphasize 
that such channels implicitly consider the presence of another negatively charged 
species in the same environment, so that can be viwed as secondary channels in 
producing complex anions. As a qualitative rule of thumb, moreover, since the 
binding energy of the extra electron to M is usually less when compared with the 
ionization energy of the other species entering the reactive collision, the first two 
reactions are often endothermic so that they can only occur if and M contain 
sufficient energy (kinetic and/or internal). An important exception to this asser- 
tion is provided by Rydberg negative ions, where the extra electron is contained in 
a very diffuse orbital having thus a very small binding energy: in this connection, 
the Rydberg anions play an important role in allowing electron transfer reactions 
close to the zero energy value. Before treating the electron attachment problem 
more in detail, it is necessary to briefly recall the main general aspects concerning 
the stability of negative ions, thus introducing the electron affinity (EA), a founda- 
mental parameter which is necessary in order to understand the reaction kinetics 
of all the above processes. 

The possibility of an atom to form a thermodynamic ally stable anion is expressed 
by the electron affinity; by analogy with the previous definition, the same property 
for a molecule M is expressed by the molecular electron affinity (EA) which im- 
plicitly refers to the parent undissociated thermodynamically stable molecular 
anion, M~. Formally the molecular electron affinity defines the energy difference 
between the neutral species M and the associated undissociated anion M~ when 
both considered in their respective ground states. By the above definition, it is 
possible to relate the neutral species and the anionic but only when each of them 
is in its own most stable geometry; this is to say that for molecules for which the 
extra electron is capable to induce a considerable geometry change between the 
neutral and the associated anion, it is necessary therefore to make a distinction be- 
tween the previously defined EA (which is better identified as adiabatic electron 
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Figure 1.1: Hypothetical Born-Oppenheier potential energy curves for the neutral 
molecule M and the associated anion M~ illustrating the difference between the adia- 
batic electron affinity (AEA) and the vertical detachment energy (VDE). See main text 
for details. 

affinity, AEA) and the vertical detachment energy (VDE), as indicated in figure 1. 
The importance of the VDE, when experimentally available, should be found in 
that such a number implicitly refers to the Franck-Condon transitions which take 
place in the photodetachment process from molecular anion M 

M- + hv — >M + e- (1.4) 

so that, in the qualitative sense, the phisical difference between the AEA and the 
VDE is reminescent of the same concepts in the case of adiabatic and vertical 
ionization energies as occur in a photoionization process. By convention, the 
electron affinity of M is taken as positive when the ground state of M~ lies below 
that of M, so that it refers to a stable molecular anion when compared to the parent 
neutral, and obviously it is considered negative in the opposite case, that is to say 
if the ground state energy of M~ is higher than M . Thus, as implicitly indicated 
by figure 1, a positive EA value for a molecule M means that the corresponding 
undissociated anion can exist as thermodynamically stable and having the extra 
electron in a bound state. 

Many molecules, when in the gas-phase undergo a direct collision with an incident 
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free electron as in equation 1.3, are not indeed able to produce a thermodynam- 
ically stable negatively charged species; it is however a concrete possibility (the 
most frequent) that the impinging electron remain temporarily 'bound' to the tar- 
get molecule: in such a case the complex between the electron and the molecule 
is indicated as transient negative ion (TNI) due to its metastable character in the 
sense of a limited lifetime. The crucial point that must be emphasized at this point 
is that the primary electron attachment is formally a molecular electronic transi- 
tion which in any case, therefore regardless the electron affinity sign, initially does 
create a negative species unstable with respect the ejection of the extra electron 
itself. So, if the EA is positive, the electron attachment leads to an electronically 
excited state of the anion which can have time enough to give rise to a transition 
to the thermodynamically stable ground state: the last step depends in fact on how 
efficient the subsequent relaxation mechanisms are, since the kinetic energy car- 
ried by the colliding electron can be dissipated either into the molecular degrees of 
freedom or radiatively by photon emission provided the lifetime is long enough. 

When considering a species with negative EA, on the contrary, the ground state 
negative ion is unstable with respect to the extra electron autodetachment, so that 
the stable negative ion cannot be produced. 

Figure 2 pictorically summarizes the differences between a species having a posi- 
tive and one having a negative EA, illustrating the situation for a neutral molecule 
M, its cation M"*" and the corresponding anion M~ in terms of the total energy 
(upper half), as well as of the binding energy (lower half) of the electrons. As 
such figure clearly shows, regardless the EA sign, the neutral and the cation are 
characterized by the presence of a large number of virtual (unoccupied) MO lo- 
cated below the dissociative energy limit (hence associated to bound states) and 
that thicken as the energy approaches the corresponding ionization limit. Conse- 
quently, both the neutral and the cation are in principle able to capture an inci- 
dent electron in a MO above the corresponding LUMO but below the dissociation 
limit. In contrast to the above observations, the same figure shows a situation 
completely different for the anion M . The latter is in fact characterized mainly 
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Figure 1.2: Total energy and molecular orbitals schematic representation for a neutral 
molecule M, the associated cation M'^ and anion M~ when considering both a positive 
and a negative electron affinity. 
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by a less value for the binding energy of the extra electron, as well as by a reduced 
number of possible bound states, if they exist: in this connection, when the EA is 
negative, the figure pictorically shows in fact the extra electron within an excited 
MO which is not bound since it is located above the ionization limit. Qualita- 
tively speaking, this feature is susceptible of the following explanation: the extra 
electron bound in M~ can be considered as interacting with a neutral 'core' (i.e. 
M) and therefore, increasing its distance from the neutral molecule, it 'feels' an 
induced-dipole potential like 

nr) = (1.5) 

a being the polarizability of the neutral M. An analogous (and qualitative!) ap- 
proach involving both the neutral M and the cation M+ therefore leads us to 
recognize that the extra electron should be viewd as 'bound' to a single and a dou- 
ble positively charged core, respectively: in these cases the potential 'felt' by the 
electron is clearly Coulombic 

y(r) = -- (1.6) 
r 

where it is sufficient to think about the hydrogen atom problem to immediately 
justify the larger number of possible bound states as well as their becoming thicker 
as the energy approaches the ionization limit. 

1.2 Collisions of electrons with molecules. TNI for- 
mation and resonant states 

The present section focuses on the interaction of a free electronic particle with a 
gas-phase molecule under the so called under collision conditions: we are thus 
considering here conditions which require a mean free path lenght large enough 
to ensure that secondary collisions (either with other atomic/molecular species or 
with other electrons) have very low probability to take place. The above con- 
straint, as we shall see also in the next chapter, turns out to be reasonably satis- 
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fled when referring to different astrophysical environments within the interstellar 
medium (ISM) where the density of matter is low enough to allow its application 
as a realistic approximation. In this framework, the interaction of a free incident 
electron on a molecular target can be divided into two classes, namely direct 
scattering and resonant scattering. As the name itself suggests, the former refers 
to the situation where the impinging electron with kinetic energy +/i^A;^/2me 
collides with the target molecule and will eventually deflected from its original 
trajectory, as depicted in the upper part of figure 3. 




Figure 1.3: Pictorical representation of direct, resonant and dissociative scattering. 

Obviously, if the energy of the colliding particle results unaffected by the im- 
pact (the moduli of the initial and final wave vectors are equals), the scattering 
is also elastic. Furthermore, the colliding electron can instead loose energy to 
some extent, which therefore results in the excitation of the internal degrees of 
freedom of the target (or in its electronic excitation) so that the scattering event 
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is inelastic. Keeping in mind that the mass ratio between the electronic particle 
and the lightest nucleus is a very small number {nie/mH ~ 1/1836), the direct 
excitation of the roto-vibrational state of the molecular target is not prevented, 
but strongly unUkely; hence the direct inelastic collisions are mainly responsible 
for the electronic excitation of the target which, in turn, could evolve in yielding 
an excited roto-vibrational state of the ground electronic state by internal con- 
version and, in general, depending on how the strength of the coupling between 
the metastable electronic state and the nuclear dynamics is. On the other side, 
the resonant collision does occur when the incident electron remains 'partially' 
bound on the molecular target for a time longer than that needed to only cross the 
region of space occupied by the molecule itself, where partially should be then 
red in the sense of temporarily. At this point it is appropriate to point out that an 
electron capture process can occur only if the kinetic energy of the colliding elec- 
tron is 'suitable', i.e. only if the incident electron kinetic energy is 'in resonance' 
with the formation energy of the transient negative ion, so that the expressions 
'transient', 'metastable' and 'resonant' are synonymous. Having thus definied a 
resonant state, it is implicit that the captured particle can decay after a finite time. 
Consequently, in the framework of a stationary representation of the scattering 
process, the time dependence of such a metastable state is reasonably expressed 
by the well known imaginary exponential function 



where the total energy E is complex, according to = E^^ - ir/2. Taking advan- 
tage of this allows one to write 



so that it is clear that the resonant state decays after a lifetime r given by 



the last espression being reminiscent of the time-energy uncertainty principle (h ~ 




(1.7) 




(1.8) 



(1.9) 
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6.6 ■ 10 eV sec ). The imaginary component T/2 of the total scattering energy 
E is the resonance halfwidth. Generally, the resonance lifetimes r, depending on 
both the coUisional energy of the electron (i.e. the modulus of its linear moment) 
and the molecular size, are susceptible to cover a large window, ranging from few 
vibrational periods (N^ : r ~ 10"^^ - lO^^^sec, E/j ~ 2.3 eV BUEl) up to much 
longer values {SFq : r ~ lO^^sec, [|3l|4l): in this connection it might be useful 
to emphasize that a free electron of about 1 eV needs of 5 ■ 10^^^ sec to cover a 
distance of 3 • 10~^° m, while a chemical bond typically has a spatial extension 
which is restricted between 1.2 ■ 10~^° m and 2.8 ■ 10^^° m. 

1.3 Resonance classification 

So far almost nothing has been said about the mechanisms by which a resonant 
state can be produced. In this regard, there are many different ways such that an 
extra electron impinging on a molecule can remain temporarily bound on it, and 
according to them it is therefore possible to make a general distinction between a 
Feshbach resonance and a shape resonance. So it is appropriate to emphasize that 
a quantitative discussion concerning a resonant scattering process is based on the 
formalism of quantum mechanics; without going into details, we should however 
keep in mind that, in the context of a quantum approach, one has to do with a wave 
packet responsible for the physical description of an electron (coUimated) beam, 
incident on a molecule, which is represented in the following form [|5l 

^!{z,t)= dk A{k) e'^' e-'^' (1.10) 
Jo 

where 'z' in turn represents the incident direction of the beam while A{k), the 
weight function for the plane waves e*^^ continuum summation expressed by the 
previous integral, is supposed to have compact support centered on a certain 
value. Therefore, in relation to the possibility of developing each plane wave in 
spherical waves 
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^ikz ^ ^iy47r(2/ + l)ji{kr)Y^{^) = ^i\2l + l)ji{kr)Pi{cos^), (1.11) 

/=0 l=Q 

where '/' labels each partial angular component of the plane wave, it is possible 
to qualitatively understand the basis mechanism of a resonant scattering process: 
such a process is always describable as a constructive interference between the 
target molecule and one or more partial waves associated to the wavepacket which 
properly represents the dynamical evolution of the colUding beam. Having so 
established a plausible definition, which has a general nature resulting effective 
as long as it is understood as qualitative and heuristic, let me now add somtething 
little more specific from the theoretical point of view, always remaining in the 
scope of a general description. 

1.3.1 Feshbach resonances 

It is now essential to consider the electronic configuration of a generic tempo- 
rary negative molecular ion. If the extra electron is able to remain trapped on the 
molecular target, so producing a resonant state, but without introducing any varia- 
tion into the electronic configuration of the target molecule itself, we are speaking 
about a 'one' particle resonance: such a process can therefore be imagined as that 
one where the TNI is ideally obtained simply adding an electron in one of the 
virtual (thus available) molecular orbitals of the neutral species involved in the 
collision. On the other side, if the extra electron does cause a variation of the 
target electronic configuration while remaining temporarily bound to it, we would 
be dealing with a so-called 'two-particle' resonance, also known as 'core-excited' 
resonance. In the latter case, the energy of the metastable species (generally indi- 
cated as M*^) can be located either above or below the energy of the correspond- 
ing electronically excited (but not negatively charged) species M*. Depending 
on whether the energy of the excited species M*~ is greater than those of M* 
or not, the negative electronically excited state is an open or a closed-channel 
core-excited resonance, respectively; hence, the closed-channel resonance, also 
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identified as electronic Feshbach resonance, cannot directly decay into the ground 
electronic state. As figure 4 clearly shows, such a possibility is necessarily ac- 
companied by a modification in the electronic configuration, so that electronic 
Feshbach resonances have usually long lifetimes. 



Energy 

A 



M 



M" 



Open Channel 

Closed Channel 
(Feshbach) 



Shape 



Neutral 



NI Resonance 



Figure 1.4: Representation of M, M*, and M*^ (as displayed by increasing en- 
ergies) to schematically illustrate the differences between a core-excited open channel 
resonance and a Feshbach electronic resonance. The figure also make possible a qualita- 
tive comparison with a shape resonance. 'NF = Negative Ion. See text for details. 

However, always in the qualitative sense, the colliding electron can also induce 
a coupling with the nuclear dynamics which in turn is sufficiently strong so that 
the energy transfer from the electronic to the nuclear degrees of freedom can actu- 
ally participate in dissipating its excess (kinetic) energy amount, thus preventing 
the ejection of the extra electron and eventually leading to a nuclear Feshbach 
resonance. 

At this point, what else about the general mechanism responsible of the electron 
'capture' in a Feshbach state? When the electron projectile induces an electronic 
excitation into the target molecule, providing an electronic Feshbach resonance, 
the extra electron is simply detained within the molecular potential field generated 
by the excited molecule: having loose too much energy to cause the electronic 
excitation, in fact, it cannot escape since now it has no sufficient energy available 
so that the target remains in the excited state; hence, before the electron can be 
emitted, it must reabsorb energy from the excited target. It is of course possible to 
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assume a little more technical point of view, without any pretense of completeness. 
Let us therefore immagine the total resonant wavefunction \l'n(r, R) developed 
in the following form, like in a close-coupling expansion 

vl>„(r,g,R) = 5^0,(g,R)x.(R)/n,(r) (1.12) 

i 

where r, g, R indicate the electronic projectile spatial coordinates, the ensamble 
of molecular electronic coordinates (both spatial and spin) and the nuclear coor- 
dinates, respectively. The symbols 0j, Xi ^^id /„ j refer then to the i-th electronic 
eigenstate for the molecule kept with fixed nuclei, the i-th nuclear state, and the 
wave function for the incident electron. 

Based on the above premise, it is then possible to say that a 'pure' electronic Fes- 
hbach resonance does occur when the impinging electron excites one(Q) electronic 
state, for example the j-th, for which the vibrational levels are all located above 
the resonance energy, so that the summation reported in eq. 12 should collapse in 
one term, the j-th. Consequently, is fairly straightforward to introduce a nuclear 
Feshbach resonance in a similar way, but now considering, as final achievement, 
a vibro-rotational excited level of a certain molecular electronic state which, as 
final result of the collision, can or cannot remain that one initially possessed by 
the neutral molecular target. 
A conclusive observation. 

If, when treating electronic Feshbach resonances, it is reasonably allowed (as 
moreover has been implicitly done) to consider a molecular geometry with fixed 
nuclei, using thus a Bom-Oppenheimer description of the molecular electronic 
levels, moving on nuclear Feshbach resonances such an assumption is no longer 
valid. For a nuclear Feshbach resonant state, in fact, the energy exchange due to 
the coupling between the electronic and the nuclear degrees of freedom comes to 
play a crucial role, thus excluding a Born-Oppenheimer dynamical description: it 
is just the fixed nuclei approach that prevents such a 'dialog' between the 'exter- 
nal' (electronic) and the 'internal' (nuclear) degrees of freedom. In other words, if 

'for the sake of simplicity, we are considering here the excitation of only one electronic state 
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the molecular nuclei were considered kept as fixed in a certain geometry configu- 
ration, the resonant (metastable) nuclear Feshbach negative state would be stable 
against the electron ejection, so that the nuclear Feshbach state would therefore 
be meaningless. 



1.3.2 Shape resonances 

Another mechanism, which is applied to physically explain both one-particle reso- 
nances and open channel core-excited resonances, describes the electron trapping 
by means of the shape of the effective interaction potential between the incident 
particle itself and the molecule. The electronic projectile, in fact, remains trapped 
by means of, or better, within the interaction potential because of its shape: ac- 
cording to this statement, in fact, one speaks about shape resonances. To get a 
qualitative idea of that, it is sufficient to assume that the colliding electron expe- 
riences an attractive potential in a spatial region centered in the molecular center 
of mass which is in turn surrounded by a region of repulsive potential whose ex- 
tension is spatially limited, thus vanishing for large distances. Depending on the 
electron kinetic energy, the electron can tunnel the external barrier; in such a case, 
then, the combined action of the latter (which is now an obstacle for the electron's 
escape) with the internal attactive region can prevent the electron re-ejection (au- 
todetachment), at least for the time necessary for the negative transient species to 
be experimentally observed (for example, by electron energy loss spectroscopy or 
electron transmission spectroscopy). At this point is legitimate, even necessary, 
to look for the origin of such a potential: its shape is in fact responsible of the dy- 
namical trapping of the incident electron. For the sake of symplicity, we can thus 
suppose the target molecule to have a spherical polarizability, so that we come to 
the following approximated expression for the effective interaction potential 

T. / N l{l + l)ti' 

K//(r) = -^ + A_^, (U3) 

where the first term was already introduced (thus not accidentally; eq. 1.5) in 
order to qualitatively describe the attractive interaction between the external elec- 
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tron, in a given l-th state of angular momentumi, and the 'neutral, soft' internal 
core in a molecular anion. Thus, following the present heuristic model, it is just 
the combination of the negative (attractive) dipole-induced potential with the pos- 
itive (repulsive) branch, associated to the angular momentum of the electron, that 
globally results in producing the potential repulsive barrier. 




r 



Figure 1 .5: Representation of the radial effective interaction potential responsible for the 
dynamical trapping of an incoming electron in a given state of angular momentum / > 0. 
e indicates the extra electron kinetic energy. 

Of course, including only the polarization term in the total interaction potential 
Veff{r), provides an understandble but not physically realistic picture by which 
one describes the capture process. As we will discuss more in detail in chapter 
2, in fact, the dipole-induced term (which moreover usually dominates at very 
low-energy) should be added with the exchange interactions as well as with an- 
other term which accounts for the electron correlation; however, for the sake of 
simplicity, it might be sufficient here to take in consideration only the indicated 
term. More in detail, as it is illustrated in figure 1.5, for short distances from the 
molecular center of mass it dominates the repulsive Coulombic electron-electron 
interaction, which gives rise to the rapidly increasing branch of repulsive poten- 
tial close to the zero distance value. Covering the remaining radial distance, then 
going away from the target center of mass, that figure clearly shows the presence 
of an attractive well (just before the centrifugal barrier) which in principle allows 
the metastable state (a bound state placed at positive energy) to exist for a limited 
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time. 

It is now possible to write down the following general statement: if the target 
molecule has a virtual MO, energetically accessible and characterized in terms of 
its spatial symmetry by a suitable value of electronic angular momentum I, then 
an incoming electron that is in turn characterized by the same angular component 
I can be temporarily captured by tunneling the above barrier, being thus located 
within the attractive potential region; in other words, since the above point of view 
is strongly reminiscent of wave mechanics, we can also say that it could take place 
a constructive interference involving the l-th partial wave of the spherical waves 
summation (eq. 1.11) describing the colliding electron and the above suitable vir- 
tual MO. Of course, the energy of such a virtual MO, generally speaking, do not 
exactly matches the resonance energy, i.e. the electron kinetic energy for which 
the metastable anion is formed. The reason should be found in the fact that even- 
though such a description seems to be physically reasonable, on the other hand it 
is not considering the electron as a dynamical perturbation to the molecular poten- 
tial field, so that such a virtual MO could or also could not exist before the capture 
process itself takes place; in this sense, the metastable state should be viewd as the 
electronic dynamical molecular response to the impact of a free electron with 
a suitable kinetic energy: hence the energy spectra of the virtual MOs appears to 
be modified, the corresponding energies being consequently shifted. 

Moreover, the resulting TNI constitutes an example of discrete state embedded 
in the continuum: according to this, such a species is also called quasi-bound 
state. Obviously, resonant states like these cannot exist in principle for an s- 
wave attachment process, unless the effective interaction potential deprived from 
the centrifugal term is however able to give rise to a similar behaviour, i.e. in 
producing some barrier behind which an attractive well may exist. 

As figure 1.6 pictorically shows, the radial scattering wavefunction behaviour is 
susceptible of substantial variations depending on the colliding kinetic energy: 
only when the collision energy matches the resonance value (e = e^, central inset 
of figure 6), the incident particle is represented by a wavefunction which has a not 
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Figure 1 .6: Behaviour of tlie radial scattering wavefunction as the kinetic energy crosses 
the resonant value (e = e^ ). It is implicitly assumed that the interaction potential is not 
affected by little energy variations around the resonant value. See text for details. 

negligible squared modulus inside the centrifugal barrier region, then justifying 
its dynamical trapping and at the same time allowing to the extra electron to be 
released by the same tunneling mechanism, but in the reverse direction. When 
such a condition is not satisfied (upper and lower insets in the same figure), the 
scattering wavefunction get to be 'simply' reflected by the angular barrier: in 
this sense the formation of a TNI provides a growth of probability density of 
the extra electron just within the spatial region corresponding to the attractive 
well. Depending of the attractive potential branch deepness, the spatial region 
where the dynamical trapping does occur can be shifted closer the molecule or 
not, resulting then in a metastable species where the extra electron is less or more 
diffuse, respectively. As we will see in the case of ortho-benzyne metastable anion 
(chapter 3), this feature qualitatively provides useful information about the anion 
stability in the sense that the less diffuse the resonant charge is, the more stable 
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the TNI: in such a case, in fact, the subsequent dynamical evolution toward the 
stable anion by radiationless processes does appear to be in general helped. 
To bring to a conclusion the present subsection, let us now consider the case of a 
finite-ranged electron-molecule interaction potential, that we can generally write 
as 



Veffir) = Vir) + ^-^L±}K_ (1.14) 



In the realistic case of a finite barrier (i.e. not too much close to the molecular 
target), the electron can be trapped inside the attractive well, but as we have pre- 
viously discussed it cannot be trapped forever. From the pure theoretical point 
of view, a rigorous non-relativistic quantum mechanical description demostrates 
that to each scattered l-th partial wave there corresponds a phaseshift the 
behaviour of which is indicative of the possible resonant capture occurrence: for 
each partial wave involved in a resonance, in fact, as the collision energy crosses 
the quasi bound metastable state value, the associated phase shift rises through 
the Tc/2 rad value up to about vr. At the same time, the corresponding partial 
wave cross sections a/, which in general provide the transition probabilities for 
the TNI formation per unit time, per unit target scatterer and per unit relative flux 
of the incident electrons with respect to the molecular target, pass though their 
maximum values given by |f (2/ + 1), so that when a shape resonance dominates 
the l-th partial- wave cross section around the collision energy Er, the well known 
Breit-Wigner resonance formula [|56l does apply(|^ 

where T refers to the full width at half maximum which, according to the energy- 
time uncertainty as expressed in eq. 1.9, provides an estimate for the resonance 
lifetime. 



^which is a function of the electron collision energy or, equivalently, of the associated linear 
momentum k 

^provided the resonance is reasonably narrow so that variation in can be ignored 
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1.3.3 Low-energy scattering and bound states 

At very low energies, partial waves for higher / are, in general, unimportant. This 
point may be obvious classically because the particle cannot penetrate the cen- 
trifugal barrier; as a result the potential inside the attractive well has no effect. 
As previously introduced, the effective potential for the /-th partial wave is given 
by eq. 14 so that, in terms of quantum mechanics, unless the potential attractive 
region is strong enough to accomodate I ^ quasi-bound or bound states near 
E = eV, the behaviour of the radial wavefunction is largely determined by the 
centrifugal barrier term, which means that it must resemble the free radial wave- 
function ji{kr). A little more quantitatively, in the framework of the partial waves 
method, it is possible to estimate the behaviour of the phase shift Si using the 
integral equation for the corresponding /-th partial wave, obtaining ^ 



where fi{k,r) under integration on the r.h.s. is the scattering radial wavefunction. 
If fi{k, r) is not too much different from ji{kr), and l/Zc is much larger than the 
potential range, it can be shown that the right hand side would vary approximately 
as /c^' if one now supposes that 6i is small, so that we are implicitly dealing 
with an interaction potential not so strong, the left hand side must consequently 
vary as 6i/k. It therefore follows that, for low-energy incident electrons (k ~ 0) 
on a neutral molecule (the interaction being reasonably finite-ranged), the partial 
phase shift goes to zero roughly as 5i ~ this means that, even excluding the 

partial wave contribution corresponding to / > 0, the s-wave (Z = 0) scattering can 
be important. 

Zero energy scattering: s-wave capture, bound and virtual states 

To gain a little deeper insight in the last statement, we can start from the radial 

equations describing the scattering as provided by the formalism of partial waves: 

^according to one of the properties of the spherical Bessel functions, for very low incident 
energies (so that fc ^ 0) it results j/(fcr) w rii+'i')!! 
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ue{k,r) + 



i2 2m 



ui{k,r, 



0. 



(1.17) 



Focusing on the s-wave case, for very low-energies one obtains 



rf2 



■Uo{k,r) 



2m 



V{r)uo{k,r) = 0. 



(1.18) 



In the region where r » R^, Ry representing the range of the potential interaction, 
the previous equation becomes 



Uo{k,r) = 



(1.19) 



therefore, literally admitting as solution a straight line 



Uo{k,r) = c ■ {r - a) , r » R^ (1-20) 

where 

a = lim (1.21) 

fc-5.o+ k 

defines the scattering length Q. In what follows, the crucial point rests on the 
simple geometrical meaning for a: the scattering length, in fact, is the intersection 
of the asymptote of Uo{k,r) with the r axis. Neglecting for simplicity sake the 
case of a repulsive potential, it follows that the scattering length can assume either 
positive or negative values so that, depending on the attractive potential's depth 
(see figure [L7| ), the sign of a is indicative of the presence(a > 0)/absence(a < 0) 
of a bound state. 

Thus, once again, the dipole induced interaction (hence, the neutral molecular po- 
larizability) plays a foundamental role in the sense that for vanishing energies it 
can yield an s-wave bound state: in other words, an s-wave attachment not nec- 
essarily produces a stable negative ion. Keeping in mind that the conditio sine 
qua non in order to have a stable negative ion is represented by a positive elec- 
tron affinity value, the above simple analysis also suggests that depending on the 
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Figure L7: Behaviour of the scattering length for various attractive potentials. From 
left to right, as the interaction potential becomes more and more attractive, the panels 
illustrate three cases: (a) a < for which a bound s-wave state cannot exist, (b) a = oo 
for which the potential is almost able to support a real bound s-wave state (virtual state) 
and (c) a > correponding to an attractive interaction strong enough to support a bound 
s-wave state. 

molecular polarizability, a molecule can provide not efficient s-wave attachment 
and, at the same time, efficient higher-order partial waves resonances that occur 
for energies higher than eV but still astrophysically relevant. 



As further shown by figure 1.7 it can also happen that a = oo: this peculiar case 
corresponds, as we will discuss in detail about the coronene molecule, to a zero 
energy resonance In such a case, on the basis of the geometrical meaning of a, 
the interaction potential is nearly strong enough to support an s-wave bound state 
so that, depending on the nuclear motion, such a virtual state can subsequently 
evolve toward a bound state or not. 
An observation to conclude the present section. 

In the framework of the partial waves method, under the hypothesis according 
which 5o{k) ~ when k -^■ 0^, it is possible to write down the following 
relations for the i =0 partial amplitude Q: 



lim Ti=o = lim -e^-^^^^^sin^ofA;) = - sm6o(k) ~ lal, (1.22) 
fc-5.o+ fc-s>o+ k k 



On the basis of the relation 
^also called virtual state 
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da 

dn 



(1.23) 



which represents the link between theory and experiments, it follows that 



Therefore, for a neutral molecular target without permanent dipole, a very large 
value of the integral cross section at very low energies, is indicative of the possible 
presence of a zero energy resonance. 

1.4 Negative ions in the interstellar medium 
1.4.1 General overview 

With the expression 'molecular plasma', one implicitly refers to a (partially) ion- 
ized gaseous medium which contains rather equal amounts of positively and neg- 
atively charged atomic and molecular species as well as neutral molecules. Nowa- 
days we have learned that matter in Universe is mostly in the state of plasma, so 
that it is not so surprising to encounter also different kinds of molecular plasmas 
in the huge seemingly empty interstellar space of the known cosmos. Between the 
typical examples of a molecular plasma, without going too much away from our 
little planet, one is provided by the ionosphere surrounding the Earth and some 
other planets/satellites (like Titan's upper atmosphere [Si): in this portion of the 
atmosphere, both atoms and molecules are ionized mainly by the high-energy UV 
and X-ray photons which, in turn, are produced either from the Sun or from other 
intense radiation sources located outside the Solar System. Another example of 
a region of the Universe which can be globally considered as a molecular plasma 
is then represented by a stellar atmosphere, which little more in detail consti- 
tutes a kind of high-temperature molecular plasma (thus usually characterized by 
a greater degree of ionization). An interstellar cloud and a circumstellar enve- 
lope are regions of the interstellar medium (ISM) characterized by rather dense 




(1.24) 
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gaseous and solid matter and which usually contain a variety of different atoms 
and molecules. Since the interstellar/circumstellar clouds often contain electrons 
and ions (ranging in different amounts and depending on the environmental astro- 
physical conditions), they also could be viewd as examples of immense molecular 
plasmas. Thus, more in general, the whole interstellar space, seemingly almost 
empty, is filled with low and mainly very-low density matter, sometimes very ten- 
uous and tipically with non-uniform mean local density, that often is in the form 
of a molecular plasma due to the ionization caused by the interstellar radiative 
field and by the cosmic rays (very high-energy electrons and protons). 

For tens of years it was a common thought that negative ions and especially molec- 
ular anions were species which did not deserve so much attention due to their dif- 
ficulty in surviving in an environment rich of ionizing radiation. In other words, 
even though the cosmic -ray penetration as well as the ionization caused by the 
presence of UV photons, the latter produced by vicinal stars and acting at the 
outer edges of a given interstellar gaseous cloud, do produce an initial assortment 
of positive ions and free electrons, it had been concluded however that the subse- 
quent formation processes of molecular negative ions through electron attachment 
were scarce, so that they should only constitute a negligible contribution to the 
ion-molecule interstellar synthetic processes. 

Historically, however, the possible existence and role of negative ions in the inter- 
stellar medium had been the object of a few intriguing investigations PilTOl. 

The first who really suggested and revealed in advance that 'large' molecular neg- 
ative ions could actually exist in dense interstellar clouds was E. Herbst [fTTI . 
which surmised that the large electron affinity of carbon chain molecules and hy- 
drocarbon radicals with more than 4-5 atoms would have been responsible to yield 
rather high radiative electron attachment rate coefficients, thus resulting in anion- 
to-neutral ratios on the order of a few percent. To be precise, in such a pioneer- 
ing work he speculated that for specific linear carbonaceous radical species like 
C4H, C3N, C5N, C7N, CgN and other oxigen-containing species [[TT]|. with elec- 
tron affinity (EA) values greater than 2 eV, the corresponding abundance ratios 
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[M~]/[M] could be high enough (0.01-0.1) so that the corresponding undissoci- 
ated molecular anion might be detectable in such an astrophysical environ- 
ment, provided the associated microwave spectral emission/absorption lines were 
determined by laboratory measurements. 

The confirmation of this hypothesis was achieved many years later following the 
laboratory measurements of the rotational spectrum of the CeH molecular an- 
ion tfT3ll . The above measurements, in fact, from already exsisting astronomical 
spectra [fT2l . allowed the verification of the presence of CgH" in the molecular 
envelope of the C-rich evolved star IRC-i-10216, with an abundance of 1-5% that 
of the neutral parent species. The ensuing successfuU astronomical search for the 
same molecular negative ion in TMC-1, the dense molecular cloud in Taurus 1, 
determined in this dense region an anion-to-neutral ratio of about 2.5 % [fT3ll . 

Nowadays, we know that six different anions do exist in the ISM: CqR , CiR , 
CgH^, CsN^, CfjN^ and CN^, each of them being a carbonaceous polyyne-like 
linear-chain structure which includes periferic atoms like H or N. 

The rotational spectra of C^R , CsH and CsN^ have since been measured in lab- 
oratory [El [IS with subsequent detections of C4H- |Il6l, CgH" O and C3N" 
[[T5l in the molecular envelope of IRC-i-10216, where CgH" was observed also in 
TMC-1 [[TSl . Prompted by the detection of several carbon-chain molecules with 
large abundances in the protostar L1527 EOll . the CgH" anion was first searched 
for and then confirmed to be present also in this source [|19l , in turn followed 
by the observation of €4!!" [|2T|. A survey of galactic molecular sources by 
Gupta and coworkers [|22l detected once again CeH^ in two further sources, the 
pre-steUar cloud L1544 and the protostar L1521F, therefore suggesting the likely 
ubiquitousness of molecular anions in the interstellar/circumstellar medium and, 
in particular, of CeH^. The next to last known anion, the CsN^, was observed by 
Cernicharo and coworkers [23], which have attributed a series of rotational lines 
observed in the envelope of IRC-i-10216 to this new anionic species. Last but not 
least, [|24l the CN~ ion was detected, with unexpectedly high abundances, in the 
molecular circumstellar envelope of IRC-i-10216. 
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On the other side, polycyclic aromatic hydrocarbon (PAH) molecules in different 
states of charge (including also anions) and hydrogenation, which are thought to 
be an important component in the ISM [i25ll26ll27ll28l since some of their spectro- 
scopic features are observed in the infrared emission bands [|29ll30ll3Tll32ll33ll34]| . 
are currently included in several astrophysical models even though no specific 
PAH responsible of the above feature has been unambiguously identified [|35l. 
Their role, generally speaking, is found when considering the charge balance, the 
chemical and the dynamical evolution (ambipolar diffusion, UV screening effect, 
general degree of ionization, presence of free electrons and their repercussions on 
both the attenuation and reflection of radio-waves) within several regions of the 
ISM. In this connection, currently there are few experiments and/or theoretical 
investigations that focuse on their behaviour under low-energy electron collisions, 
and despite the fact that some consolidated hypothesis exist that point toward their 
actual presence in the ISM Il36ll371 . together with some good agreement of a few 
astrochemical models involving PAH anions [|25l l26ll . they are not constrained 
enough to be explicitly (i.e., as individual species) considered in astrochemical 
models: beyond the unresolved problem about their sizes and abundances as func- 
tions of the global physical conditions of a given astrophysical region, no realistic 
ab-initio estimates of their electron capture rate coefficients are currently avail- 
able. 

As previously mentioned, Herbst [[TT| surmised that, subjected to specific struc- 
tural constraints and specific conditions inherent their electron affinities, spe- 
cific molecular negative ions could be efficiently synthesized, in dense interstellar 
clouds, by rapid radiative electron attachment processes 

M + e- — vM- + hv (1.25) 

where 'M' refers to the neutral molecular species. During the previous years, 
the importance of radiative association processes for the interstellar synthesis of 
molecules began to be rapidly explored with enthusiasm [|38l [39l l40l |4T| . and a 
similar approach was then followed in order to provide reasonable estimates for 
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the formation rate coefficients of molecular anions associated with appropriate 
neutral species P0ll4Tl[mi54ll55l . The first step of the original theory, as illus- 
trated with full particulars in POl . is to assume that such a process does occur in 
two distinct stages, where the free electron and the neutral molecular target 
M form a long-lived metastable complex (Af^)* which in turn either dissociates 
back into the reactants themselves or is stabilized by the emission of a photon, as 
reported by the following reaction scheme 

M + e- ^ {M-y (1.26) 
{M-y M + e- (1.27) 

(1.28) 

(M^)* -H M- + hu 

The symbols kj (cm'^/sec), k^ (sec^^) and krs (sec~^) refer to the rate coefficient 
for the formation of the metastable species (or transient negative ion, TNI), the rate 
coefficient for the backward process by means of which the TNI get destroyed by 
the ejection of the extra electron (autodetachment) and the rate coefficient asso- 
ciated with the radiative stabilization, respectively. The overall rate coefficient 
for the whole radiative association process, kuad, is thus given by the following 
relation 

[Kh + Kyg ) 

if the metastable species is assumed to be at steady state [|40l[TT]| . In the customary 
limit according which k^ >> krs, the last expression reduces to 

kf 

kRad = y-Ks] (1.30) 
kb 

in the unlikely event that k^s > h, k^ad 'saturates' at kj, so that in general the 
value of k^ad cannot exceed the rate coefficient for temporary electron capture kj 
[|4T1[TT]| . Supposing the thermal equilibrium, the capture and auto-ionization rate 
coefficients {kf and k^, respectively) must be related by the detailed balance prin- 
ciple, according which the probabilities for the forward process and for its inverse 
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are equal. Following this assumption, and assuming that no vibrational excita- 
tion on the neutral species 'M' need to be considered, it has then been possible to 
achieve to an analitical approximated expression to get the rate coefficient for the 



whole negative ion formation process, as reported in eq. 1.31 



where knad is expressed in cm^ sec~^, h is the well known Planck's constant, 
K the Boltzmann's constant, T the temperature, qm- and qm are the electronic 
degeneracy factors of the neutral and anionic species (M and M~ respectively), /ie 
refers to the reduced coUisional mass and {EA) means the vibrational density of 
states of the negative species at an internal energy equal to the electron affinity 
of the neutral. As clearly shown by the above formula, the rate of attachment 
strongly depends on the density of vibrational states of the anion and more in 



general, even if not directly deducible from eq. 1.31 in the framework of the 
phase-space theory where the angular momentum is conserved, one assumes that 
once the metastable negative ion, the TNI, is formed by electron attachment, all 
the vibrational states in all accessible electronic states BOl [55ll can be rapidly 
formed with restrictions inherent only the electronic spin, therefore intuitively 
suggesting the importance of internal conversion. Following this approach, when 
considering only s-wave electrons (^S), one easily account for the production of 
the and '^S states as the possible products associated to a neutral open-shell 
molecule ^S, the same multiplicity ripartition being valid for a neutral open- 
shell state which can thus produce ^11 and '^11. 

On the basis of direct detections (when available), or equivalently, following the 
results provided by several models where negative ions are surmised to play an 
important role in different regions of the ISM, one could say that including anions 
in an astrophysical model means to also realistically (as much as possible) ac- 
count for their initial formation: in this sense, a pure s-wave attachment coupled 
with the only radiative stabilization path to account for the ensuing dissipation 
of the extra energy content carried by the colliding electron might be little re- 
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ductive. When looking at the resonances as doorways to form negative ions, in 
fact, our well consolidated ab-initio approach enables us to take in consideration 
also the non-spherical angular contribution, therefore going beyond the s-wave 
attachment and allowing us to provide further useful information. In line with 
that, once the TNI is formed at energies that of course must agree with the envi- 
ronmental conditions for the astrophysical region under investigation, we find also 
crucial to investigate with a captivating and not so much computationally demand- 
ing approach the TNI evolution: this will be described in sec. 1.4.3. Generally 
speaking, beyond the fact that the ab-initio realistic estimate that we provide for 
the metastable anion formation rate coefficient kf{T) (see sec. 1.4.2) already con- 
stitutes a very important parameter requested by astrophysicists, our results can be 
helpful in predicting which are the likely products of the electron collision given 
a specific astrophysically relevant neutral molecule; at the same time, the com- 
parison with the available experiments enables us to gain a deeper insight in the 
nanoscopic mechanisms that are foundamental in additionally suggesting which 
are the chemical species that might be introduced in an astrochemical model. 

1.4.2 The metastable anion formation rate coefficient, kf{T) 

From a purely theoretical point of view, the rate coefficient for reactive processes 
between two colliding species is defined by (cr„(t>)t>), i.e. by the thermal aver- 
age of the product of the state dependent reaction cross section cr„ (f ), n denoting 
the internal states of the species involved in the binary collision, times the rela- 
tive collision velocity, v. The above definition can be also applied in the case of 
metastable anion formation processes by electron resonant capture as well as for 
dissociative electron attachment processes, where therefore the reactant species 
are the free incoming electron and the neutral molecular target [^4211 . 
By assuming that the electrons, having kinetic temperature Tg, and the target 
molecules described by their internal temperature T m are in thermal equilibrium 
corresponding to the temperature T = Tg = Tm, it therefore follows that the general 
expression for the metastable anion formation rate coefficient kj is given by 
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1 v-^ _I!Vl 

kf{T) = Q^J2^~"'"Mv)v) (1.32) 



1 ""fj f 

Y^e'T^ V ar,{v) f{v, T) dv 



QvQ 

where c'"„(f ), as a function of the relative electron velocity, is the integral cross 
section (ICS) for the electron-molecule collision for a given initial vibrorotational 
state n, e„ is the vibrorotational energy of the neutral target molecule (from the 
bottom of the potential), QvQr is the product of the vibrational and rotational 
partition functions and f{v,T) is the Maxwell-Boltzmann velocity distribution 
function 

/ \ 3/2 2 

f{v,T)=47ri--^j v'e-^. (1.33) 

Note that the electron velocity is usually much higher than the molecular trasla- 
tional velocity so that the relative collision velocity v is simply given by the elec- 
tron velocity; in line with that, it also follows that the reduced mass could be ap- 
proximated to the electronic mass, fie. However, to be precise, when the electrons 
and the molecules are described by two different mean kinetic temperatures, a 
convolution of two Maxwell-Boltzmann functions could be also use, as indicated 
in P3l . According to our theoretical approach to calculate the elastic electron- 
molecule ICS (see chapter 2), we will assume the target molecule with its nuclei 
fixed at the equilibrium geometry and, as our integral cross sections are rotation- 
ally summed, once the variable of integration is changed to the collision energy. 



the above equation (eq. |1.33[ ) becomes: 

This is the expression which will be used for the estimate of kf{T) in the present 
work. 

Moreover, since in a partially-ionized plasma like an interstellar cloud or a cir- 
cumstellar envelope the electron transport is primarily governed by momentum- 
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transfer cross sections (MTCSs) for elastic scattering of the present electrons, we 
find important to also provide realistic estimates of the elastic differential cross 
sections (DCSs) for electrons scattered off gas-phase molecules: despite the di- 
rect and useful physical information given by the DCSs, in fact, their knowledge 
enables us to estimate the MTCS which is defined as 

I'll 7 

amiE) = 2n (1 - cos 6) —{E) sin 9 d9 (1.35) 
Jo d\ I 



It therefore follows that eq. 1.33 expressing the metastable formation rate coef- 
ficient can be employed with our computed elastic (rotationally summed) ICS as 
well as with our computed MTCS, thus providing realistic upper and lower limits 
for the real kf value. 



1.4.3 The dynamical evolution of the metastable anion: pseudo- 
ID approach. Investigating the IVR occurrence 

Once the metastable negative species is formed, in order to follow the correspond- 
ing post-attachment evolution, one should need in principle a very detailed and 
computationally demanding treatment based on complex kinetic and/or statistical 
method [|44ll45l . We instead use a much more simplified treatment, but well con- 
solidated [|46l|471|48l, which is qualitatively indicative of the possible occurrence 
of the intramolecular vibrational redistribution (IVR) toward the stabilization of 
the full-size undissociated negative ion or to its fragmentation. 
The IVR, i.e. the process by which the vibrational energy initially localized in 
a particular mode is redistributed among all the vibrational modes of a given 
molecule, is an ubiquitous and foundamental process which has strong influence 
in the kinetics and in the final outcome of many chemical reactions, ranging from 
the unimolecular decomposition to protein folding. From the theoretical point 
of view, since the molecular vibrational modes cannot be assumed perfectly har- 
monic, they in principle can interact so that the unimolecular reaction under in- 
vestigation does occur when the superposition of these modes causes the reaction 
coordinate, on the corresponding potential energy surface, to reach the top of the 
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reaction barrier. In this connection, one should note that in the case of a metastable 
anionic state (with positive EA), the complex transient species essentially needs 
to reduce its internal energy content, since the required energy to reach the top of 
such a barrier is provided by the colliding electron. Following the Fermi's golden 
rule, the IVR rate, given the initially excited state, strictly depends on the root 
mean squared value of the coupling strength (y"^) as well as on the density of 
vibrational states, p 

T='^{V,'y{Eres))p{Eres), (1-36) 

Eres being the real resonance energy component. While predicting the density 
of bath states to a sufficiently good approximation can be currently accomplished 
[|49ll . even a rough estimate of the average coupling strength requires in principle a 
very accurate knowledge of the involved (complex, for metastable anions) poten- 
tial surface and thus extensive calculations that account for presence/absence of 
intermediates transient states [[50ll . large amplitude motion [[5T| as well as extreme 
motion states [l52ll . which can change rather unpredictably the expected size of 
couplings by several orders of magnitude. 

When considering the electron attachment process where the excess energy brought 
in the neutral target by the resonant electron gets transferred into the internal en- 
ergy of the molecular bonds, one is interested in gathering if the breakupof the 
target with the ensuing formation of stable anionic fragments takes place or, con- 
versely, if the creation of a long-lived full-size (meta)stable anion with a little 
residual energy occurs, the latter being in turn a good candidate for the subse- 
quent radiative stabilization by which the stable anion is eventually produced in 
the highly diluted, practically coUisionless, ISM. In this framework, we note that 
Herbst and coworkers [[54l 1551 assume an efficient attachment of s-wave elec- 
trons at zero energy, calculating the radiative electron attachment rates for lin- 
ear molecules like CN3 [[53 and HaC^/HC^H with n = 4, 6, 8 [|51. However, 
higher energy astrophysically relevant resonances still need to efficiently dissipate 
their excess energy before the autodetachment occurs, so that it seems reason- 
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able to consider as a very likely possibility the energy transfer within the molec- 
ular network as the molecular vibrational modes couple with the metastable elec- 
tron. It therefore follows that, according to a simplified pseudo monodimensional 
(pseudo-ID) picture for the electron attachment mechanism, one is searching for 
a possible dominant fragmentation path which finally occur via the beakup of 
that particular bond and the stabilization of a molecular anion. As will be dis- 
cussed with full particulars in connection with the ortho-benzyne and dicyanogen 
TNIs (see sections 3.1.6 and 3.3.3, respectively), we carried out several calcu- 
lations using different molecular structures for which the resonant electron den- 
sity maps showed, at the initial molecular equilibrium geometry, the presence of 
nodal planes and of partial localization along specific bonds of the resonant elec- 
tron density. The ensuing pseudo- ID potential energy curves, as cuts over the real 
potential energy surface associated with both the neutral N-electron target and 
the metastable (N-i-l)-electron TNI, follow a selected stretching deformation and 
provide the energy balance as 

= Eres{R) + En(R) - E^iRe,) (1-37) 

where E^ is the computed electronic energy (at the same level of accuracy of 
the scattering calculations at the equilibrium geometry) for the neutral molecule 
at a given deformed geometry identified by R, and Eres (R) is the real part of the 
computed resonant electron (complex) energy over the same range of molecular 
geometries. The corresponding widths, F^eslR). of the resonances associated with 
the (N-i-l)-electrons states are also given at each computed molecular geometry via 
the well known relation [|56l 

^cornp/ex(R) ^ ^^^^(r) ^ Zr,,,(R). (1.38) 
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Chapter 2 

The theoretical method 



2.1 Total and effective Hamiltonian 

In order to take up the theoretical treatment of any many-body system from the 
quantum mechanical point of view, as the corresponding Schroedinger equation is 
in principle responsible to describe all the system properties, the first step consists 
in writing down the associated interaction Hamiltonian in a suitable form so that 
the eigenfuctions for the states of interest can be either numerically or explicitly 
determined. In the electron-molecule scattering process the particles participating 
to the coUisional event are the incident electron, the electron of the molecule (ng) 
and the molecular nuclei (N^). Therefore the Hamiltonian of the system in prin- 
ciple shall depend on the 3+3ne+3N„ degrees of freedom. However, with some 
approximation it becomes possible to separate the motion of the molecular bound 
electrons from that of the other particles; in such an approach, thus, the wave- 
function determination is then reduced to the determination of a wavefunction 
with 3N„+3 degrees of freedom, thereby reducing the complexity of the scatter- 
ing problem. At the same time, in principle, such a procedure allows the nuclear 
dynamics during the collision event to be also expressly included. 
The total electron-molecule interaction Hamiltonian expressed in atomic units is 

Htot = + H^oi{r,^) + y,;f (re,r,R) (2.1) 
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with 



1 ^ 



y-(r„ r, R) = + - (2.2) 

1=1 ' ' j=l ' ■' ' 



and 



1=1 j=l J i=l j=l ' ' ^' 



i=l j=i+l ' ' j=l k=j+l 



Tn + 71^01. (2.3) 



In the above equations, r represents the positions of the molecular electrons, is 
the position of the colliding electronic particle and R refers to the positions of the 
nuclei. For the sake of simplicity, all the interaction terms are written in atomic 
unit. It is also possible to consider that the previous equations are written in the 
molecular center of mass which, on the basis of the huge difference between the 
free colliding electron weight and the molecular weight, can be qualitatively seen 
as coincident with the total system (molecule + incident electron) center of mass: 
in this connection we can reasonably neglect the spatial motion of the total system, 
focusing instead on its internal dynamical evolution. Moreover, in Hmoi the spin- 
orbit interaction has been neglected as well as other interactions involving the 
electronic and nuclear spin (like the spin-spin and the hyperfine structure term). 
If we now denote with 0(r, a, R) the wavef unction of the electronic and spin parts 
(the label a refers collectively to the molecular electron spins) of the molecule 
with eigenvalue £'n(R) within the Bom-Oppenheimer (BO) approximation, one 
can write the eigenfunctions of the total electron-molecule stationary Hamiltonian 
(eq. 1) as 
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*(re, r, a, R) = ^ |0„(r, a, R))V'„(re, R) 



(2.4) 



The above close-coupled expansion is exact, and also includes the continuum scat- 
tering states of the system (ipn)- For electronically elastic collisions, it is further 
possible to neglect in eq. 4 the contributions coming from the excited states (0„, 
n^^O) of the molecule, so that 



*(re,r,a,R) ^ |0o(r,(T,R))^/;o(re,R). 



(2.5) 



When writing the molecular Hamiltonian Hmoi as in eq.3, then the eigenvalues 
equation Htot^=E^ with the total Hamiltonian of eq.l applied to eq.5 will give 



I0o) 



A, 



(|(/.o(r, a, R))V'o(re, R)) = EroMi^o 



A, 



eMl 



where (R) refers to the ground molecular electronic state energy within the BO 
approximation. Multiplying now from the left for the ground electronic molecular 
function (0o I and integrating with respect to the position variables of the molecular 
electrons (as indicated by the label r) one obtains 



+ (^olTivl^o). + Tat - 5](0o| Vr|0o). • Vr - - + {MVin"^ \<l>^) 



1=1 



(2.6) 



having esplicitated the free electron kinetic energy as = (Etot - ^o'(^))- 
This is the complete integro-differential scattering equation, where also the nu- 
clear dynamics is susceptible to 'feel' the incoming electron: since the molecular 
ground electronic state BO wavefunction is given by 0o(r, cr, R) = cpf (r, a)xj,i,(R), 
Xj,u(R) being the nuclear wavefunction, the second and the fourth term on the left 
hand side implicitly contains the adiabatic corrections to the BO potential surface. 
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It obviously follows that using the complete close-coupled expansion as in eq.5, 
the same terms will contains also the non-adiabatic corrections. One should note 
at this point that the BO approximation, previously introduced to treat indipen- 
dently the electronic and the nuclear molecular degrees of freedom, is the most 
common and natural one which is encountered in any coUisional problem. It fol- 
lows that, when treating the scattering problem including the nuclear dynamics 
inside the scattering equation, as in eq.6, the BO electronic energy ^o'(^) ^^^^ 
indeed provide the effective potential energy for the nuclear motion. 
Furthermore, when making the approximation of eq.5, for the scattering wave- 
function ^ the most important terms which have been neglected are 



(0o|l^,r(re,r,R)|0n). 

They represent what is classically described as the polarisation interaction. Hence 
the eq.6 takes in consideration the projectile interaction with the molecular elec- 
tronic cloud in its unperturbed ground electronic state, and does not consider 
the deformation induced by the approaching free electron. So, identifying the 
(0o|Knt^|0o) term with the static interaction, the neglected terms consequently 
represent just the polarisation interaction 



Kir = J2(^o\Vir{re, r, R) = J2 ^^nf'"" (2.7) 

As shall be discussed in detail when treating the single center expansion method, 
this interaction can be brought back into the eq.6 using a modellistic form; for 
the same reason, it is foundamental to point out now that although we are deal- 
ing with fermions (the molecular electrons and the colliding particle), we have 
not suitably antisymmetrized the expansions in eq.4 and/or eq.5. Such a physical 
constraint cause the complete scattering equation (eq.6) to contain also non-local 
potential terms whose meaning is thus justified on the basis of the Pauh princi- 
ple; on the other side, such constraint would cause the same scattering equation 
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to be rather impossible to be computationally resolved: for the same reason as 
above, also the exchange interaction contribution will be introduced using a suit- 
able optical potential (which will be described later). In other words, if in the 
Schoredinger equation with the total Hamiltonian given by eq.l we would expand 
the total wave function in terms of all possible target states (rotational, vibra- 
tional and electronic), and then the antisymmetrizing operator would be applied 
to ensure the validity of the Pauli principle, the ensuing coupled integro differ- 
ential equations will automatically include all types of acting forces without any 
approximation. However, on the other hand, such an expansion would be not 
entirely feasible since in this way the equations will contain non-local and often 
complex terms that would be either very complicated to handle numerically or 
too much demanding computationally. Therefore, the way we follow to bypass 
this obstacle consists in resorting to semiempirical method thus using a DFT ap- 
proach for the short range forces and a local (energy-dependent) model optical 
potential for the exchange interactions. For the purpose of the present section, in 
conclusion, if one writes 

V^fir^, R) = V^Zire, R) + V;^f{r,, R) + (r^, r, R) (2.8) 

then the complete 'exact' scattering equation including the adiabatic corrections 
takes the following form: 



Ar 



N„ 



+ (0o|^iv|0o)r + Tiv - ^(0o|Vr|0o)r " Vr 



i=l 



2 



+ v;;f(re,R) 



^0 = 
(2.9) 



2.2 Fixed nuclei approximation. 

The most simple approximation that is possible to use for the full Hamiltonian of 
eq.9 is the one where all the nuclear motions of the target molecule are neglected, 
so that the total FN Hamiltonian appears in the following simplified form: 
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ttFN 



2 

Thus the scattering equation is given by 



(2.10) 



2 



-^ + ^.nf(re,R) 



V'o'^^^O (2.11) 



2 2 

However, this approximation is good only when the energy of the projectile is 
high enough that in the spatial region in which there is interaction, the electronic 
projectile spends effectively little time so that the molecule does not manage to 
move during that time. In other words, such an approximation should be consid- 
ered strictly valid for direct and non-resonant scattering. In this sense, then, the 
FN approach is expected to be reliable, provided that the kinetic energy of the 
incident electron is appreciably larger than the molecular rotational energy dif- 
ferences of interest. On the other side, when not excited, complex and 'heavy' 
polyatomic molecules are considered as target, the latter thus being characterized 
by relatively slow nuclear (vibrorotational) motions, due to the problem complex- 
ity, it still makes sense to describe the collision event within the framework of 
a FN description. Moreover, for very low-energy incident electrons (E < 1~1.5 
eV), the colliding particle usually has not the kinetic energy needed to deeply pen- 
etrate the molecular electronic cloud, so that its coupling with the molecular nu- 
clear dynamics can be viewed as partially negligible, as we will in fact argue for 
the electron-coronene case (see chapter 3, section 3.2, 'The C24H12 molecule'). 
However, to be precise, for very low incident energy, the collision time Tc may 
exceed the vibrational period Ti, but it might still be such that Tc ~ r^, being 
the molecular rotational characteristic period. Since the above condition might be 
considered as satisfied also at thermal energy (for low-lying rotational states in 
fact the energy spacing can be of the order of 10^^ ~ 10~^ eV), one should thus 
expect that rotationally inelastic collisions play a role in the scattering process; 
more in general, therefore, one should expect that vibrationally (or vibrorotation- 
ally) inelastic collisions play a role in scattering and also in resonant processes 
involving large molecules, where this actually takes place but expecially in the 
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differential cross sections (DCS) behaviour: according to the previous qualita- 
tive argumentations, we thus simply surmise that for the purposes of the present 
investigation such inelastic contributions, even though (very) likely, are not so im- 
portant, as then confirmed by our findings for the coronene molecule (see chapter 
3, sec. 3.2.4 and 3.2.5). 

When considering polar molecular target, the FN approximation also fails. The 
detailed description for the above breakdown would deserve more complex and 
detailed analysis being thus beyond the purpose of the present section as well as 
of the whole PhD dissertation. It will however be considered again, without enter- 
ing the complex theoretical details, when the conclusions for the ortho-benzyne 
molecule shall be summarized (see chapter 3, sec. 3.1.9). For the moment it is 
sufficient to emphasize that, once the T matrix is introduced to calculate the rel- 
evant scattering cross sections, it shows a very slow convergence in the partial 
waves number. It is possible to glean a reliable justification for this, using a sim- 
ple and qualitative physical picture. High / partial waves are mainly 'restricted' 
to the long-range spatial region due to their intense centrifugal barriers so that 
their phase accumulation (Q) is essentially due to the contribution associated to the 
electron-permanent dipole interaction. The latter statement, in fact, is easily con- 
firmed when keeping in mind that polar molecules are characterized by the pres- 
ence of a permanent dipole that has with the incoming electron a long-range inter- 
action (V^[°'''^ ~ 1/r^) stronger than those given either by the static quadrupole 
(^nt"*"^"'"^ ~ 1/r^) or associated with both the isotropic and anisotropic polar- 
izabilities (i.e. the induced dipole: V^'^^ ~ 1/r^). In conclusion, when the dipole 
is not allowed to rotate, it is possible to show that the associated DCS diverges 
in the forward direction, causing then the ICS to be also divergent, even though 
logaritmically. This means that, as discussed in section 3.1.3 in connection with 
our findings for the ortho-benzyne molecule, the larger partial waves number the 
more rehable the ICS values are. 



the phase shifts Si, see chapter 1 
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2.3 Single determinant expansion. The SE, SECP 
and SMECP approximations. 

The purpose of the previous sections was to derive the total hamiltonian for the 
electron-molecule scattering problem and consequently to introduce the FN ap- 
proximation thus separating the nuclear and the electronic dynamics during the 
collision event. Now we move to separate the scattered electron coordinates from 
the molecular bound electrons: this passage is complicated by the fact that the 
1 electrons involved in the present case are undistinguishable fermions so that they 
are affected by specific restrictions in their wavefunction which, asaccording to the 
Pauli principle, have to be suitably antisymmetrized with respect to the exchange 
of each couple of electrons. The final aim is to write down the scattering equa- 
tions when the static-exchange (SE), the static-exchange-correlation-polarization 
(SECP) and the static-model-exchange-correlation-polarization (SME CP) approx- 
imations are made, respectively. This will enables us to illustrate in section 2.4.6 
the meaning of the adiabatic-static-model-exchange-correlation-plarization (AS- 
MECP) approaches which have been also used to obtain some of the findings 
about the molecules involved in the present work. 

For simplicity sake, we now omit the esplicit dependance of the scattering wave- 
function from the nuclei position ^ so that, introducing the antisymmetrization 
operator A, eq.4 can be rewritten as follows 

^(re,r) =i^|0„(r,a))7/;„(re) (2.12) 

n 

It is now possible to make the following distinction, depending on the number 
of terms which in the previous summation are retained: when only the first one 
is considered (so that only the molecular electronic ground state is taken in con- 
sideration for the collision process) we obtain the SE approximation, otherwise 
the SEP. In what follows the associated scattering equation will be derived, at 

^having assumed the FN approximation to be valid, it is implicit for what follows that the 
scattering wavefunction depends parametrically on the nuclear geometry R 
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the same time introducing the second important constraint that is crucial for our 

theoretical model: the single determinant expansion. 

For electronically elastic collisions, only the first term in the expansion provided 
by eq. 2.12 can be considered. When the molecular electronic ground state wave- 
function 0^^'^' is expressed at the Hartree-Fock level as a single Slater determi- 
nant, one can write 

*(re,r) =i||(^i(l),(^2(2),...,(^„,(ne)||^oK + l) (2.13) 

where the symbol 1 1 . . . 1 1 refers to the Slater determinant, having thus indicated the 
i-th one-electron spin-orbital as (pi{i). It then follows that the esplicit form for the 
A operator is given by 

A = , y {-lyv, (2.14) 

V being the permutation operator acting on each possible couple of electrons. The 
FN-SE equation for the scattered electron wavefunction ?/^o(^e) is then obtained 
by projecting the stationary FN scattering equation H^qj,^!^^ = E^qj."^^^, con- 
taining the FN Hamiltonian (eq.ll) and the SE expansion (eq.l3) for onto 
the bra { ||<^i(l)</?2(2)...</7ne(^e)|| | associated to the closed shell electronic 
molecular ground state, so that the successive integration on the bound molec- 
ular electronic coordinates provides (in a.u.) 



■lvl-lkl + Vstat{re,R) 



- E / ^rVKr')rr^^o(r') ■ v,{r,) (2.15) 
j=i J l^e r I 



with 



K,„,(r„ R) = E / -(rOr— ^^'.(r') - ^ — ^ (2.16) 

1=1 I e I I e k\ 



and 
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k 



e 



2 [e^St'"' - eo{R)) 



i/l^||^i(l)^2(2)...9:'neK)|| = 6o(R)||^i(l)^2(2)...9:'„.K)|| (2.17) 

On the basis of the adopted approximations, it is easy to recognize that the SE scat- 
tering equation does contain only the static (eq. 2.16) and the exchange interaction 
(given by the sum on the r.h.s of eq. 2.15). The SECP approximation level is then 
obtained from he previous one, but now, besides the static and the exchange inter- 
actions, one takes in consideration also the dynamical effects given by the molecu- 
lar response to the colUding electron by introducing an optical phenomenological 
local potential Vcp(re,r) to include the molecular polarization; consequently, the 
scattering equations take the following form 



-^V^^ - ^kl + Vstatire, R) + Vcp{re, r) 



= E / drV,{r')-^Mr') ■ ^A^e) (2.18) 

3=1 l^e r I 

To conclude the present section, we note that the previous are clearly non-local 
equations; it is then possible to further simplify them by replacing the sum of non- 
local terms on the r.h.s. with another local energy-dependent potential Wexi^] £coii, 
Ip). The latter will be described Uttle more in detail within the framework of the 
single center expansion (see sec. 2.4.4). The resulting SMECP equations finally 
represent our local differential scattering equations: 



Vr, - ^^e + Vstatire, R) + Kp(re, r) + Wex(r; ScolU Ip) 



(2.19) 



2.4 The Single Center Expansion 

In the present section, we shall illustrate in some detail the guidelines of our 
method by means of which we computationally solve the scattering equations 
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describing the elastic resonant collision process; the method is based on the gen- 
eral single center expansion (SCE) approach adapted for electron- polyatomic 
molecules elastic scattering. Generally speaking, however, the scattering pro- 
cesses which occur in the collision of electrons with polyatomic targets have been 
object of intensive theoretical studies in the past years, so that in order to solve 
the final scattering equation several methods do indeed exist. As previously men- 
tioned, our method is based on the coordinate representation of the interaction 
forces, so that a little detailed description for the local model exchange as well as 
for the static and both the short-range (correlation) and long-range (polarization) 
interaction terms will be illustrated. 

2.4.1 The SCE bound state wavefunction 

The first step in the SCE method consists in obtaining the bound state wavefunc- 
tion that will be later included into the scattering equations. One therefore needs 
to interface with an initial quantum chemistry code (like Gaussian, [HI) that is 
thus employed to generate the Single Determinant description (at the Hartree- 
Fock level) of the target electronic wavefunction. As mentioned above, in this 
way we obtain the closed shell electronic wavefunction for the neutral target as a 
single Slater determinant of one-electron bound spinorbitals, indicated as (fi{r; R) 
with i = 1 : Ue, or equivalently of ne/2 two-electrons multicenter MOs labelled 
as <^i(r^;R) with i = 1 : ne/2 and ^ = 1 : 2. Each of the above functions is 
in turn susceptible to be expressed as a truncated linear sovrapposition of opti- 
mized known functions gk{r;R) (Gaussian Type Orbitals, GTOs), depending on 
the choosen basis set, according to 



A=l k=l v=l 

where i refers to the i-\h multicenter molecular wavefunction, and the /c-th primi- 
tive gaussian function for the A-th atomic center is 




(2.20) 
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(2.21) 



N{a, h, c; 7) being the normalization coefficient 




3/2 




-,1/2 



J\f{a,b,c;-f) = 



(2a-l)!!(2&-l)!!(2c-l)!! 



(2.22) 



The C^^ is the GTO coefficient of the u-th GTO at a given molecular geometry 
as concisely expressed by the argument (Ra). The index u labels each primitive 
gaussian function within the subgroup that belongs to a given contraction coef- 
ficint d^''^, so that the contracted gaussian function G^''^ is 



The single center expansion method can be now introduced. Since we want to 
generate the full electron-molecule interaction potential as a function of the elec- 
tronic density p(r) of the target molecule in a suitable way, the crucial step con- 
sists therefore in expanding both the bound electron functions and the continuum 
scattered function around the center of mass 



In the above formulas, i labels a specific mutlicenter MO, which contributes to 
the bound electrons density p(r) in the polyatomic target, while the p/i indices 
refer to the ;U-th component of the p-th relevant irreducible representation (IR), 
respectively. The h index labels a specific basis, for a given angular momentum £ 
value, for the p-th IR that is under consideration. 

At this point, a crucial passage which fully deserves to be pointed out consists 
in the construction of the symmetry-adapted, generalized harmonics indicated as 




(2.23) 




(2.24) 



(2.25) 
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Xfl^{9,(l)). Ror a given £, the associated 2i + I spherical harmonics Y[''{6,(j)) 
form a basis for a (2£ + l)-dimensional IR of the full rotational group, and at the 
same time they constitute the basis of a reducible representation of the molecular 
point group to which the considered molecule in the choosen frozen geometry 
belongs: the molecular point group is in fact a subgroup of the full rotational 
group. Therefore, the generalized harmonics can be defined as basis 

for IRs of the molecular point group so that they can be reasonably expanded as 
linear combinations of spherical harmonics 

XZ{0,ct>) = Y,lfl^Yr{e,ct>) (2.26) 

m 

where the 6^^^ coefficients, choosen so that the X functions transform as the /i-th 
component of the p-th IR of the molecular point group considered, are described 
and tabulated in [|2l. Moreover, since the X'^l^{9, 0) functions satisfy the following 
orthonormality relations 

J xf/*{e, (j))x^i:{e, <f)) sine de d<p = v^/^m'^w (2.27) 

then the 6^^^ coefficients constitute a unitary transformation between the Xs and 
the Fs, in the sense that V£ it can be shown 



X] ^ih/r^^ihm " ^pp'^l^l^'^hh' (2.28) 
m 

KL'Kim" = ^ra'm" (2.29) 

p,lJ,,h 

Thus, writing the multicenter electronic wavefunction which describes the target 
neutral molecule as an antisymmetrized product of multicenter MOs ipi as in eq. 
2.20, it is then possible to generate each of the radial coefficients u^j^'\r) for the 
bound molecular electrons by the following numerical quadrature 

= E E E E / / ^i^(^) 

A=l k=l u=i m=-e-^^ 
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2.4.2 The SCE one-electron density 

Having obtained the radial coefficients u'^^ir) from eq. 2.30, each bound one- 
electron wavefunction 



r 
eh 



is then expanded around the center of mass in terms of Xf}^ functions so that the 
one-electron density function p(r) can be written as usual in the following form 



p(r) = j\ (iet||(^i(r)(^2(X2)<^3(X3)---¥'ne(Xne)ll f dX2 dxs ...dXn, 

Tie 

- 2-J2\^iir)f (2.31) 

i=l 

where the numerical factor 2 is due to the sum over the spin, and the sum is over 
each i-th doubly occupied orbital. At this point it should be noticed that the p(r) 
function belongs to the total-symmetric (Ai) IR. Therefore, once the quantity p(r) 
is obtained as above explained, it can be expanded in terms of symmetry-adapted 
total-symmetric general harmonics: 

Pi^) = I E PMm)ir)Xi\^^{9, cf>) (2.32) 

eh{m) 

where 

p,^ = 2 • E / sin(^) dd / d(t> ^i{T) ■ (^,(r) (2.33) 

The knowledge of these coefficients, in principle, enable us to numerically es- 
timate all the interaction terms between the incoming electron and the target 
molecule. Let us see how. 



2.4.3 The SCE static potential 

The static potential exerted by the molecular electrons and nuclei on the surround- 
ing molecular volume is given by the sum of two contributions 
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Ktat(r, {%}) = / (2.34) 
./ |r — s| ^|r — Rjl 



K/(r) + K«c(r,{R,}) 



the first of which is the repulsive electron-electron contribution, and the second, 
with the sum over the number of atomic centers N^, refers to the nuclear attractive 

effect on the colliding electron. The substitution of eq. 2.32 and 2.33 in the above 
formula (eq. 2.34), after using the following generalized form of the well known 
spherical harmonics sum theorem. 



3^ = E E ')^T{Or, ^r)XZ{es, ^s) (2.35) 

' p/i (.h 



with 



enables us to obtain the electronic contribution to the static potential expressed as 

ye;(r,{R.}) ^Y.^-''^^^'')^thi&rAr). (2.37) 

where the potential radial expansion coefficients Vei,£h{f) are obtained via radial 
integration 



47r 

yei/h{r) = —— I ds pih'ye{r,s) (2.38) 



2i + l 

It therefore follows that the single center expansion of the nuclear static contri- 
bution Vnuc{^), using once again the above expressions 2.35 and 2.36, is given 
by 



Vnucir, {Rj}) = VnucMr, {Or, 4>r) (2.39) 

Eh 

with the radial coefficients assume the expression 

VnucMr,{^j}) ^ ^^J2^Mr,Rj)Xi^{0n^,<pR,)- (2.40) 
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2.4.4 The exchange interaction: the Free Electron Gas Exchange 
(FEGE) Model Potential 

To circumvent the difficulty of having to handle and solve integro-differential ra- 
dial scattering equations, the necessary exchange effects can be suitably modelled 
through the use of simpler, energy dependent potential forms which can be added 
to the exact static potential discussed in the previous section within the framework 
of the SCE approach. One of the most widely used for polyatomic molecular tar- 
gets, and that is also used for all the scattering calculations of the present thesis, 
is the Free-Electron-Gas-Exchange (FEGE) introduced many years ago by Hara 
O. The actual form of this exchange potential derives from two approximations 
in the integral exchange terms of the exact-exchange scattering equations (see eq. 
2.15, 2.18 and 2.19 in sec. 2.3). First, the molecular electrons are treated as a 
free-electron gas, with a given charge density p(r, R), determined for particular 
electronic state of the target under investigation (initially the ground electronic 
state). Second, the distorsion of the continuum function is neglected so that the 
colliding electron is treated as a plane undeformed wave. The resulting FEGE 
potential is then given by the following form: 

VpEGEir, {R,}) = --Kpir, {R,}) 

TT 

where the wavevector up to the top of the Fermi surface is given by the usual 
free-electron relation 

Kp{r, {R,}) = [37rV(r, {R,})] (2.42) 

and T] holds the ratio between the actual wavevector k for the scattered elec- 
tron and the one at the top of the Fermi surface for the electron gas, rfir) = 
k{r)IKp(r, {Rj}), where the local momentum is given by 

k{r)= [2{e + Q + Kl{r,{^,})\''\ (2.43) 



1 1 

2^ ~ Ar] 



7] 



In 



(2.41) 
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The collision energy e is in turn given by the initial-channel energy and Ip refers 
to the first ionization potential for the neutral target molecule. The final FEGE 
exchange is then obtained by means of the integration scheme already outlined 
and used in the above subsections for the bound state wavefunction and for the 
electron density: 

VpEGEir, {R,}) = VFEGEAr)Xth- (2-44) 
2.4.5 The correlation-polarization (Sep) potential 

Another important type of interaction between the impinging electron and the tar- 
get molecule comes from the description of the response function of the target to 
the incoming charger projectile. Due to the fact that this interaction depends on the 
distance of the incoming electron in the sense that it varys from a purely polar- 
isation effect at large distances up to including exchange-correlation interaction 
with the bound electrons whithin the molecular volume, the model correlation- 
polarization potential Vcp that is implemented in our SCE treatment of polyatomic 
molecules makes a distinction between a long-range region and a short-range of 
interaction. Within the former, the perturbative polarisation effects are dominant, 
the potential is correctly susceptible of a local analitic form and is adiabatic on the 
incident electron's velocity; conversely, within the latter, both non-adiabatic and 
non-local effects play an important role. It therefore follows that one can write 
the overall correlation-polarisation interaction as 

Kp(r,{R,}) = Kor.(r,{R,}), r<r, (2.45) 
= V^oKr, {Rj}), r>rc 

where the connecting spatial factors Tc are usually obtained from the crossing radii 
of the lower coefficients in the SCE expansion 



v;,(r,{R,}) = ^v;,,,,(r)x/,^ 

ih 



(2.46) 
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The short-range correlation term 

The short-range contribution, Korr^ rest on a DFT treatment a full description 
of which has been given elsewhere [4], so that here we only briefly outline its 
leading features. As the relevant quantities reported so far, the Vcorr potential 
must be expanded over symmetry adapted spherical harmonics 



V,orr{r, {Rj}) = J2 ycorrMr)X^i:{d, 0) = ^ [^c(p(r, {R,}))] (2.47) 



where Ec is the correlation energy which depends on the one-electron density 
p(r). Apart from the one-electron density, it can be shown tha the Vcorr i^) poten- 
tial also depends on the first and second functional derivatives of p(r): accord- 
ingly, unlike the other potentials so far described, the K;orr(r) requires a further 
single center expansion of the bound state wavefunction. 

The long-range polarization term 

Within the outermost asymptotic region we use a long-range potential which, de- 
pending mainly on the static electrical properties of the neutral molecular target 
(i.e., dipole and higher moments and associated polarizabilities) is susceptible of 
the following general form for a given molecular geometry {Rj}: 



Let us first focus on the case of a permanent dipole moment and of its associated 
polarizability a. When only the symmetrical spherical polarizability is known 
(often referred to as the average dipole polarizability aav = 1/3 (aa^. + ayy + azz)) 
it is customary to assume that aav = dxx = c^yy = c^zz while a^y = Oi^z = oiyz = 0. 
In this case the polarizability tensor is written as 




(2.48) 



a, 



av 











(2.49) 
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On the other hand, if all of the six terms of a are known, then the polarizabiUty 
tensor becomes 





2^xy 


2 ^xz 


2^xy 


(Xyy 


2^yz 


2 '-'■xz 


2^yz 


C^zz 



A common case is when only the ckq and a2 polarisabilities are known: in this 
case the polarisation potential, for a linear molecule, has the simple form 



V,oi{r,{^,}) = - |;|[3co«'(^) - 1] (2-51) 

where the angular dependence is provided by the Legendre polynomial P2[cos{9)], 
and a XX = oiyy = (ckq — 1/2q;2) while = (ckq + q;2), the other off-diagonal terms 
of the polarizability tensor being zero. Besides this special case, one has to take 
into account all the components of the dipole static polarizability so that the long 

range polarisation potential expressed in atomic units assumes the form 

where = x, y, z for i, j = 1, 2, 3, respectively. The last equation is that one 
implemented with the options to select all the different polarizabilities reported 
above, as a function of the available data. Once the Vpoi{x, y, z) is calculated and 
expanded over the symmetry adapted spherical harmonics belonging to the Ai 
irreducible representation 



V(r, {R,}) = E ^voiAr)Xt^ <P) (2.53) 

eh 

the matching point Tc is found selecting the matching direction along the line 
connecting the molecular center of mass and a given atom or, equivalently, along 
the line connecting the center of mass to a generic point with cartesian coordinates 
X, y, z provided by the user. 
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2.4.6 The ASMECP level of approximation 

In order to examine the mechanism and the qualitative characteristics of a possible 
low-energy, one-electron resonance as well as its evolution along the lower part 
of the molecular potential energy curve according to our pseudo-ID model for 
the resonance energy redistribution during the temporary anionlifetime, we need 
a model which is simple enough to be computationally attractive and which at the 
same time includes sufficient details of the full scattering problem to reproduce the 
essential features of the physics involved. Thus, we look at the low-energy res- 
onances also by using a simple, purely local model potential thatwe have called 
the adiabatic static model-exchange correlation-polarization (ASMECP) potential 
[|5l. We start by noting that the standard, symmetry adapted angular momentum 
eigenstates, X^/^*, do not form the most compact angular set for the e-molecule 
scattering problem: an alternative basis expansion is provided, in fact, by the an- 
gular eigenfunctions obtained by diagonalizing the angular Hamiltonian at each 
radius r. These distance-dependent angular eigenstates are referred to as the adi- 
abatic angular functions Zl^{9, 0, r): at each radial value, they are linear combi- 
nations of the symmetry-adapted 'asymptotic' harmonics discussed before 

Zne, 0, r) = 5^ X^i:{e, <P)CeHAr) (2.54) 
eh 

where the expansion coefficients are solutions of the matrix eigenvalue equation 



E [vrh^ehir) + See'S,J-^^j C^.^r) = Vr^{r)Ceh' ,k{r) (2.55) 

The eigenstates V^^''^(r) now form an adiabatic radial potential for each index 
k over the selected range of the electron-molecule distances (for the meaning of 
^e'h' ehi^) ^1- 2.58 in the next section). The spatial extent of the resonant- 
wave function can be determined from the well and angular momentum barrier 
of such adiabatic potential terms and the physical mechanism for the resonance 
is that of a trapped electron tunneling through the potential barrier. In order to 
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avoid the nonadiabatic coupling terms between adiabatic curves, we actually em- 
ploy a piecewise diabatic (PD) representation of the potential whereby the radial 
coordinate is divided into a number of regions so that the i-th sector is defined 
as rj_i < r < rj+i, with tq = 0. In each radial region we average the coupling 
potential Viih',ehir) over r and the resulting averaged potential is diagonalized as 
to yield a set of angular functions Zl'^-{9, (p). Then, in the i-th region 



m eq. 



2.55 



the scattering potential is transformed into the new representation in which it is 
nearly diagonal. The resulting equations are solved using the full scattering po- 
tential in each region with the further approximation of ignoring the off-diagonal 
couplings in that region: to solve the radial equations using the PD approach 
requires matching of the radial functions and their derivatives at the boundary be- 
tween radial regions. The transformation of the radial functions from one region 
to the next is given by the transformation matrix ^, defined by 



^tP^'^ = J:C^'cI,^.'- (2.56) 

ih 



When the size of the angular momentum eigenfunction basis used is larger than 
the size of the diabatic angular basis set, the transformation matrix Vk,k' is not 
in general unitary: we however accomplish the unitarization of U^*^/^*^ using 
simple Graham-Schmidt orthonormalization on the columns of that matrix. 

It is now possible to conclude the present short but important section anticipat- 
ing that the above (ASMECP) treatment, being easier in terms of computational 
cost, usually finds the resonances (located as maxima in the integral cross section, 
each associated to vr-jump in the corresponding eigenphase) slightly shifted from 
those given by the full coupled-channel solutions [|5l which will be derived and 
discussed with full particulars in the next section (see eq. 2.57). However, their 
nature and general features remain unchanged and of the same physical signifi- 
cance, while allowing for a reduction of the computational effort. 
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2.4.7 The scattering equations; S and K matrices 

To derive the equations describing the radial motion of the scattered electron, as 
previously introduced, we make two approximations. Firstly, we assume that the 
fixed nuclei approximation, where the main symmetry molecular axis is taken 
to be stationary during the collision and pointed toward the incoming electron, 
can be made; secondly, we adopt a single centre expansion of both the bound 
molecular orbitals and the scattered electron orbitals about the centre of mass. 
Having in the previous sections already discussed both of them, now we move 
to write down and to solve the relevant radial scattering equations the solutions 
of which enable us to numerically calculate the K matrix (and consequently both 
the T and § matrices) that in turn provides the scattering amplitudes and thus the 
scattering cross sections. 

Once all the interaction potential terms (Vgtau Vfege, Vcp) as well as each of 
the N+1 one-electron wave functions (both the bound and the free particles) are 
recasted according ot the SCE method, their substitution into the eq. 2. 19 followed 
by multiplying the whole equation by X^,^, * and finally integrating all over the 
angular coordinates provides the following ensamble of coupled radial differential 
equations for the radial scattering function 



where i and j (the same for k) label a couple of angular channels iihi, ijhj among 
the ric coupled angular channels, so that 



^ffir) = {Xr \Vstat{r, {R,}) + V^EGEir, {%}) + V;,(r, {R,})| ) 



sin Or dOr dct>r *(^., 0.) [Ktat + VpEGE + Kp] ^fl^r, ' 



E / dWrXr*{er.ct>r)Xl^{erAr)X^/{erAr) (2.58) 
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For Tie coupled differential equations (for each state p/i) there will be a set of 
independent solutions which satisfy the 'internal' boundary condition /(O) = 0: 
this is the reason according which we introduce the additional label j by which 
we therefore write the solutions as ffj^. The § matrix is defined by looking for 
solutions with the asymptotic form ^ |71 

fPf{r) ^ e-'^''-^'^^6i, - §f;e+'('=^-t^^). (2.59) 

We can form a general solution which is a linear combination of these inde- 
pendent solutions 

rtc 

^r = E<^'w (2.60) 

from which we can then construct the total function for the scattered electron 

-^w = E E -^r(o^r.,(^r, 0.) (2.6I) 

i=l p/i 

The af^ numerical coefficients in eq. 2.60 are choosen such that the outgoing 
(r — 7- +00) boundary condition is satisfied, namely 

T{r) ^ e^*^"-^ + fM^o ■ r) (2.62) 

r 

where the fpNiM-o ■ r) is the FN scattering amplitude, and ko, r are the initial and 
final directions of the electronic projectile. Taking advantage from the following 
expansion for the plane wave 

e^k-r = AnJ2i: .^'^Xr(k)X-(r), (2.63) 

£ h,pfi 

obtained starting from the expression for a plane wave as a function of spherical 
Bessel functions and Legendre polynomials in which one introduce the spherical 
harmonics sum theorem coupled with the asymptotic form for the spherical Bessel 
functions j/(fcr)(|^, one can easily deduce the explicit form for the Oj^ coefficients. 
In fact, by writing the asymptotic total scattering function J^(r) as 
^je{kr) ^ sin(fcr — for r +00 
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J^(r) ^ J^.ncir) + J-.eatt(r) (2.64) 

one should first substitute into this expression the equations 2.59-2.61 as well as 
the asymptotic expansion 2.63 for J^j„c(r). Now, even though Fscatt{v) is not 
still known (as the scattering amplitude /FAf(ko ■ r), too), we know that e'""^^ 
and e^'^^' are two linearly independent functions as well as J^scatt{^) should not 
have any incoming wave component behaving like e"*'^'', so that by equating the 
coefficients of e^*^** on both sides of 2.64 it is possible to show that the expansion 
coefficients ^ given by 

«r = (-l)^'^|^^r^,(ko). (2.65) 

Now the scattering amplitude fpN can be obtained [8l by first substituting in eq. 
2.62 all the equations 2.59-2.61 (0) and then by taking the difference between the 
asymptotic total function J-'(r) and the expansion 2.63 for the plane wave, so that, 
finally 



/^^(ko ■T) = '£Y. ^^''''^^r(ko)Xf (r) - 5,,) . (2.66) 

As computations are usually done with real functions, we therefore introduce the 
IK matrix which is defined from the following real boundary conditions 

fffir) = sin (fcr - S,, + Kff cos (fcr - , (2.67) 
where the § and K matrices are related by the following transformation 



(2.68) 



§ = — 

I-iK 

2.4.8 The cross sections in the molecular frame 

The integral elastic cross section, often improperly called the total elastic cross 
section, a, describes the scattering from randomly oriented molecules. It is ob- 



'^remember the factor ^ as reported in eq. 2.25 
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tained by the integration of 



over all the orientations r of the scat- 



tered electron and then averaging over all the orientations ko of the impinging 
beam 



a = 



1 



dkn dr 



fpNi^o ■ r) 



(2.69) 



which is easily shovm to be equivalent to 



71 



EE IS 



(2.70) 



It is also instructive to express a in terms of eigenphases r^'^ which are in turn 
defined via the further diagonalization of each matrix corresponding to each 
contributing molecular symmetry. We could obtain a unitary transformation U 
such that AP'^ is consequently susceptible of a diagonal representation: 



/ p2<^(fco) 








(2.71) 



g2i)7^^(fco) J 



y • 

Another equivalent way to define the eigenphases r^f^ is via the eigenvalues of the 
¥J''^ matrix, which is hermitian and therefore has real eigenvalues. Likewise there 
is an unitary transformation by which the K^'* is diagonalized as 



( tanr^f (A;o) 

tan?72'^(A;o) 








(2.72) 



V ••• tan77^^(A;o) ) 

Having defined the diagonal elements of Ap^ to be given by exp(2i77f '*), where the 
Tyf are the eigenphases for each contributing symmetry \pn), it therefore follows 
that it is easy to show that 



47r 



(2.73) 



Ptj, 
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To obtain an expression for the differential cross section one needs to introduce 
then a coordinate system S' which is space-fixed (SF), the so-called LAB system, 
and another coordinate system S which is fixed to the molecular highest symmetry 
axis and that is identified as body-fixed (BE). The latter is obtained from S' by a 
rotation through the appropriate Euler angles (a, (5, 7). The previously discussed 
expression for the EN scattering amplitude (see eq. 2.66), through the symmetry 
adapted functions X^^(ko)) and Xj^(r)), contains angular factors like Y£m(ko) 
and Y^„i{^) that are in turn defined with respect the BE system, S (kg and r are 
defined with respect the S reference). The corresponding expression in the SE 
frame is obviously given by 

Y,^{v') = DL^'ic^, /3, l)Y,^ir) {2.1 A) 

m' 

where r is defined with respect to S' and r with respect to S and the coefficient 
/3, 7) are the Wigner functions S. The differential cross section for 
scattering by randomly oriented molecules is then obtained by averaging the SE 



quantity 



/5F(k'o-r; a,/3,7) 



over all the molecular orientations (a, (3, 7): 



7 \ f 

^{k,-r') = — J dadpd^ sin P\fsF{k,-r'; a,P,^)\\ (2.75) 



This angular integration can be done analytically and lead to an expansion of 
differential cross section in terms of Legendre polynomials Pf. 

^(k'o -r') = J2M^inPc{cose) (2.76) 
c 

where the description with full particulars inherent the calculations of the Ai co- 
efficients can be found in the first appendix and particularly in [fTOl[TT| . 
To conclude the present section, it is useful to emphasize the fact that a body-frame 
treatment like that one briefly illustrated in this chapter, is always appropriate as 
long as the incident electron wavefunction overlaps the molecule, that is to say in 
the so called 'inner region': in this spatial region, in fact, the colliding electron is 
moving faster than the nuclei. A qualitative explanation is found when one notes 
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that in general, due to the influence of molecular attraction, the projectile velocity 
varies during the collision in the sense that it is reasonable to say that an incident 
electron which is slow at large distances becomes faster at short distances. It 
therefore follows that the whole three dimensional space is divided into two parts. 
In the inner region the use of the BF frame impUes the Bom-Oppenheimer ap- 
proximation, which in turn regards the projection of the electron orbital angular 
momentum along the figure axis (i.e. the main symmetry molecular axis) as a 
good quantum number, so that generally either the Hund's case 'a' or 'b' is appli- 
cable to couple the incident electron orbital angular momentum to the molecular 
electronic angular momentum and to the angular momentum of nuclear rotation. 
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Chapter 3 



Results and discussion 

3.1 The ortho-benzyne molecule, 0-C6H4 
3.1.1 Introduction: metastable anions of o-benzyne 

In a carbon-rich red giant or in a proto-planetary nebula (PPN) atmosphere, the 
major carbon-bearing chemical species that is initially formed is acetylene (C2H2), 
rather than the carbon monomer, due in principle to the wide availability of hydro- 
gen. Generally speaking, specific chemical reactions, depending on both density 
and temperature conditions, could then lead to the formation of planar molecules 
of polycyclic aromatic hydrocarbons (PAHs), large aromatic molecules with five 
or six-membered carbonaceous rings. Due to their marked electronic delocaliza- 
tion, they are very stable against photodissociation, thus they could also be present 
in circumstellar environments. In this connection, the still unassigned infrared 
(IR) bands detected in different regions of the interstellar medium (ISM) are com- 
monly related to infrared fluorescent emissions from highly excited PAHs iTTSll . 
Moreover, some of these bands have been observed in planetary nebulae P31l and 
have also been detected in stars that were in transition from the asymptotic red 
giant branch (AGB) to the planetary nebula phase (PN) [28J. It seems reasonable, 
in fact, that PAHs are formed in the envelopes of C-rich stars and within regions 
spectroscopically characterized by the featureless, continuum infrared (IR) emis- 
sion attributed to amorphous carbon and to the formation of C-particulates ll25ll . 
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A similar mechanism may also be at work in circumstellar shells where R^H 
molecules could be viewed as providing the building blocks of stellar dust p2|. 

In line with the above considerations, some chemical schemes assume RA.Hs pro- 
duction in the molecular envelopes of carbon-rich stars where, just assuming the 
acetylene as the major carbon source, thay take advantage of classical nucleation 
theory in order to provide estimates of the formation rates [[39ll . 

More detailed kinetic mechanisms for gas-phase formation of R\Hs were already 
put forward in [26J, whereby the onset of their production in the ISM was pre- 
dicted to occur within the temperature range between 900 and 1 100 K, and require 
particularly dense and slow stellar winds with high initial acetylene abundances. 
Because of the presence of a fair amount of photoionized electrons within that 
stellar wind, it is also a likely possibility that very reactive, metastable negative 
ions are formed under the above conditions and therefore play a significant role 
during the R\Hs formation. 

Unfortunately, although significant laboratory spectroscopic work has been done, 
and dedicated to R\Hs in order to support the astronomical observations, their 
identification has always been elusive. Only benzene has been tentatively detected 
[|22l| towards the proto-planetary nebula CRL 618, since its z/4 bending mode (the 
strongest infrared feature for CqHq) has been observed, and consequently the col- 
umn density as well as the kinetic temperature have been determined to be respec- 
tively of 5x10^^ cm^^ and of 200 K. The discovery of benzene molecules [|22l . 
which had never been found outside the solar sistem, thus represents a funda- 
mental step for our understanding of the possibile interstellar and/or circumstellar 
synthesis of R^Hs because of its being the simplest, basic aromatic unit. More- 
over, the detection of circumstellar benzene molecules indicates that both physical 
and chemical conditions in PPN CRL 618 are right to activate a rich and diverse 
chemistry. This PPN, in fact, is an evolved object with a central BO star 03611 (T^ 
= 30.000 K) and an ultracompact H II region surrounded by a carbon-rich (C/0 
> 1) molecular envelope [f36ll . In addition, it shows optical high-velocity bipolar 
outflows as well as a low velocity expanding torus ll37l . the latter being charac- 
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terized by an intense molecular emission by means of which it has been possibile 
to realize a rich molecular inventory including a large variety of hydrocarbons 
[|37l |6l. It then follows that either the UV photons from the rapidly evolving 
BO central object, or the strong shocks associated with the high-velocity stellar 
winds (up to 200 Km s^^), could in principle drive the polymerization of acety- 
lene and consequently form benzene as the starting point for PAHs' synthesis. In 
this connection, therefore, the circumstellar(Q) chemical formation of large aro- 
matic species has been modelled not only in strong analogy with soot formation 



processes which take place in terrestrial combustion environments (see figure 3.1 



[|26l|23l), but also surmising a complex chain of reactions driven by photoioniza- 
tion and dissociative electron recombination processes (see figure |3.2[ 071). 
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Figure 3.1: The radical-based reaction scheme for PPN benzene chemistry (adapted from 
ll26l ) in which the 0-C6H4 is expressly suggested as an actual competitor to benzene, thus 
leading to the larger PAHs. 

The carbon-rich nature of proto-planetary nebula CRL 618, coupled with the 
CgHg, C6H2, C4H2 detections (221, niakes it an excellent target for additional 
benzene derivative searches, although no spectral features associated with 0-C6H4 
has yet been detected towards this object [|44l . Keeping in mind that we do not 
have as yet reliable predictions of CeHg, CeHs and 0-C6H4 abundances in PPN 
'for the outflows of AGB and post-AGB C-rich stars 



92 



CHAPTER 3. RESULTS AND DISCUSSION 




Figure 3.2: The ion-based reaction scheme for PPN benzene chemistry (adapted from 
BTl ). See text for details. 



(i.e. on the key species involved in the ionic and in the radical-based reaction 
mechanisms as those given in figures 3.1 and 3.2[ ), new observations and better 
theoretical data on both CeHs and 0-C6H4 would yield, in principle, important 
information as well as new insight for future modelling efforts. We therefore 
present here a study directed at providing one of the possible reasons for which 
the ortho-benzyne molecule might not have been detected during searches over 
CRL 618 [ill, despite of its large electric dipole moment of 1.68 D (971 and of 
its lack of hyperfine splitting features, both of which would make it a very good 
candidate for detection with infrared and microwave spectroscopy. In fact, at least 
qualitatively, one could say that a molecular system such as 0-C6H4 is likely to 
be present in the torus of CRL 618, where higher densities could shield it from 
photodissociation as well as photoionization, while at the same time allowing its 
reaction with low-energy thermalized electrons, present in that envirormient [|36l . 
to form the corresponding metastable resonant anion. The present study is fur- 
ther encouraged by the fact that three molecular anions such as C4H , CgH and 
CgH" have recently been detected (and in large quantities) in different astrophys- 
ical environment [|6T1 [H |9l [6OI [59l| like the envelope of the carbon-rich AGB star 
IRC-i-10216 [l6D[8l[60l, the latter being considered the AGB star archetype. 
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The possibility that a not so small but marginal fraction of interstellar molecular 
material might be in the form of anions was first suggested by Herbst [29J in re- 
lation to carbon chains and other radicals with electron affinities large enough to 
exhibit in principle large radiative attachment rate. This characteristic is linked to 
the fact that the higher the electron affinity of a molecule the more states are ac- 
cessible for the redistribution of the excess energy carried by the attached electron 
while, at the same time, auto-ionization is less likely to occur before the compet- 
ing radiative stabilization of the excited anions could take place. 
Thus, the question we wish to address here is the possible existence of metastable 
anions of 0-C6H4 as "doorway states" for the ensuing possible stabilization of 
different negative species which, in turn, can lead to ionic reaction chains on the 
way to the formation of increasingly more complex PAH species. 

3.1.2 Possible chemistry and role of stabilized 0-C6H4 negative 
ions in PPN atmosphere. 

Although significant laboratory spectroscopic work has been dedicated to PAHs 
in order to support the astronomical observations, their identification has always 
been elusive: none of them has been unambiguously detected and only benzene 
have been tentatively revealed toward the proto-planetary nebula (PPN) CRL 618 
[|22l . Hence, the first step to understand the possible involvement of the ortho- 
isomer of benzyne (0-C6H4), and specifically of its negative ion, in a realistic 
circumstellar gas-phase PAH synthesis is to analyze the benzene -related circum- 
stellar chemistry i.e. its currently suggested formation pathways at the physical 
conditions that characterize the proto-planetary nebula stellar phase: in this re- 
gard, the benzene (CeHg) could be considered as the corner-stone of carbonaceous 
poly cyclic aromatic structures. 

Generally speaking, as mentioned in the introduction, PAHs are considered to 
be produced in carbon-rich star outflows, as it seems not very likely that such 
large molecules could be produced by gas-phase or by grain-assisted chemistry 
within the interstellar medium, given the much lower pressure (and density) of 
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that environment. The PPN phase, which for stars of low and intermediate initial 
mass represents the link (the transition object) from the asymptotic giant branch 
(AGB) phase and that of the planetary nebula (PN), is at present both physically 
and chemically rather poorly understood, and even though its evolution is tipically 
fast as compared to the initial AGB and final PN phase lifetimes, it is however an 
established fact that investigating its microscopic (chemical) evolution could pro- 
vide a deeper insight in predicting the nature of the molecular and atomic material 
returned to the ISM as the central star approaches the end of its life. The circum- 
stellar chemistry for a PPN, or a post-AGB atmosphere, appears to be strongly 
dependent on both the radiative flux from the central star and its density, because 
circumstellar material undergoes energetic processing. Various observations [O 
have suggested that during the PPN phase a second, more violent, mass-loss event 
occurs with high velocity winds that partially destroy the remnant AGB enve- 
lope, so that over a relatively short period that envelope is depleted and gradually 
exposes the central star OTll . As the envelope thins down, the stellar radius de- 
creases, raising in turn the temperature of the star to T > 30000 K. Finally, the 
increasing UV output from the central star itself begins to ionize the remnant en- 
velope, thus leading to a PN [|3T1l . 

In this framework PPN CRL 618 represents an extreme example of an evolved 
object, at present characterized by optical high-velocity bipolar outflows, where 
the central BO star, having lost its dusty shroud, is producing sufficient UV radi- 
ation to create a small, compact nearby HII region [|36l surrounded by a carbon 
rich molecular envelope: this slowly expanding toroidal envelope (showing a rich 
molecular emission [|37l[38ll ) is a high-energy, high-density environment so that its 
chemistry could be expected to be similar to the one observed in combustions [|26l . 
Consequently, during the PPN evolution the surrounding gaseous shell expands, 
the density decreases and the material in the circumstellar envelope is subjected 
to photoprocessing from the central star, thus very likely leading to predominant 
ion-molecule chemistry P6l . 

It is for the above reasons that we can surmise as a formation mechanism for large 



3.1. THE ORTHO-BENZYNE MOLECULE, O-C6H4 



95 



molecules the likely combination of radical- and ion-molecule reactions: in fact 
both a radical-based mechanism [|26l l23l and an ion-molecule mechanism [|47l 
have been proposed for the benzene molecule circumstellar synthesis (see figures 
3.1 and 3.2| ). Moreover, it should be noted that while the former Il23l directly 



supports and accounts for the possible presence of ortho-benzyne (neutral ortho- 
benzyne), the original ionic chemical model pTl does not include that formation 
reaction. Such a reaction, however, could be realiably induced by free electrons 
and accompanied by the phenil cation (CgHjl") dehydrogenation: in fact the ionic 
chemical model involves ionized hydrocarbon - neutral acetylene reactions as well 
as hydrogenation and/or (dissociative) electron recombination processes, thereby 
implicitly accounting for the presence of low-energy electrons that we argue here 
to play an important role. Additionally, one of the benzene derivatives formed in 
the radical network, the phenil open-shell radical (CeHs), is a suggested precursor 
to PAHs in circumstellar environments Il25l l23l due to its high reactivity with 
acetylene. However, at present no PPN model includes both the radical- and the 
ion-based benzene reaction networks, thus making difficult the interpretation of 
observational results. 



Some current studies of high-energy environments indicate that, if circumstellar 
chemistry is similar to plasma or combustion chemistry, then CeHs and 0-C6H4 
may coexist with CgHg: theoretical investigations showing that benzene unimolec- 
ular decomposition during combustion actually leads to CgHs and 0-C6H4 pro- 
duction do exist [|33l . as well as experimental studies indicating that gas-phase 
benzene electrical discharges produce CeHs and 0-C6H4 [|32l . However, no spec- 
tral features associated with 0-C6H4 were detected in the MHz frequency range 
during its astronomical search toward CRL618, despite of its rather large perma- 
nent electric dipole moment coupled with the lack of hyperfine splitting. Such 
features should in principle yield stronger lines for the neutral molecule, making 
it a more likely and easier candidate for detection than CgHe, if their abundances 
were comparable. A molecule such as 0-C6H4 is however likely to be present in 
the toroidal molecular envelope of CRL618, where higher densities could shield it 
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from both photodissociation and photoionization. More in general, such a species 
could also be present in more extended regions when actually produced by ion- 
molecule chemistry. 

We therefore follow here the suggestion of p4| and make additional investiga- 
tion on possible chemical networks involving 0-C6H4. As will be shown in what 
follows, from the results involving the o-benzyne molecule at the equilibrium ge- 
ometry [21 J, we learn that this molecule is an efficient producer of compact, fairly 
long-Uved metastable anionic intermediates. In line with that, we then present 
the logical extension of that initial investigation, in keeping also with the occur- 
rence of fast radiationless electronic transitions to an energetically lower-lying 
electronic state for a similar molecule [41 J, and suggest a likely stabilization of 
the metastable negative ion via an intramolecular vibrational redistribution pro- 
cesses, the latter chiefly involving the triple bond within the exagonal aromatic 
structure. This picture, in fact, offers not only a qualitative explanation for the 
surprising absence of spectroscopic features associated with neutral 0-C6H4 to- 
ward CRL618, but enables us to surmise a direct involvement of the stabilized 
negative ion of ortho-benzyne in the complex chain of reactions eventually lead- 
ing to circumstellar R^Hs, as we shall further discuss below. In other words, we 
shall show that the anionic benzyne can provide an highly reactive intermediate 
that is very likely to disappear fairly rapidly via the ionic reactive path which, in 
turn, we shall qualitatively illustrate in the following sections. 

3.1.3 The integral cross sections 

Equilibrium structure, geometry optimizations and the potential energy of the neu- 
tral species involved in this study, as well as the starting molecular orbitals (MOs) 
used for the present scattering calculations, have been performed at the HF/aug- 
cc-pVTZ level using the Gaussian suite of programs ni4ll . These ab initio calcu- 
lations yielded for the title system a total SCF energy of -229.467905 hartrees. 
The molecular structure obtained with the above Gaussian basis-set shows good 
agreement with the earlier more accurate ab-initio calculations ifTSl : see figure 
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33] for their values. Using the above basis set the calculated wave function for the 



neutral target yield a permanent dipole moment of 1 .77 D, in good agreement with 
the theoretical value of 1.68 D [97 1. According to the molecular symmetry point 
group to which the ortho-benzyne belongs at its equilibrium geometry, reported 



by figure 3.3 the elastic cross sections and poles of the S-matrix for the scattering 
process 



o - CeH^C^) + e- ^ (o - CgH^'YC^) 

have been computed within the four irreducible representations of C2V group using 
our scattering program [fT6l: these calculations have been initially performed by 
us at the equilibrium geometry. 

For the S-matrix calculations we expanded the wave function up to a maximum 
angular momentum value of Lmux = 40; the corresponding partial wave expan- 
sion for the scattering potential was thus carried out up to /max = 80; moreover, a 
maximum of 25 partial waves (LmoxA = 25) was used by us for the corresponding 
scattering electron expansion within the piecewise diabatic potential approach. 
We note that to increase either values of the maximum angular momentum (L Max 
and LuaxA respectively up to 60 and 30) had no significant effect on the final 
results. The size of the physical box that contains the diabatized partial terms, as 
described in the previous section, was 50 A, then the overall radial and angular 
grids (r, 6*, 0) included 1760 x 80 x 80 points. The model correlation-polarization 
potential Vcp is the one described in a previous section (2.4.5) and in our earlier 
work lfT0l[TTl[T2l[T3l ; it was used in the present calculations with the total asymp- 
totic polarizability value of 64.98 by us calculated. We shall describe below in 
greater detail how this term was included in our calculations. 
In what follows, we shall describe and discuss our computed elastic (rotationally 
summed) integral cross sections for the equilibrium geometry of the title system. 
The actual energy behaviour of the dominant symmetry components which make 
up the total electron scattering elastic cross-sections is depicted, in the energy 



range of interest, within the three panels of figure 3.4 
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Figure 3.3: Computed optimized equilibrium geometry and bond values ( in A) for the 
ortho-benzyne molecule. See text for further details. 

We report there the partial elastic cross sections for the a2 (top panel), bi (mid- 
dle panel) and 62 (bottom panel) contributions; no maxima in the integral cross 
sections are found in the totally symmetric (ai) component, the latter not being 
shown here. We are initially interested in revealing the possible existence of low- 
energy shape resonances in the 0-C6H4 (^S), features which are qualitatively due 
to short-range effects leading to dynamical trapping of the impinging electron. 
However, the ortho-benzyne molecule, due to its triple bond, has a large perma- 
nent electric dipole moment as well as a large spherical polarizability having their 
origin in the large number of vr electrons in the aromatic pseudo-exagonal ring. 
Thus we cannot in principle discard the long-range polarizability contributions on 
the resonance parameters, and in the present case we indeed expect they shall play 
an important role in the scattering process. 

Our long-range polarization potential contribution to the Vcp potential is obtained 
by using a potential form which asymptotically agrees with the value given by the 
static polarizability a of the molecule under investigation. This can be done either 
by assuming a single polarization center or by partitioning the static polarizability 
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Figure 3.4: Computed elastic integral cross sections of the three molecular symmetry 
contributions. For the definition of the various curves shown see the main text. 
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a on different centers where the individual atomic polarizabilities are subjected 
to the constraint that the total molecular polarizability is reproduced in the long- 
range region. 

We have actually performed both such calculations including the long-range part 
of the electron-molecule interaction for all the elastic partial cross sections in the 
three contributing symmetries and also increasing the physical box up to 50 A. 



The corresponding results are reported in figure 3.4 the continuous lines repre- 
sent the elastic partial integral cross sections computed without using the 0-C6H4 
(^E) polarizability and at the same time referring to a physical box with a radius 
of 15 A, while the circles and the diamonds refer respectively to a calculation 
which is now including either the spherical or the multicenter polarizability (the 
latter weighed and scaled by the MuUiken charges on each atomic center) as well 
as considering the larger physical box of 50 A. We can immediately recognize 
that both the 02 maxima and the higher energy maxima of h\ and 62 contribu- 
tions appear at the energy values for which we also found the maxima without 
the dipole-induced polarization corrections. The situation is instead different for 
each of the very low-energy features visible in the partial cross sections of h\ and 
62 symmetries: when the polarizability is introduced in the calculation, they both 
disappear below threshold thus suggesting that they are actually artefacts created 
by the limited size box employed initially. All further calculations were therefore 
done employing the larger integration "box" mentioned before (50 A). 



3.1.4 The cr* and tt* resonances 

We discuss below our results on the computed resonance parameters, i.e. the 
eigenphase sums and the diabatic potential curves, reported respectively in figures 



3.5 and 3.6 We are in fact interested in characterizing the metastable anions 
obtained from resonant trapping of the impinging low-energy electron because of 
their possible role as reaction intermediates in the circumstellar environment. 
In figure 3.5 we show the behaviour of the eigenphase sums together with their 
first derivatives for each symmetry in which we find a maximum in the corre- 
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spending partial cross-section. The four eigenphase sums were each fitted to a 
Breit-Wigner formula, using the familiar functional form [|40l 



5{E) = a + b{E- Eres) + c{E - Eres) + arctan 



(3.1) 



_2{E — Eres) _ 

in order to extract the pure resonant component from the background and to obtain 
the resonance width from the lorentzian fit of the corresponding delay matrix. It 
is clearly seen that the four isolated resonant contributions in the three different 
symmetries 02, &i and 62, show a sharp tt jump when the collision energy of the 
charged projectile crosses the value for which each of the corresponding cross 
sections shows a maximum: these findings confirm that the four maxima in the 
elastic partial cross sections correspond to pure shape-resonances. The energy 
dependence of the lower diabatic potential components for each resonance are 



shown in the four panels of figure 3.6, where angular momenta are reported up to 
/ = 3 (61 and 62 contributions) as well up to / = 4 (for the two 02 resonances). 
We find it important to point out that the diabatic potential curves of the bi and 
62 symmetries are such that the corresponding impinging electron energies are in 
both cases degenerate with the minimum value needed to overcome their barriers; 
this particular feature is of importance in controlling the width of a resonance. 
The two resonances of 02 symmetry, as clearly shown in the two upper panels of 



figure [X6| are instead due to the tunnelling across a barrier generated by the / = 3 
(/ wave) and I = A (g wave) angular momenta, features that therefore explain 
their being located at higher energies. 

By examining the contour plots of the resonant wavefunctions of the three irre- 
ducible representations (IRs) in which we find shape-resonances (which we do 
not report here for the sake of brevity), we have discovered that three of the four 
resonances have n* nature (respectively the two of 02 symmetry and the one of hi 
symmetry), while the remaining 62 resonance has a* nature. 



The &2(<7*) is now explicitely shown in figure 3.7 and compared there with the 
nearest virtual orbital (MO) obtained from the HE calculations on the target molecule. 
It shows that the distribution of electronic charge associated with the impinging 
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Figure 3.5: Computed eigenphase sums within each of the scattering contributions for 
the three symmetries discussed in the main text. All values are in eV. 




Figure 3.6: Computed diabatic potential energy curves employed for the calculation of 
each partial elastic integral cross section. See main text for details. All energies are in eV. 
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= 4.88 



Figure 3.7: Left-hand panel: shape of the real part of the scattering wave function (w.f.) 
for the a* resonance computed using the piecewise diabatic model interaction (color on- 
line). Right-hand panel: HF/aug-cc-pVTZ computed virtual MO of the ground electronic 
state of the 0-C6H4 molecule for the 62 symmetry where the a* resonance is located. All 
energy values in eV. The lifetime r in units of 10^^^ s. 

electron is displaced asymmetrically, mainly on the carbon atoms of the triple 
bond where a nodal plane clearly exists perpendicular to that bond: this strongly 
suggests its antibonding behaviour across this CC bond. Thus, according to these 
findings, the &2(o"*) resonance at the equilibrium geometry could be reasonably 
considered a good candidate for adding antibonding behaviour to the triple bond 
in the (N-i-1) electron system and therefore for being responsible of a possible 
ring-breaking process guided by the non-adiabatic coupling between this a* reso- 
nant electron and the nuclear dynamics which then would chiefly involve the triple 
bond. The following section shall further elaborate on this point. 
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3.1.5 0-C6H4 versus CqHq: a comparison between their equi- 
librium metastable anions 

From a structural point of view one of the main differences between benzene and 
ortho-benzyne, both considered in their equilibrium geometries, lies in the lack for 
the latter of a vicinal couple of hydrogen atoms to which corresponds the presence 
of the triple bond. The ensuing appearance of a large permanent electric dipole 
moment (1.68 D) in the ortho-benzyne, in spite of the 0-C6H4 and CqHq having 
aromatic rings, turns out to be responsible, at least qualitatively, for many of the 
differences that we shall find below between the behaviour of the resonant states 
of these molecules. 

Let us briefly summarize the properties of the two systems: 

i) the benzene molecule in its equilibrium geometry belongs to the Deh symmetry 
group, while the ortho-benzyne belongs to the C2vgroup; 

ii) benzene is characterized by resonances in seven distinct irreducible represen- 
tations which correspond to symmetries for which there are no occupied MOs 
(«2g» b2g, e2u) as well as to symmetries where occupied MOs of the target ground 
electronic state exist (aig, e2g, &iu, eiu) [1901; 

iii) according to the diabatic potential curves and to the eigenphase sums be- 
haviour as functions of the impinging electron energy, the first group of computed 
resonances [90] can be viewed as actual metastable negative ions of CqHq; 

iv) for the resonances which belong to an IR that also contains occupied MOs 
(with the exception of aig contribution), the corresponding diabatic potential curves 
indeed exhibit their lowest angular components with an / value greater than zero 
and located at rather low energy, thus having a barrier too low to support a shape 
resonance. In fact, each of their cross sections also shows one low-energy struc- 
ture which could be qualitatively related to the top of the centrifugal barriers for 
the lowest contributing angular momenta ||90l . 

In order to make a more efficient comparison between the benzene and the ortho- 
benzyne resonances, we can use the correlation table C2V ^ Deh shown in table 
3.1. A first important difference between these molecules can be found in the 
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fact that no resonance features appear in the totally symmetric contribution (ai) 
for the ortho-benzyne, while for benzene two one-electron resonances, both with 
symmetry aig (for which also bound MOs exist) can be identified. 
Each of the partial cross sections in the bi and 62 symmetries for 0-C6H4 shows, 
as discussed before, that high energy maxima do not undergo variation of their 
energy location once polarization effects are included. Moreover, their eigenphase 
sums as energy functions show a sharp vr jump, thus indicating the dynamical 
trapping of the incident electron due to a nonadiabatic coupling between the / = 2 
and I — 3 angular momentum contributions (see Table 3.2). 



Deh 


C2V 


Aig 


Ai 


A2g 


A2 


Big 


Bi 


B2g 


B2 


Eig 


Bi + B2 


E2g 


Ai + A2 


Alu 


A2 


A2u 


Ai 


Blu 


B2 


B2U 


Bi 


Elu 


B2 + B1 


E2U 


A2 + Ai 



Table 3.1: Correlation table C2v ^ Deh between the symmetry groups associated with 
benzene and benzyne. 

Through the above correlation table is also possible to connect each of the four 
ortho-benzyne resonant states in table 3.2 with those found for the CqHq molecule 
(reported in table 3.3). In general, however, with the exception of the 02 (0-C6H4; 
Bres = 8.0 eV) — > a2g (CqHq; Ej-es = 21 eV) connection, which is geometrically 
rather simple to do, the remaining three metastable states for the ortho-benzyne 
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IR 


Eres{eV) 


T{eV) 


^dominant 


b2 


1.05 


0.27 


2/3 


bi 


2.29 


0.61 


2/3 


a2 


3.33 


0.63 


3 


a2 


8.04 


0.81 


4 



Table 3.2: Computed resonances for the 0-C6H4 target. See text for further details. 

molecule do not appear to be comparable with any of those found for the benzene 
molecule in the calculations of ref [|90l . In other words, the removal of a vicinal 
couple of hydrogen atoms from the aromatic ring and the consequent presence 
of the large permanent electric dipole moment due to the triple bond, make the 
resonant metastable states for the ortho-benzyne very distinctive and it is then 
difficult to compare them with those found for the benzene, as they are in general 
different not only in the geometrical shape of the real parts of the resonant wave 
functions but are also associated with different dynamical mechanisms, as we shall 
further see below. 



IR 


Eres{eV) 


T{eV) 


^dominant 




4.66 


1.50 


3 


b2g 


9.02 


2.26 


4 




12.25 


4.75 


5 


aig 


12.32 


6.78 


4 


e2g 


15.46 


7.04 


4 


blu 


16.52 


12.81 


5 


a2g 


21.49 


4.82 


6 


Clu 


21.96 


14.63 


5 


blu 


24.48 


5.31 


5 


aig 


25.73 


10.09 


6 



Table 3.3: Computed resonance parameters for the benzene molecule as reported in ||90l . 
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There are also additional differences that have to do with the global behaviour of 
the metastable anions o-C6H4^ and CeHg^, both in their ground electronic states, 
as summarized in the tables 3.2 and 3.3. With the exception of the e2u contribu- 
tion, the benzene molecule exhibits resonances located at markedly higher ener- 
gies than its ionization potential (9.24 eV, [fTTl ). Furthermore, the corresponding 
widths allow one to obtain for each benzene resonance a lifetime which is shorter 
than that associated with their nuclear vibration dynamics. The ortho-benzyne, 
on the other hand, appears to be somewhat more "reactive" in forming metastable 
anions: the ones calculated and discussed here correspond to energies which are 
well below its ionization potential (9.02 eV, ifTSll ) and, as reported in the third col- 
umn of table 3.3, the ortho-benzyne resonance widths are smaller, corresponding 
to greater lifetimes which then become comparable with the times for vibrational 
motions within the bound molecule. 

On the basis of the above findings we can therefore suggest that in the energy 
range of a few eV, that correspond to kinetic temperatures up to about 10.000 K for 
the circumstellar atmosphere of CRL 618 [|36ll . the ortho-benzyne molecule could 
actually enter the complex chain of interstellar reactions by forming metastable 
anionic states more efficiently and less transiently than benzene. 

3.1.6 Pathways to ortho-benzyne stabilization: computational 
findings 

In the previous sections [|2T]|. we have investigated the quantum dynamics of low- 
energy electron scattering from gaseous 0-C6H4, at its equilibrium geometry. We 
have shown there that electron-benzyne single collision events give rise to reso- 
nant scattering states, since 0-C6H4 shows four resonances under 10 eV one of 
which is around 1 eV, a fairly long-lived a* resonance with clear antibonding 
character across the triple bond. Such a metastable state, because of its longer 
lifetime and of its relatively small excess energy, has been considered by us a 
likely candidate which could play a crucial role in ring-opening paths, in specific 
bond cleavage (each of the CH bonds or specific CC bonds), or also in radiation- 
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less stabilization by intramolecular vibrational energy redistribution (IVR). By 
assuming the formation of neutral ortho-benzyne in an astrophysical environment 
like a circumstellar atmosphere of a carbon star where there are also free 
electrons as suggested in [|471 . a ring-opening path would provide a negatively 
charged intermediate which could then initiate condensation reactions with the 
acetylene molecules of the same environment; furthermore, on the other hand, a 
stabilization of the whole negative ion would lead in principle to a highly reactive 
anionic species that can in turn react with carbonaceous cations (like C2H^ and 
C4H^), participating in a second aromatic ring closure toward the formation of 
the naphthalene structure. 

In this section we thus want to focus on the investigation of such evolutionary 
channels of the metastable o-CgHj by analysing its behaviour when the nuclear 
structure is deformed: in other words, we want to consider the possible elemen- 
tary process where the excess resonant electron energy could couple with, and 
transfer into, the internal energy of the network of molecular bonds. Whenever 
one specific bond can be reasonably considered as more directly responsible for 
energy transfers to the vibrational degrees of freedom, by selecting different one- 
dimensional (ID) pictures for the electron attachment process one should observe 
a corresponding shift in both the resonance energy position (the real scattering 
energy) and width (the imaginary energy contribution). Since for a polyatomic 
molecule like the 0-C6H4 system the dynamical coupling of the resonant electron 
with its vibrational degrees of freedom is a multidimensional process, to account 
for it as realistically as possible we carried out several scattering calculations us- 
ing partially optimized molecular (nuclear) structures in which various specific 
molecular bonds were selected (and deformed) as cuts of a larger multidimen- 
sional molecular deformation path. 

The calculations previously presented [|2T]| have shown a partial localization of 
the extra electron as well as a clear presence of nodal planes across the triple and 
the single bond respectively. Thus, one possible TNI evolution is given by the 
ring-breaking pathways which, in a simpler ID model, can be seen as lengthening 
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of either the triple CC bond or the single CC bond on the opposite side of the 
molecule. A different pathway to stabilization is represented by the cleavage of 
each CH bond located respectively in the ol and /3 positions with respect to the 
triple bond, a process that can cause the loss of neutral hydrogen atom (forming a 
highly reactive closed-shell radical anion CgHg ) or the detachment of a negative 
atomic ion (leading then to a highly reactive radical CgHs). We have therefore 
computed the corresponding pseudo-ID potential energy curves associated with 
both the neutral N-electron target and the metastable (A^ + l)-electron negative 



ion (figures 3.8 and 3.9 1, setting the energy balance by the following formula: 



-E'Tot(R.) = ERes(R.) + En(R) — EN(Req). (3.2) 

In this expression Ej^i^^^ is the electronic energy of the neutral parent molecule 
at its equilibrium geometry, computed at HE level with the aug-c-pVTZ quality 
of the basis set, E^ilVj is the computed electronic energy (at the same level of ac- 
curacy as for the equilibrium geometry) for the neutral deformed molecule at a set 
of geometries identified by the involved single -bond coordinate R, and Ejies(R) 
is the real part of the computed resonant electron energy over the same range of 
molecular geometries. The corresponding widths r(R) of each resonance, associ- 
ated with the (A^ + l)-electron states, are related to the complex scattering energy 
via the well known relation [|40|| 



E^°r'^^(R) = ERes(R) + irRes(R). (3.3) 

According to the energy-time uncertainty, the above widths are then indicative of 
the TNI lifetime at each computed molecular geometry. 



We report in each panel of figure 3.8 the real part of our computed ID potential 
energy curves when the single CC bond and each of the CH molecular bonds are 
deformed. The black solid lines identify the Bom-Oppenheimer potential energy 
curves for the neutral parent molecule, while the red ones are the real branches 
of potential energy for the metastable anion; the vertical bars refer to the imagi- 
nary resonant energy component (r/jes(R)) for the deformed nuclear structures. 
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For each of the pathways reported in figure |3.8[ we can easily recognize a non- 
fragmenting behaviour for our metastable state: generally speaking, in fact, de- 
pending on the shape of the potential energy curve for the metastable anion the 
TNI itself can evolve either toward bond fragmentation or anion stabilization. 

From the inset in the upper panel, containing the computed real and imaginary 
(vertical bars) resonance energies as a function of CC single bond deformation, 
we see that the extra electron is relatively weakly coupled with this nuclear mo- 
tion and that it remains attached to the target molecule since the lifetime increases 
as the CC bond is lengthened. A very similar situation characterizes the CH bond 
deformations located in a and /3 positions, for which the extra electron also ap- 
pears to be uncoupled with these ID nuclear motions. In sum, our analysis of each 
of the above deformations suggests that these evolutions of the nuclear dynamics 
yield a bound metastable negative ion where the extra electron seems to remain 
attached to the target molecule for at least one vibration. 



In the case of the triple CC bond stretching (see figure 3.9 1, the resonant electron 
is seen instead to be strongly coupled with this nuclear motion, as suggested by 
the steep downward behaviour of the real scattering energies when the bond is 
deformed (upper panel) moving from a rather compressed geometry to a slightly 
stretched bond with respect the equilibrium value for the neutral. 

Furthermore, the resonance disappears into a bound state for a deformation of ~ 
0.1 A beyond its equilibrium value: in particular we are able to locate the last 
resonance as a pole in the S-matrix. Our findings, summarized in the lower panel 
of the same figure, indicate that a relatively small lengthening of the triple bond 
causes the potential energy curve associated with the metastable negative ion to 
cross the one of the neutral, and to move toward the stable anion. To support this 
occurrence, accurate ab-initio calculations on the stable negative ion (0-C6H4 ) 
indicate that the triple bond has to be stretched with respect to that of the neutral 
[l35ll . We therefore report in figure [3TT0l a comparison between the shape of the 



resonant wave function (on the right) associated with the last TNI geometry we 
found during C=C bond lengthening, and the highest occupied molecular orbital 
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Figure 3.8: Upper panel: behaviour of the real part of potential energy curve for the 
resonant ct*(62) anion (o-C6H4)~ (upper red curve) and for the neutral species 0-C6H4 
as a function of the single CC bond deformations (lower black curve). The small inset 
shows the computed changes of the resonance features (Ejies and T/jes) as a function of 
the same single CC bond lenght. Central panel: the same as the upper panel, but for the 
CH bond located in a position respect the triple CC bond. Lower panel: the same as the 
upper panel, for the /3 CH molecular bond. See main text for details. 
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Figure 3.9: Upper panel: computed variations in the resonant features {Er^s and V^es) 
as the triple CC bond is deformed. The inset shows an enlargement around the equilibrium 
value for the triple CC bond. Lower panel: shape of our computed potential energy curves 
for the resonant anionic species (upper red curve) and for the corresponding neutral parent 
molecule (lower black curve) as functions of the triple CC bond deformations. See main 
text for details. 
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Figure 3.10: Right-hand panel: shape of the real part of the resonant wave function for 
the (7* (62) resonance associated with the +0.1 A lengthening for the triple bond. Left- 
hand panel: DFT/cc-pVTZ computed HOMO of the ground doublet electronic state of 
the stable (o-C6H4)^ anion at its equilibrium. See main text for details. 

(HOMO, on the left) that corresponds to the stable negative ion minimum energy 
structure: the figure clearly shows there is no difference between the two electron 
density maps. 

We are thus encouraged to surmise that the residual amount of energy nedeed 
to form the stable anion (0-C6H4 has a positive electron affinity (EA) of about 
0.56 eV, [jlj llOOID could actually be dissipated via fast radiationless IVR. Such a 
stabilization mechanism is also indirectly supported by the calculations of such 
an internal conversion time involving the S2->Si(~20 fs) singlet electronic states 
for the pyrazine iHTTl : the pyrazine molecule has in fact the same number of atoms 
as in ortho benzyne, thus a very similar density of vibrational modes likely to be 
coupled with the excess electron during an IVR process. 

3.1.7 Estimating the rate coefficients to electron attachment 

The intense radiative flux from the central star is chiefly responsible for the ioniza- 
tion of the inner edge of the dense molecular envelope in CRL618 and therefore 
the energy distribution of the electrons created by photoionization is linked to the 
energy of the stellar photons. 
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In principle, low-energy electrons are more likely to recombine with protons (if 
available) and with molecular cations and thus could be selectively removed from 
the free electron pool. However, electrons thermalize very quickly because of the 
very large electron-electron collision cross sections, so finally their recombina- 
tion rates with molecules or atoms could be considered slow in comparison with 
the interaction among electrons, thereby allowing a local thermodyamic equilib- 
rium (LTE) condition as a good macroscopic approximation for the free electron 
ensemble. One therefore requires a single parameter (the free electron kinetic 
temperature, Tg), to characterize the free electron energy distribution, reasonably 
described by a Maxwell-Boltzmann distribution. In this framework, the rate co- 
efficient for the formation reaction of the metastable ortho-benzyne anion by free 
electron capture 

CqH^ + e- ^ CeH^ (3.4) 

could be estimated by the convolution ||24l between our computed elastic inte- 
gral cross sections and two Maxwell-Boltzmann energy distribution functions, 
one referred to Te < 10000 K Il36ll and the other describing the molecular ki- 
netic temperature, Tm = 200 K [[22|. Such calculation, for the above values of 
Te and Tm, yields a value of kf = 4.2 x 10^^ cm'^ sec^^. At this point, we find 
crucial to emphasize the fact that recently the elementary gas-phase reaction of 
Dl-ethynyl radical (C2D) with vinylacetylene (C4H4) was studied under single 
collision conditions via crossed molecular beam experiments and electronic struc- 
ture calculations Ol. The corresponding findings point toward the formation of 
the neutral 0-C6H4 which is in fact identified between all the products and for 
which the computations predict branching-ratios of up to 10 per cent |l3l; the au- 
thors consequently argue that since the reaction path has no barrier, it might be 
accessed even in a cold molecular cloud like TMC-1, where the environmental 
conditions are deeply different than those of a circumstellar molecular envelope. 
In this connection, we find intriguing to also estimate the metastable 0-C6H4 anion 
formation rate coefficient at those environmental conditions, therefore assuming 
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a much lower kinetic electron temperature, Te = 30 K; in this case the rate coeffi- 
cient appears to be increased of one order of magnitude, being kf = 1.62 x 10~* 
cm^ sec^^ 

Since we use here the elastic integral cross sections for the resonant scattering 
event, these findings represent an upper limit for the actual rate coefficient k; 
it is however interesting to see that such calculated estimates suggest fairly large 
rates originating from electron-molecule interactions [|29l . Hence the combination 
of the above values for the metastable 0-C6H4 formation rate coefficient with 
the possible fast radiationless stabilization of such a resonant anion, enables us 
to suggest that the o-benzyne (meta)stable anion might actually play a role in 
the circumstellar chemistry of CRL618 and more in general in the chemistry of 
either PPN atmospheres or cold molecular clouds, if the neutral parent molecule 
is actually produced in those environments. 



3.1.8 The stable o-benzyne anion: a quantum chemistry view 

To further strengthen the computational exploration on the possible effects of 
ring-deformations and bond-deformations on the fate of the benzyne anion, we 
have carried our quantum structural calculations on the stable negative ion of 
o-benzyne. The quantum calculations employed a Density-Functional- Theory 
(DFT) approach [|43]| using a basis-set expansion defined by the well known acronym 
B3LYP/cc-pVDZ ffT4l . generating all anionic wavefunctions as spin-polarized 
doublets; all "passive geometries" were optimized at each point of the deforma- 
tion curves described below. Energy convergence was tested within 10^^ hartrees 
and the optimization convergences were based on the maximum and root mean 
square gradients, with accuracies of about 10~^ hartrees p3|. 



The data of figure 3.11 report the energy changes upon the stretching, from the 
minimum geometry of the anion, of four different CC bonds, labelled according 
to the figure. One clearly sees there that all stretching deformations present rather 
substantial energy barrier when (o-C6H4)~ is deformed away its equilibrium ge- 
ometry. Hence, it is reasonable to argue that electron attachment stabilization 
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Figure 3.11: Computed reative energies changes of o-benzyne anion by stretching dif- 
ferent CC bonds of the molecular ring. See main text for details. 

would not be able to cause ring-opening processes as direct consequences of anion 
formation: the stable molecular anion, in fact, would need a substancial amount 
of excess energy in order to undergo such ring-breaking processes. In conclusion, 
our suggestion that electron attachment stabilization paths would still leave un- 
broken the molecular structure is further supported by the model calculations of 



figure 3.11 



It is also interesting to note here that the stable anion of the minimum structures of 



figure 3.11 is also an even more markedly polar molecule: the permanent dipole 
of o-C6H4^, in fact, is computed here to be around 4.25 D, a value which is nearly 
three times that of the stable neutral molecule of o-benzyne. 



3.1.9 Present conclusions 

The present theoretical investigation has been directed at the computational analy- 
sis of the possible role which could be played by the ortho-benzyne molecule as a 
reactive intermediate in proto-planetary atmospheres. In particular, we have stud- 
ied the quantum dynamics of the scattering of low-energy (< 10 eV) electrons off 



3.1. THE ORTHO-BENZYNE MOLECULE, O-C6H4 



117 



gaseous 0-C6H4 in order to verify its feasibility for forming long-lived negative 
ions as "doorway states" which, after dissipating the extra energy content toward 
the thermodynamic ally stable negative ion, can evolve through chemical reactions 
participating to the complex reactions chain responsible of the formation of PAHs 
species. 

We find it important to emphasize that the present theoretical investigation on this 
particular aromatic system was encouraged by the fact that 0-C6H4 molecule was 
already conjectured in previous experimental work [|26l l23l [32ll to be a possible 
intermediate during the benzene formation in several astrophysical environments. 
In this connection, therefore, in the present work we have endevoured in order 
to tentatively analyse the possibility of an actual negative ion chemistry involve- 
ment in the astrophysical PAHs synthesis, where the anionic o-benzyne, due to the 
positive EA [|2l II 001 that characterizes the parent neutral species, could be con- 
sequently viewed as a reasonable important aromatic precursor of larger species 
belonging to this family. 

All the calculations have been carried out by treating the scattering events within 
a time-independent formulation of the problem, solving the ensuing multichannel 
equations within a SCE description of the scattering states and via a model, local 
potential for the interaction between the low-energy electrons and the benzyne 
molecules, as described with sufficient detail in chapter 2. 

When considering the 0-C6H4 as frozen at its equilibrium geometry, the results 
indicate that this molecule indeed exhibits four fairly narrow resonances below 
10 eV, three of which can be classified as tt* resonances while one of them is a 
low-energy a* resonance exhibiting antibonding character across the triple bond. 
Hence, we single out this resonance as the most likely dynamical mediator state 
which could in principle lead to ring-opening processes. A comparison with sim- 
ilar scattering resonances which involve the benzene species (CeHg), indicates in- 
stead that the latter shows only fairly broad metastable states which start to appear 
above about 4.66 eV. Our conclusions, therefore, rest also on having observed 
marked differences of behaviour between the more stable resonant states of o- 
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benzyne and those appearing at higher energy in benzene, a species considered to 
be the comer-stone step for the chain of reactions leading to PAH formation. As 
discussed in the previous sections, however, this comparison should be grasped 
only in the qualitative sense. In this connection, in fact, one should note that 
computational parameters like the number of partial wave for the incoming elec- 
tron description and the number of interaction potential terms employed for the 
benzene and o-benzyne investigations are little different; due to the fact that the 
latter has been studied with the same approach but several years later, nowadays a 
larger partial wave expansion like that used for the o-benzyne can be easily used 
since it is not so much demanding from the computational point of view. More- 
over, we find it important to emphasize that the present calculations as those on 
the benzene molecule employ a realistic, but approximate, form of interaction, 
and therefore possible experimental resonances would very likely be shifted in 
energy with respect to our estimates. However, since for both the above aromatic 
molecules we have used the same theoretical approximated approach, according 
to the last consideration, we expect that the use in both cases of more realistic, 
improved model potentials would not modify our comparative conclusions. On 
the other hand, we know that the benzene molecule has a large negative electron 
affinity (-1.12±0.003 eV, flU), so we think that even surmising the formation of 
metastable negative ions for this molecule by low-energy electron resonant attach- 
ment, it seems unlikely their actual involvement in such an astrochemical reaction 
network. 

In the framework of the above qualitative comparison, we can therefore collect 
the following useful information. Since the excess energy carried by the extra 
electron has to come from the environmental medium in the atmosphere, the at- 
tachment processes to benzene do require much higher electron temperatures than 
those needed for the formation of metastable anions of benzyne. Furthermore, the 
(o-C6H4~)* temporary negative ions are seen to have markedly longer lifetimes, 
thereby allowing more easily the dynamical coupling of the extra electron with 
the nuclear degrees of freedom which play a crucial role in the ensuing possi- 
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ble molecular break-up or also in possible non-dissociative stabilization paths. It 
therefore follows that our present calculations for the equilibrium structure (sec- 
tions 3.1.3, 3.1.4, 3.1.5) support the picture of an intermediate anionic formation 
for the o-benzyne species that could then rapidly decay either by ring-opening re- 
actions, a feature which therefore could prevent the actual detection of the parent 
neutral species by observational spectroscopy [|44|. or by non-dissociative stabi- 
lization mainly due to intramolecular vibrational redistribution processes. Hence, 
we surmise that the possible presence of neutral 0-C6H4 in proto-planetary atmo- 
spheres like CRL-618, where free low-energy thermalized electrons might be also 
present, could help in producing the negative, metastable species that in turn could 
then be likely to contribute to the chain of reactions relevant for RA.Hs formation. 
Admitting in fact the above hypothesis, we find reasonable that such reactions 
could thus be easily triggered also by the anionic o-benzyne molecule than by 
benzene itself, a feature that could then explain the observation of the latter neu- 
tral species [|22l with respect to the absence of the former neutral one [|44ll . At 
this point it is useful to emphasize that, currently, two different mechanisms for 
the benzene synthesis in proto-planetary stellar atmospheres were developed, one 
involving neutral and cations [I46II and the other being radical-based [|26]| . where 
only the latter directly accounts for the possible presence and production of neu- 
tral o-benzyne; in connection with this, as already illustrated in the introduction 
(sections 3.1.1 and 3.1.2), since the above mechanisms are associated with rather 
different physical conditions, the formation of anionic o-benzyne could be viewd 
as a first step in drawing a new and more complex reaction network for ben- 
zene and more in general for RA.Hs formation in proto-planetary atmospheres: no 
chemical models, in fact, include both the radical- and the ion-based chemical 
networks, so that the lack of a single model makes the observational results quite 
difficult. 

After having shown that metastable anionic formation for 0-C6H4 is a likely op- 
tion in the circumstellar environment, a feature that would therefore suggest its 
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involvement in R^Hs formation in C-rich proto-planetary atmospheres like CRL- 
618, we consequently started to investigate further aspects of this problem as dis- 
cussed in sections 3.1.6, 3.1.7 3.1.8, where our findings about the possible dynam- 
ical evolution once the metastable anion is formed are illustrated and discussed 
with full particulars. 

The pseudo-ID quantum dynamics which we have employed to study the above 
process, as discussed in section 3.1.6, has indeed shown: 

(i) the formation of a long lived metastable anion at fairly low electron energies, 

(ii) the stabilization of the initial molecule by small ring deformations which can- 
not lead to ring-breaking effects, and 

(iii) the formation of a bound anion by small direct stretching of the triple bond. 
The last step of the quantum dynamics futher indicates the possibility of reaching 
the bound (o-C6H4)^ ground state by fast internal rearrangement of vibrational 
states of the benzyne anion. The above conjecture, proven to be realistic by our 
quantum scattering calculations, could then allow the very reactive, anionic radi- 
cal to participate in reactions with H and/or if +, and also in possible condensa- 



tion reactions with acetylene and acetylene -like cationic species (see figure |3.12[ ) 
or even with the corresponding neutrals of the same molecules, hence explaining 
its negligible presence within column density measurements. 

are 



The cationic molecular partner surmised in the coarse outline of figure 3.12 



two of the most likely abundant carbonaceous species in the proto-planetary at- 
mosphere of CRL-618 [47 1 . One should note, in fact, that A+ -i- B neutralization 
reactions have been estimated to have very large reaction rates at the conditions 
of interest {k ~ 10^'' - 10^^ cm'^ sec^^ [29], and therefore could be very efficient 
initiators of ring-addition reactions as those outlined above. 

An important observation to bring to a conclusion the present section. 

As introduced and qualitatively discussed in chapter 2 (sec. 2.2), the FN nuclei 
approximation for polar molecule like the o-benzyne molecule, is not the best one. 
In the framework of this approximation, in fact, it can be shown that the ICS di- 
verges logarithmically due to the slow convergence in the £/ matrix elements 
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Figure 3.12: Surmised schematics of the gas-phase ionic reactions involving the benzyne 
anionic intermediate and either ionized/neutral hydrogen or polyynic cations. 

which dramatically affects the large i,^' range (where the static dipole interac- 
tion dominates). More precisely, for polar molecules the divergence of the elastic 
integral cross section for vanishing energies comes from the divergence of the dif- 
ferential cross section in the the forward direction, which thus in turn causes the 
ICS to become infiniteUJ. 

From a practical point of view, this means that one should be aware of the fact that 
the larger / value, the more trustworthy the computed ICSs are. 

In the present investigation we have used / = 40 as the maximum partial waves 
number, which corresponds to 441 coupled channels considered for the ICS cal- 
culations, where the / = 60 value used for the necessary convergence tests had 
no significant effect on locating the resonances (see sec. 3.1.3). Therefore, since 
the rate coefficient for the TNI (o-CgH^ )* formation is numerically estimated as 
the integral of our computed elastic ICSs weighted by the Maxwell-Boltzmann 
distribution function 
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/ 8 K T 1 

k{Te) = \ / a{E)Ee^dE, (3.5) 

our rate coefficient estimates appear to be reasonably large and physically mean- 
ingful already when referring to ICSs computed with an / = 40 partial wave expan- 
sion. Furthermore, although the logarithmic divergent behaviour in the ai partial 
contribution to the elastic ICS, due to the Boltzmann exponential, at most we 
expect that using a greater partial waves number for the ICSs calculation could 
then give rise to little larger rate constant values, therefore again supporting our 
conclusions. 

In any case, we have planned for the near future to calculate the ICS (and also 
the momentum-transfer cross sections, MTCS) for the present molecule using an 
improved approach which was already used in the past by our group [I104L 
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3.2 The coronene molecule, C24H12 

3.2.1 Introduction: PAH anions in the interstellar medium. 
Role of C24H12 anion 

Polycyclic aromatic hydrocarbons (PAHs) are a family of planar molecules, con- 
sisting of carbon atoms arranged in a typical honeycombed lattice structure of 
fused six-membered rings, typically characterized at the external edges by C-H 
sigma bonds. 

Due to their four valence electrons, each carbon atom in polycyclic aromatic struc- 
ture forms three covalent o bonds with C (and H) atoms, thus usually leading to 
a planar structure; the remaining electron, located perpendicularly to the molecu- 
lar plane in an hydrogen-like p-orbital, overlaps with the other adiacent electrons 
(one per C atom) to form tt molecular orbitals, which can be pictorically rep- 
resented as diffuse delocalized electronic clouds above and below the molecular 
plane, respectively. This peculiar electronic structure of PAHs enables to intro- 
duce the general rule according which, qualitatively speaking, the more exten- 
sive the delocalized electronic cloud, the more stable the molecule. Moreover, 
because of the typically large binding energy of aromatic carbon atoms in their 
planar exagonal lattice, neutral and singly charged cationic PAHs are known to 
be particularly stable among the family of organic compounds: they need, in fact, 
large amount of energy to be photodissociated i^), so that they are currently 
considered to be potential carriers for at least a subset of the absorption diffuse 
interstellar bands (DIBs, UllOi 18511 ). as well as reliable candidates to efficiently 
convert the absorbed energy in infrared emission (the so-called unidentified in- 
frared bands, UIBs [|49ll20ll93ll79ll). 

Besides the fact that the history of PAHs in the Universe has been a real fascinating 
saga, nowadays these organic molecules are considered to be almost ubiquitous in 
the cosmos: free gas-phase PAH molecules in different charge and hydrogenation 
states, in fact, are commonly considered to be an important component of the in- 
terstellar medium (ISM) [gHl HH |92l, of diffuse molecular clouds [80], of dense 
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molecular cores [I108II . as well as of circumstellar environments (like the plane- 
tary nebula NGC 7027, flM ). Unfortunately, although significant and intensive 
astronomical searches have been dedicated to R\Hs, their detection has always 
been elusive. Only tentative identifications of neutral benzene [|22l . pyrene and 
anthracene II107II were proposed in the last years; the first was detected toward the 
protoplanetary nebula CRL 618 since its z/4 bending mode was observed, while the 
last two species (containing respectively three and four fused aromatic exagonal 
C-rings) are considered to be likely responsible for the blue fluorescence observed 
in the Red Rectangle Nebula. 

Through astronomical observations, laboratory spectroscopy, computational char- 
acterization as well as increasingly detailed modeling of several and different as- 
trophysical contexts, during the last 25 years it was possible to unfold their multi- 
ple role in the chemistry of the interstellar medium (ISM), and more interestingly, 
to surmise that most of interstellar processes involving PAHs requires a detailed 
knowledge of their charge distribution. 

One of the first important attempts to model the chemical evolution of PAHs and 
the modifications to interstellar chemistry that consequently would result from the 
presence of such 'large' molecules in the interstellar gas, was provided by Lepp & 
Dalgarno ([[99l). In such a study, the authors concluded that for PAHs consisting 
of molecules with 30-50 carbon atoms(|3, with fractional abundances greater that 
10~®, the corresponding negatively charged species, PAH^, would be the main 
source([^ of negative charge. 

It has then been tentatively shown that, both in dense HlOSi |5T1 [52ll and in diffuse 
[|80ll clouds, electron attachment to neutral PAHs and to nanoscopic grains can 
have significant effects not only on the free electron density: the formation and the 
possible presence of PAH negative ions in those environments have been in fact 
determined to be particularly important in relation to the ensuing recombinations 
involving positive ions, both atomic and molecular. Conversely to what usually 
occurs for 'direct' recombinations between free electrons and positive ions, one 

^thus little larger than coronene 
■^in the sense of 'carriers' 



3.2. THE CORONENE MOLECULE, C24H 



125 



should in fact note that, in general, the recombination processes involving nega- 
tively charged PAHs are less exothermic, due to the PAH's electron affinity iTTlll . 
In this connection, when considering recombinations between molecular cations 
and negatively charged PAHs, the latter could be viewed as actual 'energetic shock 
absorbers' in the sense that some of the recombination exothermicity can be effi- 
ciently stored in the internal energy of the anionic PAH species itself (due to their 
typically large number of vibrational degrees of freedom), so that less is available 
to dissociate the positive molecular ion upon neutralization. On the other hand, 
specifically in relation to positive atomic ions(|^, negatively charged PAHs could 
instead act like 'catalysts' in the neutralization reactions: beside the dielectronic 
recombination, the recombination processes of atomic cations with free electrons 
present in the surrounding gaseous medium, in fact, are mediated by the emission 
of a photon, so that they are intrinsically very slow processes. 

Accordingly, the presence of PAH anions in photodissociation regions (PDR) in- 
fluences the CO transitions [l56ll . as investigated by [1791 . 

The presence of PAH negative ions has also repercussions on the ionic (atomic and 
molecular) abundances for other important and widely diffused interstellar species 
like H+, H^ and HCO+: the importance rest on the fact that the abundances of 
these ions directly affect the ortho:para H2 ratio, which in turn has very important 
repercussions on the ensuing global chemical evolution ( [|53]| ). since for example 
it plays a crucial role either in the kinetics of specific excitation processes (like the 
reaction between ortho/para molecular hydrogen and sulfur dioxide, [|54ll ). or in 
establishing the ortho:para ratios for other species, including the deuterated forms 
ofH+ ([53J). 

As furthermore suggested by recent theoretical and laboratory studies ( 118211711 ). 
aromatic species like PAHs show strong UV absorption in two broad bands cen- 
tered at 200 nm (tt*— t-tt excitation, ~6 eV) and 70 nm (a*— t-ct excitation, ~17.7 
eV) respectively, thus they are considered to be good candidates to explain the 
217 nm (5.7 eV) bump as well as the far-UV rise of the interstellar extinction 



as distinct from the molecular cations 
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curve (ISEC, [|73l). Moreover, since R\Hs absorb efficiently dissociating pho- 
tons (A < 100 nm, E > 12.4 eV), they are considered to be responsible (at least 
partially) in controlling the atomic fraction obtainable by photodissociations and, 
in this sense, the global degree of excitation of interstellar/circumstellar clouds. 
This UV screen could be also increased if FAHs are more abundant at cloud edges 
(limb brightening effect, [[66ll671[68l[63l70l). 

In this framework, the role of the charge state of R\Hs in ultraviolet extinction 
has been the object of some recent theoretical investigations II82II . More in detail, 
the authors of [|82ll use a specific extinction model [|83l l94l where the PAH contri- 
bution to extinction [|94ll is represented by two Lorentz profiles, one for the vr— j-tt* 
electronic transition accounting for the bump feature and the other for the a^a* 
plasmon resonance whose low-energy tail contributes to the far-UV rising feature 
in the interstellar extinction curve (ISEC). 

The crucial point of that work is that it shows how a combination of classical dust 
particles and mixtures of real PAHs satisfactorily matches the observed ISECs, 
where variations of the spectral properties of PAHs in different charge states (in- 
cluding the anions), produce changes consistent with the varying relative strengths 
of both the bump and the non-linear far-UV rise features [|82ll94ll . 
Recently, the role of neutral and anionic PAHs has been also the object of some 
investigations focused on the fractional ionization in dark molecular clouds [[50l . 
In this connection, even in absence of intense external sources of radiation and of 
the background interstellar field, thus admitting a very efficient screening by the 
surrounding PAHs(|^, a not negligible amount of high-energy free electrons and 
UV ionizing radiation can be 'locally' produced due to the interactions of high- 
energy cosmic rays (CR) mainly with H2 [|57l|: it was in fact demonstrated as a 
very likely possibility that the following ionization reactions 



H2 + CR^ + e 



(3.6) 



H2 + CR^ + H + e 



(3.7) 



^and, of course, by H2 and by the solid phase represented by the dust particles 
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do produce highly energetic secondary electrons ([|58l. with mean kinetic energy 
of about 30 eV), which in turn are able to cause the excitation of the Lyman, 
Werner and the other bands of the H2 molecules present in the same medium, so 
that in the ensuing energetic cascade toward the X^S+ ground electronic state for 
the molecular hydrogen, the produced UV photons can be finally able to either ion- 
ize or dissociate other molecular species. Despite the apparent complexity, which 
depends on the abundance ratios between molecular and atomic (neutral and ion- 
ized) hydrogen, on PAH abundances and sizes, on the dust phase represented by 
PAH clusters and carbonaceous nanoparticles, as well as on the other molecular 
species (both neutral or ionized) present in the same environment, the degree of 
ionization for a dark molecular region like TMC- 1 has been tentatively estimated. 
Such a complex mechanism, in fact, has been studied in connection with the de- 
termination of the C„H^/C„H ratios, with n = 6, 4, (8), because all of these anions 
have been recently discovered (CeH^ and CgH" in TMC- 1 and in the C-rich dense 
envelope of IRC-i-10216 [I8l|59l|60l, while C4H- only in IRC-i-10216, UMl): even 
if currently there are no direct observations of negatively charged PAHs in the in- 
terstellar space, since the electron affinities of both C4H and CeH f6T\ are larger 
than those of some medium-size PAHs [|63 J including the C24H12 whose EA is ~ 
0.5 eV [|84l|89l, the above chemical model [50] has been recently investigated in 
order to take in consideration also the possible electron transfer from PAH to the 
above neutral linear polyynes for (n = 6, 4). The results have led the authors to 
surmise that, for such an astrophysical environment, the key factors in determin- 
ing the C„H"/C„H ratios (n = 6, 4) are the free electron density and the density 
of atomic hydrogen: a high fractional abundance of PAHs, i.e. a low fractional 
abundance of free electrons, yields in fact a qualitatively rather good agreement 
with the observations. 

When specifically considering dense interstellar environments like dense molecu- 
lar clouds, several observations suggest that the abundances of smaller PAHs (N^ 
< 30) are decreased |l66l[67l|68l. According to the absence of UV photons, which 
in fact are thought to be not able to efficiently penetrate those regions, Rapaci- 
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oli and his coworkers [|68l US suggested that small-size gas-phase RA.Hs could 
form larger R^Hs just in such dense regions; in line with that, therefore, smaller 
free-flying R^.Hs observed at the border of illuminated clouds may be formed by 
photodissociation of these larger species, provided an intense UV source exists in 
the surroundings. Despite the fact that many spectroscopic observational indica- 
tors about the R^Hs presence almost all over the ISM do exist, no specific types 
of R^H responsible of the previously mentioned features (UIBs, DIBs) has been 
unambiguously identified. The form and intensity of the spectra, in fact, strongly 
vary from one region to another [[651 . hence suggesting that differences in sizes 
and abundances are present; moreover, one should also note that the actual exis- 
tence of R^Hs in dark molecular clouds is difficult to be established observation- 
ally, since their fractional abundance is surmised to depend on the extent of their 
accretion onto the nanoscopic solid grain particles [i50|[ . As a consequence, the 
aggregation process responsible for cluster and nanoscopic carbonaceous grains 
formation in the above regions is still poorly understood, and the abundance of 
smaller gas-phase FAHs remaining screened by UV radiation in such dense re- 
gions, as well as the time needed for the aggregation process to occur, remain 
currently undetermined. Furthermore, since the general constraints on the R\H 
size and abundance, both as functions of the physical conditions typical of dense 
interstellar regions, are currently weak, it therefore follows that most of the de- 
tailed gas-phase chemical models of dense molecular clouds do not expressly 
include R\Hs in the sense that R\Hs are not treated as individual species. In this 
connection, the recent work of Wakelam and Herbst [[108[[ constitutes an excep- 
tion: in such an interesting work, in fact, in order to analyse the role of FAHs in 
dense cloud chemistry (like TMC-1 or L134N), the authors explicitly include free 
medium-size R^Hs with a mean radius of 4A (Nc = 30) and a fractional abun- 
dance with respect uh of 3 x 10^''; due to the fact that R^Hs may not be free in 
dense cold environments, they also include R^H clusters of radius 10 A (Nc ~ 
470) and 100 A (N^ ~ 4.7 x 10^) having proportionally lower fractional abun- 
dances (2 X 10"^ and 2 x 10"^^, respectively). In order to investigate the possible 
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chemical scenario at a fixed temperature and density, the authors of the above 
work include electron attachment processes to neutral PAHs, recombinations of 
PAH anions with smaller positive ions and photodetachment processes of nega- 
tive PAH ions triggered by UV radiation. The PAH anions formation is treated as 
a pure radiative electron attachment 

PAH + e- ^ PAH- + hu 

whose rate coefficient, according to [|76l . is surmised to be fcg = 1-3 x 10~^ cm^/s 
for the attachment of electrons to neutral gas-phase PAHs containing 30 carbon 
atoms. We find it important to emphasize that the authors themselves, originally, 
pointed out that their results depend on the above value of ke only if it is decreased 
by at least two order of magnitude. 

So far, I have summarized some between the recent and most important mod- 
els concerning neutral and negatively charged PAHs in order to provide a short 
survey(|^ about their physical and chemical role in different astrophysical envi- 
ronments. 

Currently, in fact, it is a very demanding task to realistically include PAHs in the 
chemical modeling of different regions of ISM, including dense molecular clouds 
and diffuse clouds. Due to the fact that PAH molecules occupy a middle ground 
in size between small gas-phase species and dust particles, they are consequently 
seen as an 'intermediate' phase among the molecular gas components and the solid 
dust phase of the interstellar matter. It therefore follows that in order to include 
PAHs properly, as suggested in [I108II . it should be necessary to look at them as 
individual species, but the single aromatic species that are surmised to make up 
the PAH content of interstellar clouds currently are not constrained enough from 
spectral observations and laboratory measurements. Consequently, as emphasized 
by the authors of I1108II , the inclusion of PAHs in a realistic chemical model at 
present means that one is obliged to treat them as small grains and thus follow 
^absolutely not comprehensive! 
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only their charge state. 

The present theoretical investigation on metastable C24H12 anions formation by 
electron resonant collisions should therefore be taken into account just in connec- 
tion with the above general comments: the coronene molecule, in fact, with an 
equilibrium structure of six aromatic rings surrounding a central exagonal C-ring 
and with a positive electron affinity (~ 0.47 eV, [|89ll ) whose value is comparable 
with that for two neutral tentatively observed R\HsQ, could be viewed as an im- 
portant and representative medium-size PAH among the surmised interstellar aro- 
matic molecules. Keeping in mind that the electron affinity is essentially indica- 
tive of the anion thermodynamical stability, looking for the possible presence of 
long-lived resonances in an energy range which agrees with the global environ- 
mental astrophysical conditions, might be helpful in characterize the feasibility 
of this PAH to efficiently form stable anions by electron attachment and, more in 
general, its feasibility to react with cold electrons. Furthermore, when considering 
the possible presence of C24H]^2 interstellar/circumstellar region as mainly 
produced by resonant electron collisions, some experiments [1771 [861 l96ll show that 
the undissociated anion production takes place for colliding particles having en- 
ergy very close to threshold while the vibrational excitation by electron impact 
should occur well below 10 eV. Since this energy range appears to be reasonable 
also for the general ISM conditions, (ranging from dense/dark clouds to diffuse 
clouds and hot circumstellar envelopes), we find useful to make a comparison 
between the experimental results and our findings in order to better characterize 
the behaviour of this molecule in an astrophysical environment where also free 
electrons are. 

In the interstellar space, in fact, free electrons represent an alternative way to 
transport energy; when compared to photons in the non-relativistic energy range, 
they have a much lower energy/momentum ratio, hence they can correspondingly 
carry much more momentum. In other words, although energy-exchange between 
the ionized gas, the neutral gas and the dust component might be mainly radiative 

^pyrene and anthracene, EA = 0.41 eV Ml 021 and 0.53 eV LlQSij, respectively 
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(although still depending on the environmental constraints), free electrons can also 
play a significant role because the neutral molecular gas can be heated by energy- 
transfer electron collisions. Moreover, in a partially-ionized plasma, the electron 
transport is primarily governed by the momentum-transfer cross-section for 
elastic scattering of the present electrons which is defined as 



da/dQ being the elastic differential cross section (DCS). Since experiments show 
the coronene molecule to be very "reactive" under low- and very low-energy elec- 
tron collisions [20,21], here we want to investigate the angular redistributions that 
the impinging particles undergo after the coUisional event with the title molecule. 
Therefore, in the framework of the multiple roles played by PAH molecules in in- 
terstellar or circumstellar environments, we find important to provide realistic esti- 
mates of the elastic DCS for electrons scattered off gas-phase coronene molecules. 
The angular redistributions of scattered electronic particles can in fact affect other 
physical processes: under the application of an external electric field (originated 
by radiowaves propagating through the plasma), electrons start to move away from 
the field direction where the degree of deviation of the electron trajectory is deter- 
mined, at least qualitatively, by the momentum-transfer which has occurred during 
the coUisional event. This means that the propagation (both reflection and attenu- 
ation) of radiowaves in a plasma can be considered as controlled by am [22] in the 
sense that, generally speaking, the larger the larger will be the deviation of the 
scattered electrons from the field direction, whereby the global electron motion 
would contiibute less to the conductivity of the partially ionized gases. 




(3.8) 
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3.2.2 The partial and total elastic cross sections at low energies 

We start our analysis of the electron scattering on gas-phase coronene by assuming 
that the C24H12 molecule is in its ground electronic state as well as in its equilib- 
rium geometry. The neutral closed shell target wave function was obtained using 
first an expansion with Gaussian functions centered on the nuclei of the molecule 
using a 6-31-i-i-G** basis set. The total SCF energy provided by calculations was 
—915.992787 hartrees. The equilibrium geometry for the coronene molecule ob- 



tained with the above Gaussian basis set is reported in figure 3.13 




Figure 3.13: Coronene otpimized equilibrium structure at 6-31++G** level 

According to the molecular symmetry point group to which the coronene molecule 
belongs in its equilibrium structure, the full nonspherical nature of the interaction 
potential with the incoming free electron was included in the present scattering 
calculations using the Deh symmetry point group. In order to include in the in- 
teraction potential the long-range polarizability contributions we used the model 



3.2. THE CORONENE MOLECULE, C24H 



133 



correlation-polarization potential Vcp with the total asymptotic tensorial polariz- 
ability of 105.928, 343.560 and 343.560 {Uxx, otyy and Uzz, respectively) calcu- 
lated by us. To provide the relevant scattering K-matrix elements, the parameter- 
free calculations at the adiabatic- static-model-exchange-correlation-polarization 
(ASMECP) level were carried out in the molecular frame of reference by solving 
the close-coupled equations where we have employed a partial wave expansion 
up to Imax = 60 both for the bound molecular orbitals and for the scattering wave 
functions, so that the corresponding expansion of the static contribution included 
2/max terms for both the electron-electron and electron-nuclei Coulombic interac- 
tions. For the S-matrix poles as well as the ICS and DCS calculations the size of 
the physical box that contains the diabatized partial wave terms was taken by us 
to be 50 A. 



200^ \ ' \ ' \ ' \ ' \ ' \ r 




eV 



Figure 3.14: Computed elastic (rotationally summed) integral cross sections. 



We show in figure 3.14 our computed total elastic integral (rotationally summed) 
cross sections in the 0-7 eV energy range. Keeping in mind that the reahstic. 
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Table 3.4: Computed one particle shape resonance real and imaginary energy compo- 
nents. 



albeit approximate, form of interaction with the free incoming electron, as well 
as the necessary truncation of the / expansion, could reasonably provide the main 
reasons according which the resonance energies could be located at somewhat 
different position than what expected by experiments, the perusal of the results 



shown in figure 3.14 already provide the following general information on the 
expected features of the low-energy dynamic response of the coronene molecule 
to its interaction with environmental electrons: 

i) in the 0-7eV energy window, the computed cross sections show rather marked 
structural peaks which can be identified as one particle shape resonances (i.e. 
formation of metastable C24H12 anions) in the Ei^ (E = 1 .277 and 5.86 eV), E2U (E 
= 2.40 eV), A2„ (E = 2.98 eV) Ai„ (E = 3.56 eV) and Bag (E = 4.80 eV) symmetry 
components of the D^h point group, respectively. The resonant parameters are 
summarized in Table 3.4, where we report the symmetry, the energy at which 
the metastable anion is formed, the imaginary component width (both in eV) and 
the corresponding lifetime; within the last column on the right we report the way 
we have assessed them ('BW refers to the occurrence of a clear tt jump in the 
eigenphase sum, hence the ICS lorenzian peak is fitted at the Breit-Wigner level; 
'S' means identification via the associated S-matrix pole on the unphysical sheet 
of the complex energy plane). 

Furthermore, we find it important to point out now that the likely presence of vi- 
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brational coupling effects, which are not considered by the present fixed nuclei 
calculations, could cause the above electron attachment cross sections to be both 
broader and smoother. It is also interesting to note that in the energy range con- 
sidered in the present investigation, which contains the main contributions from 
electrons having a mean kinetic temperature Tg < 10.000 K, the computed ICS 
are not particularly large: apart from the sharp resonance features (thus neglect- 
ing for the moment the possible resonant capture of the incoming charged particle) 
and the rapidly increasing branch close to eV, the background cross sections are 
in fact a little greater than lOOA^, a value which is not much different from the 
rigid sphere limit given by the geometrical size of the neutral system with a di- 
ameter of roughly 10 A. This feature is evident between 0.5 eV and 2 eV, a range 
corresponding to electronic kinetic temperatures of Tg ~ 5000 K and ~ 20000 K, 
respectively; 

ii) above 2 eV we find five relatively crowded and rather marked maxima cor- 
responding to five resonances which, with the exception of the Aiu symmetry 
contribution, are relatively broad. Their resonant character was confirmed by both 
the TT jump in each eigenphase sum and by locating the corresponding S-matrix 
poles on the unphysical sheet; 

iii) below the collision energy of 2 eV, we locate instead only one very narrow peak 
of Eig symmetry at 1 .277 eV. The resonant character of this feature is provided by 
the clear and very sharp tt jump in the corresponding eigenphase sum, which we 
have fitted to a Breit-Wigner formula, yielding a resonance width of 0.0046 eV; 

iv) although the molecule has not a permanent electric dipole moment, due to the 
highly symmetric structure of its equilibrium geometry, our calculations show a 
markedly increasing behaviour of the elastic ICS close to threshold. One can sur- 
mise here that this feature indicates that the slowest (the 'coldest') component of 
the free-electron gas present in the ISM can strongly interact with the gaseous neu- 
tral C24H12 molecule, if present, as already suggested by laboratory experiments 
[iTTl [861 and as we shall further discuss below. 

In order to gain a deeper insigth in the different contributions which make the total 
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ICS, we report all the partial ICS coming from the twelve irrreducible representa- 



tions of the D^h point group in the twelve panels of figure 3.15 Within our present 
modeling of the scattering processes it is obviously important to ensure that all the 
necessary indicators of numerical convergence (e.g. size of radial integration box, 
size of the multipolar expansion of the interaction potential, size of the partial 
wave expansion for the scattered electron, etc.) be tested to provide final, con- 



verged cross sections. To this purpose, we thus report in each panel of figure 3.15 
the partial ICSs computed using different values for the interaction box size (50 
A(blue straight line) and 80 A(maroon diamonds)), as well as two different values 
of Imax (Imax = 60 (bluc Straight line) and Imax = 80 (maroon diamonds)) for the 
corresponding partial wave expansion. As showed in that picture, the numerical 
convergence is already reached when using a box of 50 Aand assuming l^ax = 60 
in the sense that, in the present case, the final cross sections can be reasonably 
considered as converged within a few percent of their reported values. 

From a purely physical point of view, we can now add a few additional com- 
ments on the general and global behaviour of the ICSs for the coronene molecule 
under low-energy electron collisions. First of all, in connection with the Ai^ par- 



tial contribution (first upper panel on the left in figure 3.15), it is important to 
stress the ICS behaviour close to zero energy: such a threshold behaviour, which 
is distinctive of the totally symmetric component (the one including the s-wave 
scattering contribution), together with the large and negative scattering length (a 
= -17.63) obtained at the smallest energy value examined by the present calcula- 
tions (E = 10~''eV), suggests the formation of a zero energy resonance ll95l due 
to the 'strong' long-range attractive potential provided by the large polarizability 
of this molecule. However, as also shown by the other panels of the same fig- 
ure, we further see that the Ei^ and symmetry components also contribute to 
the strong rise of the total ICS in the low-energy range, since each of them has 
a well resolved maximum when the collision energy downs to 0.06 and 0.21 eV, 
respectively. We can then surmise the likely formation of a loosely quasi-bound 
state [|95l with components in all these symmetries. Such a feature will however 
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discussed more in detail in the next sections. To conclude the present section, we 
underline that our calculations indicate markedly enhanced interaction effects be- 
tween electrons and coronene close to the vanishing energy range of the surround- 
ing electrons. The implications for its possible behaviour, if present, in interstellar 
dense clouds (low Tg) and protoplanetary atmospheres (high Tg) are therefore that 
most of the efficient probing of this molecule by the surrounding electron swarms 
will be occurring via the very low-energy tail of their distributions and chiefly 
through the lower partial waves describing the interacting electrons. Before to 
analyse the angular redistributions of electrons scattered off coronene molecule, 
in the next section we will focus on the mechanisms responsible for the electron 
capture process, discussing our findings first of all for very low-energy impinging 
electrons (0-2 eV) and then for little more energetic (2-7 eV) colliding particles. 
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Figure 3.15: Computed partial ICS for the tewlve irreducible representations of Dg/j 
point group. 
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3.2.3 Threshold behaviour and resonances in electron scatter- 
ing from gas-phase coronene 

In the present section we focus more in detail on the mechanisms which could be 
active and then responsible for the dynamical trapping of the electronic impinging 
particle in each of the previuosly located shape resonances and in particular when 
the colliding energy is moving down to its threshold. We shall also discuss in 
what follows some non-resonant features associated with the maxima in the Aig, 
Ei„ and symmetry contributions. We shall in fact speculate that such features, 
even if not properly ascribed to resonances, could indeed play a role in trapping 
the incident electron, and therefore for the production of the stable bound anion: 
the lack of a permanent electric dipole moment for this molecule could qualita- 
tively justify their non-resonant character in the framework of our fixed nuclei 
approximation. 

Generally speaking, at impact energies below 10 eV, the interaction of an elec- 
tron with a molecule is mainly dominated by the existence of short lived anions 
(called 'resonances') due to the temporary capture of the incident electron by the 
target molecule. From the experimental point of view, when the lifetimes are 
'short' (e.g. 10^^^ sec ~ 10^^^ sec) their direct observation is generally rather 
challenging: they are thus observed in many cases through their decay, either by 
autodetachment or by dissociation into different fragments, one of which carries 
the extra electron. If the lifetime of the resonance under investigation lengthens, 
the ejection of the extra electron can also occur at an internuclear distance differ- 
ent from those of the ground state, hence giving rise to vibrational excitation. This 
channel is however in competition with the channel which yields a stable negative 
fragment plus one or several neutral ones via the dissociative electron attachment 
process. 

The most recent among the available experiments shows that the major resonant 
interaction takes place for low and especially vanishing energies iTTTll . where for 
higher collision energy values (1.5, 2.0, 2.5, 3.5, 4.5 7.0 and 9.0 eV) five groups 
of excited vibrational modes are observed with electronic energy losses ranging 
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from 0.065 eV to 0.375 eV confirming that several more energetic resonant states 
appear to contribute to the global interaction. These experiments also confirm 
that the formation of the corresponding not dehydrogenated negative ion occurs 
at zero energy, and at the same time they shed new light on the real behaviour of 
the C24H12 molecule under very low-energy electron collisions: if on one side the 
presence of one intense peak centered at zero energy in fact confirms the findings 
of [l86ll . the better electron energy resolution allows to demonstrate the existence of 
another peak at 0.3 eV, in sharp contrast with the previous data. The nature of this 
second, less intense, peak is not obvious: the authors argue that it can therefore 
correspond either to a particular long-lived shape resonance or to resonant state 
having a core excited character, the latter reasonably not directly coupled to the 
ground state of the anion. Moreover, since the spacing between this peak and the 
larger one at zero energy is reminiscent of a specific C-H stretch mode (0.375 
eV), a vibrational Feshbach resonance associated with this mode could also be a 
possibihty. 

Keeping in mind that elastic scattering and electron capture cross sections can of- 
ten be related since the same resonant features can in principle contribute to both 
processes, we intend now to analyse and discuss our findings by making a com- 
parison (when possible) with experiments in order to provide some useful infor- 
mation about the possible formation channels for low-energy metastable negative 
ion (C24Hj^2)* which can actually be the gateway for the production of the ther- 
modynamically stable, long-lived parent anion, the electron affinity for the neutral 
species being in fact positive: -1-0.54 ±0.1 eV according to the kinetic measure- 
ments in [l84l and -1-0.47 ± 0.09 eV according to the photodetachment experiments 
reported in [|89l . Since the theoretical modelling of low-energy electron scatter- 
ing for large polyatomic molecules is indeed computationally very challenging, 
we have employed in the present study our ASMECP model of potential interac- 
tion which has been previously shown to provide good results not only for small 
molecular objects: it has been in fact proved that such a potential is accurate 
enough for examining the low-energy attachment on the systems as large as the 
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Ceo molecule ( IfTOTTl ). 



Non resonant very-low energy features: Ei„ and partial contributions. 



The two panels of figure 3.16 report the behaviour of our computed ICS in the 
0-2 eV energy range for the partial contributions where we find a pronounced 
maximum, clearly below 0.5 eV. 




Figure 3.16: Low-energy (0-2 eV) ICS behaviour for the Ei„ and E2g symmetry. See 
text for details. 

The blue squares in each panel refer to ICS calculations including all four interac- 
tion terms (Coulombic, exchange, correlation and polarization) as previously il- 
lustrated (see chapter 2), while the triangles refer to the computed ICS in the same 
symmetry contribution but without including the correlation-polarization poten- 
tial. Yep- We can easily recognize that for the Ei„ and E2g symmetries the low- 
energy maxima, located at 0.06 eV and 0.22 eV respectively, are not numerical 
artifacts: due to the large polarizability of the neutral target molecule, originating 
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from the large number of delocalized tt electrons, such maxima could be in fact 
indicative of the true physical response of the target to the incident charged pro- 
jectile at a very low energy. In both cases, in fact, qualitatively we see that the 
colliding electron has not enough kinetic energy to come very close to the target 
molecule, causing the short-range interaction terms (correlation and exchange) to 
be less important then the long-range interaction terms. Thus, in order to provide 
an accurate physical description of the electron-C24Hi2 molecule collision, the 
asymptotic part of the long-range static interaction can be reasonably expressed 
as a sum of multipolar moments where the dipole term (the first one) is zero in the 
present case due to the inversion center that characterizes the present molecule at 
the equilibrium; furthermore, besides the 'pure' static interaction, the molecule is 
also perturbed by the presence of the incident electron which in fact induces an ad- 
ditional smaller dipole term, the effects of which is then considered by including 
our Yep potential in the ICS calculations. In other words, since such a behaviour 
in the ICS for the Ei„ and symmetries seems to be dramatically depleted when 



the long-range effects are not considered in the calculation (see figure 3.16), such 
maxima tell us that the impinging electron induces an electric dipole moment, de- 
pendent on the molecular polarizability, which is strong enough to cause both the 
above ICS to have a pronounced peak close to threshold. 

However, their not properly resonant character is confirmed by analysing the cor- 



responding eigenphase sums that we depict in figure 3.17 as a function of collision 
energy. 

As the panels corresponding to the Ei„ and E2g symmetries show, we do not have 
a Tc jump around the energies corresponding to the maxima: for both of them, in 
fact, the rapidly rising cross sections moving away from threshold correspond in 
fact to a slow eigenphase rise. The above low-energy features in both the partial 
ICS here under investigation are then followed by a Ramsauer-Townsend-type 
structures which occur in the energy window where the associated eigenphase 
sum crosses zero (for the E2g symmetry, to be precise, there is a little shift in 
energy of about 0.9 eV). In order to glean further physical information about this 
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Figure 3.17: Low-energy (0-2 eV) eigenphase sum for the Ei^ and £2^ symmetries. 

couple of maxima, we have computed the corresponding diabatic potential curves 
that we report in the insets of figure [3TT8 



It is thus clear from them that the energy of each maximum corresponds to the top 
of the corresponding barrier for the lowest angular component in the computed 
diabatic potential curves: / = 1 for the 0.06 eV peak in the Ei„ symmetry and / 
= 2 for the 0.22 eV peak in the E2C, symmetry, respectively. Additionally, in both 
symmetries it is possible to recognize a rising eigenphase sum close to threshold 
and which has a maximum immediately after the top of each of the above potential 
barriers is reached, as shown by the vertical black dashed lines within the insets in 



figure 3.17 We can therefore conclude that the analysis of these two features leads 
us to say that they do not appear to be the result of some actual resonant structure 
since none of the required indicators does exist. Qualitatively speaking, and using 
a semiclassical picture, we find it possible to speculate that these low energy rising 
eigenphase sums might be the quantum analog of the classical orbiting, where 
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Figure 3.18: Computed lowest angular momenta diabatic potential curves for the Ei 
and symmetries. 



there is an increase of time delay that in fact takes place when the energy of the 
colliding electron equals the energy value associated to the maximum of its own 
centrifugal barrier. This effect is reminiscent of that in a previous investigation 
[|90ll that focused on the benzene molecule in its equilibrium geometry (which 
belongs to the same point group of the molecule here under investigation); in 
the present case, however, such non-resonant maxima are more pronounced and 
markedly stronger in intensity if compared with those that characterize the same 
partial cross section contributions for the CeHg molecule. Thus, according to 
the experiments [TTTll . we surmise that even if such maxima are not effectively 
describing an actual dynamical trapping of the electron, the markedly enhanced 
static polarizability for the coronene molecule could indeed make them participate 
in the formation of the parent anion at very low energy via such temporary states. 
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Aig contribution threshold behaviour: zero energy resonance. 

The totalsymmetric partial contribution also deserves a detailed analysis: as we 
shall discuss in what follows, it plays in fact a key role in the formation of the 
resonant metastable anion at vanishing energies. 



We show in figure 3.19 our computed Ai^ cross sections (left panel) and the cor- 
responding eigenphase sum in the energy range between 10""^ eV and 2 eV. 




Figure 3.19: Left panel: behaviour of the Ai^ elastic (rotationally summed) ICS. The 
blue squares are referred to an ASMECP calculation; the red diamonds show the ICS 
behaviour when omitting the Yep interaction potential from the calculation (ASME calcu- 
lation). Rigth panel: Computed ASMECP eigenphase sums for the same symmetry and 
for the same energy window. 

The presence of a marked and very rapidly rising threshold value (for 10^^ eV our 
computed ICS value is 1091 A^) in the totally symmetric component, followed by 
a Ramsauer-Townsend minimum located at 0.06 eV as the corresponding eigen- 
phase sum moves down crossing the zero value, are features which demonstrate 
that the dominant contributing channel is the s-wave channel. In particular, as for 
the case of electron-Ceo molecule low-energy scattering HlOlll . we note a rising 
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eigenphase sum at threshold which is clearly shown between 0.0001 and 0.01 eV 



in the right panel of figure 3.19 such feature is indeed indicative of a virtual state 
[EQl. 



To better characterize the likely presence of a virtual state, we show in figure 3.20 



the behaviour of the Ai^ adiabatic scattering wavefunction for very low-energies 
as a function of the radial coordinate r. 




Figure 3.20: Behaviour of the Ai^ adiabatic scattering wavefunction for very low- 
energies as a function of the radial coordinate r (0 < r < 15 A). 

As that figure clearly shows, the adiabatic Ai^, resonant wavefunction becomes 
increasingly flattened, smoother and less prominent as the collision energy ap- 
proaches zero. The corresponding scattering length becomes negative and very 
large (a = -17.63) when computed at the lowest energy value (0.0001 eV) here 
considered. When using a computationally more demanding partial wave expan- 
sion, i.e. varying l^ax from 60 up to l^ax = 80, we obtain a = -19.53. According 
to the geometrical interpretation of a, this strongly suggests the existence of a 
zero energy resonance (a — > oo for an s-wave virtual state) as the intermediate 
physical situation between an attractive potential which for vanishing energies is 
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actually able to support an s-wave bound state (a > 0), and an attractive potential 
which instead is not deep enough to produce an s-wave bound state (a < 0) in the 
energy range close to zero. The definitive confirmation that we are thus dealing 
with a zero energy resonance comes from the S matrix pole for the Ai^ scattering 
symmetry, that we locate at an energy of -0.037 eV on the unphysical Riemann 
sheet (associated to a pure imaginary negative linear moment of -0.052 15z a.u.). 

These findings provide an intriguing and interesting explanation for the experi- 
mental scattering data reported by [r77l[86l . raising the possibility that the bound 
anion observed close to zero energy may have its origin in a virtual state scat- 
tering feature. This argumentation is in fact only seemingly in contrast with the 
experimental findings, since in scattering of this nature the zero energy resonant 
negative ion does not form and indeed cannot be accessed, unless some changes 
in the nuclear geometry do occur. In other words, we can surmise that the phys- 
ical process leading to the zero-energy resonant anion is an s-wave dominated 
collision where the interaction, in the framework of the fixed nuclei approach, is 
not able to support a bound state. However, the large polarizability of the target 
molecule allows us to suggest that there is time enough during that special com- 
plex formation for nuclear motion to take place in order to access the bound state, 
the latter being in fact located in a classically inaccessible regime. In this sense, 
the virtual state scattering can thus act as a doorway state for the formation of 
negative ions like it does for Ceo electron scattering [I1091I101I ). Given the large 
number of closely spaced vibrational modes for the C24H12 molecule, we find it 
very likely that such a virtual state may decay via coupling with the dense mani- 
fold of excited roto-vibrational states in a sufficiently short time to compete with 
the autodetachment process, as indeed shown by the experiments [1771 . The nega- 
tive ion resulting by virtual state scattering is in fact a metastable species, even if 
for a small amount of energy, and thus its lifetime should be long enough so that 
the excess energy could be efficiently dissipated either by photon emission or via 
intramolecular vibrational redistribution (IVR). 

Hence, the coronene molecule represents an example of medium size R^.H species 
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where virtual state scattering plays an important role for its parent negative ion 
formation. 

To conclude the present theoretical analysis of the threshold electron scattering 
behaviour for the coronene molecule at its equilibrium geometry, we report in 



figure 3.21 the computed squared modulus of the lowest angular contributions (/ 
= and Z = 2) of the totalsymmetric adiabatic scattering wavefunction when the 
incident energy is 0.5 eV (upper panel), as well as the diabatic potential curves 
that correspond to the above partial angular contributions at the same energy. 
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Figure 3.21: Upper panel: squared modulus of the lowest angular momentum contribu- 
tions (dotted-dashed line: / = 0; red straight line: / = 2) of the totalsymmetric adiabatic 
scattering wavefunction when the incident energy is 0.5 eV vs the radial coordinate r. 
Lower panel: diabatic potential curves that correspond to the above partial angular con- 
tributions computed at the same collision energy and depicted in the same radial domain. 

When the incident projectile has an energy of 0.5 eV, we can tentatively note that 
the / = 2 angular component ([^, included in the A\g scattering wavefunction, 
is further able to yield the same feature previously discussed for both the 

'^a d22-like angular wavefunction which belongs to the Aig symmetry; see the character table 
of the Dg/i point group of symmetry. 
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and symmetries, even if in the present case for an higher energy and with a 
comparatively reduced intensity. The Ai^, ICS shows in fact a peak at ~ 0.5 eV 
which is much broader when compared to those of Ei„ and symmetry (see the 



top right panel in figure 3.15). Furthermore, such a collision energy corresponds 



to the top of the / = 2 centrifugal barrier (lower panel of figure 3.21 1, and the 
squared modulus of the totalsymmetric scattering wavefunction, computed at 0.5 



eV (see upper panel of figure |3.21[ red solid line), intriguingly exhibits just a 
small probability density mainly located at the outer edge of the region occupied 
by the target C24H12 molecule. However, since in the present case this effect is 
not caused by the lowest angular momentum contribution (/ = 0), we find such 
a probability density to be negligible so that the / = 2 partial wave contribution 
within the Ai^ symmetry should not be very effective for the low-energy trapping 
of the incoming electron. 



The low-energy shape resonances. 

Having discussed our findings for the marked but non-resonant features which 
characterize the Ei^, and Aig scattering symmetries, we now pass to anal- 
yse more in detail the six one-particle shape resonances introduced in a previous 
section, starting from that one at the lowest energy. To this purpose, we report 



in figure 3.22 the shape of the wavefunction corresponding to the very narrow 



Eig shape resonance that we locate at an incident energy of E = 1.277 eV. Within 



the two panels of figure 3.23 we further show the associated diabatic potential 
curves for the lowest angular contributions (/ = 2, 4, 6; lower panel) involved in 
the dynamical trapping of the electron and (upper panel) the radial behaviour of 
the resonant wavefunction for each of the above angular contributions. 



It is then clear from the two panels of figure 3.23 that in such a very narrow reso- 
nance the dominant dynamical trapping of the impinging electron arises from the 
/ = 4 centrifugal barrier, where the tunnelling would have to span a rather large 
radial region (4 < r < 7.5 A). Excluding for the moment the possible decay into 
a bound stable state, in most elastic one-electron resonances originated by an an- 
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Figure 3.22: Spatial sliape of the low-energy Eig (E = 1.277 eV) resonant wavef unction. 

gular momentum barrier tliere are two possible mechanisms for their decay: the 
resonant state can either give rise to a direct decay since its probability density 
can tunnel backward through its own dominant centrifugal barrier, or it can decay 
due to a dynamical coupling with other partial waves having lower / values. In 



the present case, as shown pictorially in the lower panel of figure 3.23 the prob- 
ability densities associated with the lowest angular contributions {1 = 2 and / = 4) 
can provide such a coupling, although a possible nonadiabatic coupHng between 
the / = 4 and / = 6 partial resonant waves, occurring just behind the / = 4 barrier, 
can play a role in preventing the rapid ejection of the extra electron, thus provid- 
ing an additional qualitative explanation for the very long lifetime that we find 
to be associated to this resonance. It is furthermore interesting to note that this 
resonance belongs to the same symmetry of the C24H12 HOMO-1 as well as in 
the LUMO (j^; in a qualitative sense, the above Ei^ resonance should then be ex- 
pected to be overlapped (i.e. interacting) with the lower-lying bound states which 
are more likely to participate in the formation of a stable bound anion. In other 
words, according to our equilibrium geometry analysis, we surmise that such a 
very long-lived metastable anion could play a role in the formation of the neg- 
ative parent species by dissipating the fairly small excess energy content via a 



'where, making a comparison between our resonant 3D map and the 3D map for the Eig 
LUMO calculated by 11771 at DFT level but with the same basis set, we can tentatively note some 
resemblance 
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Figure 3.23: Upper panel: computed diabatic potential curves for the low-energy Eig 
(E = 1.277 eV) resonance. Lower panel: squared modulus for the adiabatic resonant 
wavefunction lowest angular contributions. Black-dashed line: / = 2; red circles: / = 4; 
green circles: / = 6. 



radiationless intramolecular vibrational redistribution process. 



In figure 3.24 we collect the real part of the wave functions for the remaining res- 
onances, all of them occurring for energies which are greater than 2 eV; the above 
figure, from top to bottom, contains the panels ordered from the lower resonance 
energy (of the E2„ resonance, E^ = 2.26 eV) moving up to the higher one (Eig, E^ 
= 5.79 eV). 

On the basis of these findings, to conclude this section, we tentatively try to make 
a comparison with the available fairly detailed experimental data Il86ir77]| . We can 
thus summarize the main points of contacts in the following way: 

i) in the considered energy range, which is the one that is here astrophysically 
relevant, our calculations do not show the presence of any a* resonance, as directly 
proved by the real 3D maps of the resonant wavefunctions that we collect and 



depict in figures 3.22 and 3.24 
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ii) the resonance ^Ai„ located at E^ = 3.62 eV, due to its relatively long lifetime 
(r2^2„ ~ 253 -10^^^ sec, see table 3.4), could be however involved in the exci- 
tation of some CH stretching modes: we can tentatively surmise it by looking at 
the spatial shape of the corresponding metastable wavefunction, depicted in the 



central inset of figure 3.24 since only in such a resonant state, in fact, the ex- 
cess electronic wavefunction seems to efficienly 'cover' the periferal CH groups 
although located above and below the molecular plane. Furthermore, we find it 
interesting that this resonant feature occurs at an energy which falls almost at the 
onset of the energy range where the experiments [TTTll show the uch = 0.375 eV 
mode to be quantitatively excited; 

iii) in general terms, the presence of several tt* resonances up to about 7 eV for the 
coronene molecule is not surprising: besides the fact that the associated electron 
affinity is positive, [IMllHSl, the experiments indeed reveal that at incident energies 
higher or equal than 2.0 eV the variation of the relative intensities of different vi- 
brational bands indicates the occurrence of several resonant excitation processes. 
It would also be helpful [fTTl to investigate the coronene molecule by means of the 
electron transmission spectroscopy, since the findings of such experiments might 
enable us to have a better comparison with the present results, thus clarifying the 
role played by its metastable anions in the more general astrophysical context in 
which such electron-driven reactions have been studied. 




Figure 3.24: Spatial shape of tlie resonances located over 2 eV; see table 3.4. From top 
to bottom: ^E2u, ^A2„, ^Ai„, ^B2g and ^Eig 



154 



CHAPTER 3. RESULTS AND DISCUSSION 



3.2.4 The scattered electron angular distributions at low ener- 
gies 

The previous subsections have described in detail the behaviour of the elastic, 
partial and total, ICS for an energy range of the impinging electrons which is as- 
trophysically relevant. Now, the above results further suggest to be of interest 
knowing how the angular distributions would behave within the same range of 
energies, even if, due to the lack of specific laboratory data, we cannot make a 
direct comparison of our results with measured angular distributions for gas phase 
coronene molecule. However, the present computational experiment already al- 
lows us to draw the following features from the computed behaviour that we report 



in figures 3.25 and 3.26 
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Figure 3.25: Computed differential elastic (rotationally summed) cross sections in the 
'low-energy' domain. See main text for details and discussion. 

i) at the lowest collision energies considered here {E = 10~^ eV), the scattering 
is dominated by the dipole-induced interaction which, according to the results 
previously discussed, is however not able to support a bound state; as shown by 
figure 3.25[ the corresponding DCS is found to be very little angle-dependent and 
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Figure 3.26: Computed differential elastic (rotationally summed) cross sections in the 
'high-energy range'. See main text for details and discussion. 



largely isotropic, thus confirming once again a behaviour chiefly due to "s-wave" 
scattering at those low energies; 

ii) in the energy range between 10^^ eV and up to 0.5 eV (see figure 3.25), with 
dominant contributions from Ei„ and symmetries to the corresponding ICS 



(see figure 3.15), our findings show a progressive increase of the forward peak, 
as well as a marked reduction of the backward scattering as indicated by all the 



panels of figure 3.25 



iii) as the collision energy moves across the narrow Ei^, ICS maximum associated 
to a strong resonance of that symmetry, the computed angular distribution at E = 



1 .277 eV given by the top row of panels in figure 3.26 becomes much flatter (from 
9 = 60° up to = 120°). Such a feature occurs as soon as the energy leaves the 
very narrow resonance region but resumes again up to at least 1 .5 eV; 



iv) the data from figure 3.26 show that when the colhsion energy increases to 
higher values (which are still astrophysically relevant), we see that the DCS, as 
a function of the scattering angle 9, shows an increase of the forward peaks at 
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smaller angles while the backward scattering becomes largely negligible. 
As a general trend in the energy range of interest, we thus find as dominant the 
presence of the forward scattering behaviour, a feature that involves a rather large 
geometrical cone (Omax ~ 60°). Furthermore, the forward scattering appears to 
be dominant also when the collision energy matches the Eig resonance position 
(E = 1.277 eV), where using a qualitatively semiclassical picture, the 'orbiting' 
of a trapped electron in a shape resonance should distort the angular distribution 
generated by the induced dipole on the target molecule. 



3.2.5 Computed momentum-transfer cross sections and the rate 
coefficients 



In figure 3.27 we report again, as a remainder, the ground-state total elastic ICS 
(solid line) together with our computed momentum-transfer cross sections (MTCS; 
black right-oriented triangles), both as a function of the collision energy. The red 
circles refer to the elastic ICS computed by numerical integration of the differen- 
tial elastic cross sections daldVL: 

PIT 

a^nt = 271 / —Ck-f)sined9 (3.9) 

^ 

while the blue solid line reports the ICS that were calculated directly from the K 
matrices, i.e. without the recoupling operations necessary to express the DCS in 
terms of K matrices [I104II . The overall trend for the MTCS is seen to follow the 
ICS profile but producing consistently smaller values, as expected from the be- 
haviour of the DCS discussed before. The rate coefficients can in turn be obtained 
by integration over the energy distribution of the scattered electrons present in the 
astronomical environment 



8 K T 1 

k{Te) = \l TfTT^ I (j{E)EeT^dE (3.10) 



vr/x {KTfJ, 

as already discussed in our introduction. They are computed using both the present 
ICS {Oint) and MTCS {am) values, and are summarized in table 3.5 as a function 
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Figure 3.27: Comparison between computed elastic integral cross sections and momen- 
tum transfer cross sections. The blue straight line reports the ICS directly computed via 
K matrix elements; the red circles and the black triangles refer to the ICS and MTCS 
computed by numerical integration of the differential elastic (rotationally summed) cross 
sections. See text for details. 

of the kinetic electronic temperature, Tg. In that table the first column of data 



refers to a numerical quadrature using the integral cross sections of eq 3.9 while 
the second column shows the data obtained from the momentum-transfer cross 
sections, that for elastic scattering are defined as 

am = 271 (1 - cosO) — sine d9, (3.11) 
daldVt being the elastic differential cross section (DCS). 

An interesting comparison is provided by the Langevin rate coefficient, which is 
given by the following expression 



27rey/^, (3.12) 
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Te 


{cm? ■ sec~^) 


A;^(Te) = (a„(^)v^) 
{cm^ ■ sec~^) 


10000 


5.34-10-9 


3.24-10-9 


5000 


3.94-10-9 


2.67 -10-9 


2000 


3.22 -10-9 


2.53 -10-9 


1000 


3.02 -10-9 


2.49 -10-9 


500 


2.88 -10-9 


2.25 -10-9 


250 


2.68 -10-9 


1.98 -10-9 


150 


2.43 -10-9 


1.83 -10-9 


100 


2.16 -10-9 


1.75 -10-9 


50 


1.76 -10-9 


1.62 -10-9 


30 


1.56 -10-9 


1.55 -10-9 



Table 3.5: Computed electron-molecule rate coefficient using the cross sections from 



eq 3.9 (left column) and those from eq 3. 1 1 (right column) 



a and yU being the mean polarizability of the molecule entering the collision and 
the reduced coUisional mass respectively. After substituting a = 264.35 bohr^, 
such an expression returns the value = 4.57- 10-9 cm^ sec~^, which is not 
dependent on the temperature T^.. 

We find that some comments should be now made on the above results: 
(i) first of all, our computed rate coefficients (both kint{Te) and kmiT^)) follow 
the same trend in the sense that they both decrease as the Tg temperature moves 
down to the lowest value considered (Te = 30 K); (ii) they reach almost the same 
value (~ 1 .55/1 .56- 10^'^ cm? sec^^) when Te = 30 K, according to the fact that both 
ICS and MTCS decrease close to each other in the very-low energy range, (iii) 
The increasing of Te up to 10000 K causes both kint{Te) and km{Te) to increase, 
although /cmCat Te = 10000 K) is lower than the Langevin ki rate coefficient value: 
the latter thus appears as an upper limit to the present modelling. Hence, in this 
framework, the Langevin rate coefficient turns out not to be 'accurate' enough, 
expecially at low electronic temperatures (Te < 150 K). 



3.2. THE CORONENE MOLECULE, C24H 



159 



Since our DCS as well as our ICS are both rotationally summed, we cannot sep- 
arate the rotationally elastic and inelastic contributions to the final rates. Fur- 
thermore, within the FN approach we are neglecting the vibrational degrees of 
freedom; consequently, given the present collision energy range and the initial in- 
ternal temperatures of the 'large' molecule of interest in the present study, (e.g. 
see also ref. [|9TI ). little vibrational excitation effects should be present as in fact 
confirmed by the results of llTTll . In other words, since the vibrationally inelastic 
processes require stronger penetration of the molecular volume by the impinging 
electron, they should correspond to different angular distributions for the scattered 
electrons: however, we expect that for the present case the latter would provide 
only small variat ions. We are in fact here interested in the global behaviour 
of the ICS and MTCS in providing a reasonable estimate for the TNI formation 
rate constants which could be thus viewed as restricted with respect to our esti- 
mates. Our computed features for the scattered electrons would therefore lead us 
to suggest that gaseous coronene molecules are indeed relatively 'transparent' to 
free electrons in the considered energy range, a conclusion also confirmed by our 
computed DCS, which indicate negligible electron scattering flux for values of 6 
larger than 60°. 

3.2.6 Present conclusions 

In this work we have endeavoured in order to analyse in some detail the behaviour 
of a very typical polycyclic aromatic system which is considered to be part of the 
possible PAH species existing in different regions of the ISM like diffuse clouds 
[l65l . dense molecular clouds [I108II and in circumstellar atmospheres UlUL 
In particular, having discussed in a previous section ([^ that another possible can- 
didate as a PAH precursor, the ortho-benzyne molecule, does indeed exhibit strong 
interaction with the electrons present in the astronomical environments under con- 
sideration, we have deemed to be now important to move to a much larger partner, 
the gas-phase coronene molecule, which has been often considered important in 

'"and discovered in our earlier published work 11211 [STl 
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early R\H studies [78], in order to try to understand its possible behaviour, at 
the molecular level, as a deflector of low-energy electrons present in those astro- 
physical contexts and therefore understand better its possible role as an energy 
absorber from the environmental electrons or within a low-temperature molecular 
plasma produced in those astronomical surroundings. We note moreover that one 
important point that makes this comparison intriguing is provided by the electron 
affinity for such a couple of molecules; despite in fact their marked difference 
in size (thus in the number of normal modes) as well as the presence of a rather 
large permanent dipole that characterizes the o-benzyne (1.68 D, [i97l '). the EA 
value is similar: about -1-0.5 eV for C24H12 (-1-0.54 ± 0.1 eV according to the 
kinetic measurements as provided by [l84ll and -1-0.47 ± 0.09 eV as measured by 
photodetachment experiments reported in [[891 ) and -1-0.56 eV for 0-C6H4 [[21 I100L 
respectively. This suggests that thermodynamically both of them should be char- 
acterized by almost the same capability to form the low-lying parent not dehydro- 
genated metastable negative ion; according to this, we have consequently studied 
their kinetics in forming the metastable anions, as provided by our estimated cor- 
responding rate coefficients. 

The calculations of the present study are indeed the first example of such com- 
putational investigations on a Targe' ([^ member of the R^Hs family, one which 
reasonably can be considered to be present in several different astronomical envi- 
ronments ( [|108[[80l ) where, together with larger R^.Hs, plays a role in contributing 
to both their chemical and physical properties ( [|108l[T02l[82l[80ll ). Tackling such 
a large species with reliable ab initio methods is furthermore meant to provide a 
link with earlier studies on smaller members of that large family of molecules, 
and therefore to clarify to us the relative importance of at least two emblematic 
systems: benzyne and coronene; at the same time, one can also speculate that 
such an endeavour could be useful to better extrapolate which is the behaviour 
of polycyclic aromatic molecules interacting with low-energy electrons, moving 
from small-size PAHs (benzene, benzyne) toward large-size PAHs (with more 

"where from the computational point of view, large means demanding 
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than 30-40 C atoms): the coronene, in fact, as a medium-size PAH without per- 
manent dipole, can be viewed as representing one possible ideal link between the 
previous two PAHs classes. 

The present findings come from a realistic, quantum treatment of low-energy elec- 
tron scattering and provide both integral and momentum-transfer cross-sections 
as well as angular distributions of the elastically scattered electrons over a broad 
range of collision energies, i.e. covering the range of the likely energetic of elec- 
trons in dense molecular cores and in protoplanetary atmospheres. 

The detailed search of the possible formation of threshold, metastable anions for 
the coronene molecule was carried out and discussed under the ligth of the most 
recent available experimental data: our findings point toward virtual state scatter- 
ing (zero energy resonance), occurring for vanishing energies, as the likely door- 
way to the production of the parent metastable not dehydrogenated negative ion. 
This allows us to surmise that the nonadiabatic couplings of the trapped electron 
with the nuclear motion during the near zero-energy electron-molecule encounter 
could be particularly efficient to access the stable bound state. Conversely to 
what is thought to occur for the benzene and naphtalene molecules, due to their 
negative EA, for the present molecule we are thus led to confirm an active in- 
volvement of the virtual state scattering in forming the corresponding resonant 
anion. Thus, more in general, the coronene molecule would constitute a reliable 
first example where the virtual state scattering indeed plays a crucial role in the 
complex interstellar processes responsible for the production of stable negative 
ions, which currently are in turn understood to be foundamental both in modify- 
ing, even substantially, specific molecular abundances HlOSi l99l [8OII and even the 
global physical conditions ll82l II 0211 in several astrophysical environment. In this 
framework, the coronene molecule can therefore be viewed as a relatively large 
polycyclic aromatic species that, on the basis of the results obtained in this study 
about its capability to form resonances in the 0-7 energy range at the equilibrium 
geometry, can form by different paths bound stable anions at low- and very low- 
energy in competition with autodetachement channels. From the astrophysical 
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point of view, this means that the C24H12 molecule, when considered either in a 
diffuse or in a dense and cold molecular cloud, should be able of soaking up a 
rather large amount of free electrons but only at selected energy windows and not 
over the broad spectrum of energy distributions of the environmental electrons. 
This therefore confirms our previous conclusions about the small rate coefficient 
for the metastable anion formation and on its 'relatively transparency' under low- 
and very low-energy electron scattering. Its marked features at vanishing energies, 
however, allows us to suggest that it could additionally strongly interact with dif- 
ferent atomic Rydberg species, if present in the same environment, thus actively 
taking part in charge-transfer chemical processes at fairly low energies. 
Furthermore, the present results provided by the angular redistributions of the in- 
cident electrons, definitely show that for all the contributing IRs of this molecule, 
the scattering is dominated essentially by forward ejection of the electrons and 
by negligible intensities in the angular regions beyond ~ 9 = 60 ^. Hence, the 
corresponding momentum-transfer cross sections, which are important indicators 
of electron deflection gas-phase molecules, turn out to be smaller than the cor- 
responding integral cross sections calculated here for the same system with the 



same quantum method (as depicted in figure 3.27 1. In conclusion, the correspond- 
ing rates of electron-molecule interaction are seen to be smaller than those orig- 
inated from the aint cross sections and therefore, in a statistical sense, indicate 
the substantial inefficiency of the coronene molecule for deflecting electrons. We 
thus argue here that more in general the coronene molecule is likely to have a lim- 
ited role in disturbing radiowave propagation through interstellar plasmas since it 
should not participate in producing conductivity inhomogeneity. 

To conclude, we want to underline that all the present findings can only be had 
from a realistic treatment of the quantum collisions between coronene molecule 
and environmental electrons, as performed in the present study. 
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3.3 The carbon nitride molecule, NC2N 

3.3.1 Introduction: role of metastable dicyanogen anions in the 
negative ion chemistry of Titan's upper atmosphere 

Titan has been the object of very substantial scrutiny since the arrival of Cassini- 
Huygens at the Saturn system in July 2004. Titan has been found to be the only 
satellite of the Solar System with an extensive atmosphere largely composed of 
N2, with abundant minor components like CH4 and H2 II112II . The presence of 
negative ions in its atmosphere had originally been considered only within the 
cosmic-ray-induced ionosphere which lies below 200 Km, since the pre- Voyager 
models [|113[ II 1411 considered their formation, by three-body electron attachment 
in the nighttime ionosphere, strongly reduced in the daytime because of the large 
photodetachment rates HUSII . 

It therefore came as a surprise when the Electron Spectrometer (ELS), one of 
the sensor of the Cassini Plasma Spectrometer (CAPS), revealed the existence of 
numerous negative ions in Titan's upper atmosphere Hl 16[ [TTTll . and furthermore 
negative ions have been detected on all Titan encounters when the spacecraft alti- 
tude was low enough and pointing conditions were favourable [jllSl . 
In order to understand the presence of negative ions in Titan's atmosphere, the set 
up of an ionosphere chemistry model for the formation of lighter-mass negative 
ions was recently presented and discussed [I119II . The dominant presence of CN^, 
followed by C3N and further down C5N" had been the result of the CAPS-ELS 
measurements at the altitude of about 10^ Km HUSH , where ions with a polyyne- 
like structure (C4H~ and CgH ) were also found to be fairly abundant. The mod- 
eling included a description of the suprathermal electron production from solar 
flux photoionization that presented two broad peaks around 1 eV and between 4 
and 6 eV, with the flux tapering down by several orders of magnitude beyond 10 
eV [I119II ; it also included a description of the neutral constituents distribution, 
producing at the end an abundance of the same lighter-mass anions found by the 
observation UllSL It did conclude, however, by indicating the desperate need of 
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knowing data on the radiative recombination of molecular radicals after the pri- 
mary electron impact processes. 

In the present study we shall therefore focus on one possible source for the for- 
mation of the most abundant anion observed in the Titan upper atmosphere, the 
CN~ anionic species, i.e. on the role of electron resonant attachment to NC2N and 
the possible pathways leading to the stable formation of both CN~ and (NCCN)^, 
another observed anion [I119L The latter has also been observed as a trace con- 
situent neutral species within the upper atmosphere of this large satumian moon 
[|1201I121]| . hence the existence of some amount of dicyanogen, which can be ei- 
ther depleted after electron attachment by the decay into CN~ or survive as the full 
anionic (NCCN)^, is an intriguing possibility that we shall explore in the present 
study. 

The present analysis is also connected with our previous published work on a 
different anionic species of polyyne-like structure, the dicyanoacetylene NC4N 
ni22|| . where our calculations demonstrated the likely hood of stabilizing the above 
anion by either radiative emission or by internal vibrational rearrangement (IVR) 
paths, possible mechanisms after the formation of a series of transient negative 
ions (TNIs) intermediates. It was found there that, besides the NC4N^ anion, 
other species like CN^ and CaN^ (mentioned before as being the most abundant 
anions in Titan's upper ionosphere [ |1211 ^ could also be produced via the above 
stabilization mechanisms and have indeed been observed in laboratory experi- 
ments [fT23l . 

3.3.2 Electron scattering at the equilibrium geometry 

The neutral molecule N-electron wavefunction, at its equilibrium geometry {Rnc 
= 1.16 A and Rcc = 1.39 A), was initially obtained at the Hartree-Fock (HE) level 
by using a 6-31 l-i-i-G(3df,3pd) basis set expansion within the Gaussian 03 suite 
of codes [|T4]| . The computed polaryzabilities were = 30.37 a.u., 0:2 = 19.59 
a.u.: the experimental values are not available for comparison. The employed 
range of integration for the scattering wavefunctions was extended up to 100 A, 
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employing a multipolar expansion of the bound molecular orbitals up to L^a^. = 
60, as well as for the continuum electron, extending the computed multipoles of 
the e^ -molecule interaction up to Xmax = 120.The total grid of integration within 
the (r, 9, 0) space is given by 1344 x 84 x 324 points. The results for the two 



resonances of vr* symmetry are shown in the panels of figure 3.28 
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structure 
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Figure 3.28: Spatial representation within the molecular space of the metastable excess 
electron for the two computed resonances of vr* symmetry. The energy positions, Er, and 
widths, r^, are also reported. 

We report in that figure the real components of the complex, continuum one- 
electron wavefunctions associated with the dynamical trapping (shape resonances) 
of the incoming electrons by the target molecule. The formation of TNIs is seen 
to occur at two distinct energies, shown in the figure, and to have different widths: 
the lifetimes before detachment correspond to 2.7 x 10~^^s for the lower energy 
resonance and to 9.2 x lO^^^s for the resonance at higher energy. According to 
what has been surmised [I119II about the energetics of the electron flux in the Ti- 
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tan's ionosphere at 1015 Km, we see that the resonant state at lower energy falls 
within one of the maxima of electron flux of about lO^cm^^s^^eV^^, while the at- 
tachment at higher energy would occur within a flux range reduced by nearly two 
orders of magnitude. The metastable anion of n„ symmetry also suggests that the 
extra density from the (N+1) electron is mostly located onto the C-C single bond 
and on the end N-atoms, as expected from chemical intuition. The correspond- 
ing lowest unoccupied molecular orbital (LUMO) exhibits a very similar spatial 
shape, although located at a very different energy (~ 5.4 eV), and is also a vr* 
orbital. The metastable anion is therefore qualitatively described as binding the 
extra electron into that vacant orbital, although the more rigorous scattering calcu- 
lations place the vr* resonance at a much lower energy in the electronic continuum. 
To help with a pictorical view of the quantum chemical calculations for the vir- 
tual orbitals, we report in figure 4.27 the spatial probabilities of the wavefunctions 
for the highest occupied molecular orbital (HOMO) and for the first three virtual 
orbitals of vr*, n* and a* symmetries. 

We clearly see the strong spatial similarities between the resonant wavefunctions 



(real part) shown in figure 3.28 and the LUMO and virtual02 orbitals from quan- 
tum chemical calculations; the eigenvalues, however, are very different from the 
actual resonant positions obtained from our calculations. All wavefunctions in- 
deed show the charge accumulation along the C-C single bond for the extra elec- 
tron attached to the target molecule in the tt* orbital. 

Another interesting result could be had by allowing the single C-C bond to stretch 
by about 30% of its value up to 1.8 A: the scattering calculations now reveal the 
appearance of an additional a* resonance very near threshold (at about 0.12 eV) 
and with a rather long lifetime of about 4 x 10^'^s: the data reported by figure 
13. 301 describe in more detail the behaviour of this resonance. 

The spatial features of the resonant wavefunction associated with the cr-type ex- 
cess electron clearly indicate that now the additional charge is chiefly localized 
onto the terminal C=N groups of the linear system. As one further stretches the 
C-C bond, however, we find that the resonance disappears into a bound state, i.e. 
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Figure 3.29: Computed molecular orbitals (MO) for the NC2N at its equilibrium ge- 
ometry. From top to bottom: highest occupied MO (HOMO), lowest unoccupied MO 
(LUMO), 2nd virtual MO (virtual02), 3rd virtual MO (virtual03). 



the present resonance seems to correspond to a compound state of very-long life- 
time and with very little residual energy that can rapidly dissipate through a linear 
IVR process and thus form a bound anion. We can surmise, then, that it consi- 
tutes a possible candidate for the formation of a negative ion of the full molecule: 
(NCCN)~, with a slightly stretched geometry. 



It is also interesting to note here that the sequence of MOs from quantum chemical 
calculations (see the data in figure 3.29 1 shows also the presence of a a* MO after 
the two TT* virtual orbitals, but with a very different energy collocation. 
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Figure 3.30: Computed resonance of symmetry obtained by stretching the single 
C-C bond. The upper panel reports the real part of the scattering wavefunction, while the 
lower panel show the location in energy of the strong resonant peak. 

3.3.3 Stretching and bending effects on the tt* resonances. 

In order to follow more in detail the effects of bond deformations on the computed 
resonances, we carried out a series of calculations for the vr* resonances depicted 



in figure 3.28 The first set of calculations involved the symmetric stretching of 



the C-C bond, leaving both C=N bonds at their equilibrium values. 



Figure 3.3 1 reports the findings for both resonances: we show there the resonance 
locations Eres(R) and indicate at each grid point the associated immaginary com- 
ponent, i.e. the width of each metastable state. 

The calculations clearly indicate the dramatic difference of behaviour between the 



two resonances. If one examines first the vr* resonance, we see from figure 3.28 



that the excess-electron wavefunction has no nodes across the bond in question 
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Figure 3.31: Computed behaviour of the resonance energy positions as a function of 
the stretching of the C-C bond (in A). The vertical dotted line shows the values at the 
equilibrium geometiy. Solid line: ^Hg metastable anion; dashes: metastable anion. 

and therefore it stands to reason that, in this case, the resonance gets destabilized 
as the nuclei move apart and its lifetime decreases. In other words, any electron 
attachment into this TNI is seen to reduce its stabilization upon bond lengthening 
and to favour instead the electron detachment path: no stable anion can form by 
the single bond stretching, contrary to what occurred for the a* resonance. 

If we now turn to the vr* resonance, the situation appears as entirely different: from 



figure 3.28 we see first that the excess-electron wavefunction exhibits a marked 
nodal plane across the C-C bond, a sign of strong antibonding character. The 
stretching of the bond, therefore, now stabilizes the resonance and lengthens its 



lifetime, as clearly seen from the data on the top panel of figure 3.31 In other 
words, the formation of a vr* TNI, although occurring at higher energies and 
therefore likely to be caused by the less intense component of the suprathermal 
electrons produced by solar flux photoionization [I119II , is stabilizing the negative 
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ion upon bond stretching rearrangements. 

This difference of behaviour can also be seen even more clearly from the potential 
energy curves of the resonant molecular states reported in figure 3.32[ obtained 



according to eq 1 .37 
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Figure 3.32: Computed potential energy curves of the two resonant states created in 
the ^Ug (solid line) and ^XI^ (dashes) symmetries. The C-C bond distances are varied 
along the range shown. Vertical line: location of the equilibrium geometry of the neutral 
molecule. 



The data of figure [332| indicate that both resonances, as expected, do not lead to 
dissociative attachment (DA) processes since the considered cuts of the overall 
potential energy surfaces (PESs) do not lead to molecular dissociation. On the 



other hand, the computed lifetimes of the associated resonances (see figure 3.31 ) 
tell us that bond stretching stabilizes a negative ion for the tt* resonance while 
leads to electron detachment for the tt* resonance. Thus, we could surmise that 
our calculations indicate already two possible ways to stable anionic formation af- 
ter electron attachment: (i) via the a* resonant intermediate near threshold and (ii) 
via the vr* intermediate formation at higher energies. Both stabilization processes 
appear to possibly occur via IVR couplings for the initial excess energy of the 
electron. It is also interesting to note that, when the resonant electron wavefunc- 



tion is represented at the largest stretching value shown by figure 3.32 it clearly 
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exhibits a strong localization of the extra charge density over the two CN regions 
with no additional density along the bond. Since the physical IVR process is 
really occurring on a multidimensional complex PES, then one could argue that 
more complicated molecular deformations after stretching could turn out to be 
dissociative deformations and therefore lead to final (CN)^ detachment. 

With regard to the stable molecular anion, earlier and highly accurate quantum 
chemical calculations [112511 indicated that the most stable bound anionic species 
is the one given by the trans configuration of the N=C-C=N, with a C-C=N 
angle distorted from 180° down to about 150°, yielding a vertical electron affinity 
of 0.65 eV, not far from the experimental indication of 0.58 eV H 12611 . In order 
to see the effects on the tt* resonances of deforming the neutral molecule close 
to the ground state geometry of the bound anion suggested in ref I1125II . we have 
carried out scattering calculations in which the angle mentioned before is distorted 
from the hnear configuration {9 = 180°) towards the value of the stable anion {9 
= 154.7°) for the trans configuration [I125II : the bond distances were kept at Rcn 
= 1.124 A and Rcc = 1-393 A during the bending, as indicated by the quantum 
chemical study [I125L 



The results of the calculations are shown in the two panels of figure 3.33 for both 



vr* and vr* resonances. 

The data on the upper panel report the removal of the vr-degeneracy as the sym- 
metry lowers to the -D2/1 symmetry: the a* resonance component acquires now 
the a* -like character existing in the linear configuration, while the &* component 
now follows the original 7r*-like character of the linear geometry. The calcula- 
tions confirm here what we had already surmised from the stretching effects seen 



in figure 3.29 but now with an inverted behaviour: the higher resonance moves 
little in energy when the C-C bond does not change and only the orientation of the 
C=N group changes, leading instead to a less stable compound state with strongly 
reduced lifetimes in both components. In other words, we see that the bending de- 
formation preferentially causes detachment of the resonant electron attached at 
the higher energy (~7.7 eV). 
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Figure 3.33: Computed resonance positions and widths (bars at the grid points) as a 
function of changing the angle 9 of the trans configuration (in sert). See main text for 
details. Upper panel: deformation effects on the 5-type resonant components, a* and b*, 
when the degeneracy is r emoved. Lower panel: same evolution for the u-type resonance 
components. 



On the other hand, the data in the lower panel clearly show that the a* -like com- 
ponent of the lower resonance is dramatically affected by the bending motion, 
moves down to threshold energy and in our calculations becomes almost a bound 
state around 9 = 160°, i.e. not far the stable anion geometry (9 = 154.7°) given by 
the quantum chemical calculations [|125l . This shows that electron attachment to 
the partner molecule, via the lower energy vr* resonance of its equilibrium geome- 
try, can evolve into a stabilized negative ion that undergoes efficient IVR coupling 
mechanism via the bending degree of freedom: this possibility could be examined 



more pictorically from the data reported in figure 3.34 



We report in that figure the real energies of the (N-i-l)-electron metastable anions 



associated with the two resonances of figure 3.33 3. We further show the energy 
dependence on the same angle 9 for the case of the ground state neutral molecule. 
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Figure 3.34: Computed electronic energies of tlie (N+l)-eIectron metastable states re- 
ported in fig. 4.31b {Au, Bu) and of the (N)-electron neutral molecule {Ag) as a function 
of the deformation angle 9. The reference energy is that of the linear neutral {0 = 180°) 

In spite of the fairly simplified description of the quantum electronic state of the 
neutral molecule (given at the Unrestricted Hartree-Fock (UHF) level only), the 
general physical picture comes up rather clearly: as the system is deformed to- 
wards the optimal geometry of its stable anion [I125II we see, in fact, that the 
(j*-like resonant state moves down in energy and tends to become a bound state 
crossing the real energy curve of the corresponding neutral partner. By extrapo- 
lating the scattering data to negative energy we see that stabilization of (NCCN)^ 
starts around 155°, i.e. very close to the structural value for the bound anion II125II . 
If one further considers that the estimated adiabatic electron affinity is around 0.3 
eV [I126L then we see that our extrapolation reaches that value around 150°, a 
failry realistic estimate from scattering calculations of the bending angle of 154° 
of the stable anionic species given in [11251 . 

A further way to qualitatively confirm the stabilization of the anionic molecule 



when following the bending path depicted in figure 3.33 is to look at the spatial 



features of the resonant excess electron at the last geometry of our scattering cal- 



culations, i.e. at the 9 value of ~160°. Those features are given by figure 3.35 and 
describe the a* -like component of the w-type resonant state. 



174 



CHAPTER 3. RESULTS AND DISCUSSION 




Figure 3.35: Computed real part of the wavefunction for the 6* component of the reso- 
nant electron in symmetry at the bending angle value of ~160°. See text for details. 

We clearly see that the excess, resonant electron is distributed over all the atoms of 
the system, it has 'V„" character with respect to the plane of symmetry in the bent 
molecule and it does not exhibit any "bond rupture" feature perpendicularly to the 
C-C bond. In other words, we are producing a metastable precursor, at nearly zero 
energy, of the bent configuration of the (NCCN)^ bound anionic species. 
It is also interesting to note here that we classify under the TVR' heading a broad 
variety of processes where nonadiabtic coupling effects are induced by nuclear 
motion. For instance, the use of Renner- Teller nonadiabatic effects to explain dis- 
sociative recombination processes as discussed in [I127II for the HCO+ -i- e^ system 
is another instance where our nuclear-induced nonadiabatic energy redistribution 
is seen to be at work. 

3.3.4 Present conclusions 

The present work has investigated, at the nanoscopic level, the quantum dy- 
namics of the stabilization paths which could lead to the formation of stable, 
bound anionic molecules of (NCCN) after the attachment of environment elec- 
trons. The latter electrons span the energy range of the observed distributions of 
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free electrons produced by the solar photons in the Titan's atmosphere I1119II and 
the quantum description of the scattering processes between such electrons and 
the title molecule reveals that, at specific energies, the collision events cause for- 
mation of resonant compounds, i.e. of metastable anionic complexes temporarily 
bound within the continuum energy of the electrons. In particular, our present 
calculations found that a resonance is formed at very low energy and with a 
very large lifetime, indicating its likely stabilization as a bound state by internal 
vibrational redistribution (IVR) processes of the small excess energy. 

Additionally, we found two more resonant anionic states: a state around 1.3 
eV and a^Ug one around 8 eV. The former resonance exhibits a longer lifetime 
than the latter and, upon stretching the C-C bond, they follow different evolutions. 
Thus, the resonance of vr* character stabilizes into a very long-Uved metastable 
state with strongly reduced residual energy of that excess electron, while the one 
of TT* character moves up in energy and appears to follow preferentially an electron 
detachment path: its lifetime decreases upon bond stretching. In other words, the 
stretching of the single C-C bond provides an efficient nuclear coupling of the 
excess electron motion only in the case of the vr* resonant electron, where it can 
lead to the formation of a stable (NCCN)^ molecule. 

We have further investigated the effects of bending deformations on the above 
resonances and carried out calculations in which the C=N groups are bent sym- 
metrically from the linear configuration and move the molecule to a geometry 
suggested by earher calculations [I125II to be the ground state of its bound nega- 
tive ion. 

Our calculations find now that the reduction of the symmetry removes the de- 
generacy of the TT-resonances and that each component follows a different path 



(see our data in figures 3.33 1, 3.33 5). We see there that the (T*-like component 
CBu) of the 7r* resonant electron is now dramatically lowering its energy loca- 
tion and lengthening its lifetimes. In comparison with the potential energy curve 
of the neutral molecule (^^4^), our scattering calculations indeed suggest that the 
metastable potential energy component (real part) for the '^Bu anion now crosses 
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the neutral potential around a bent structure of 155° (see figure [334]), and there- 
fore stabilizes the bound anionic configuration suggested by the earlier quantum 
chemical calculations. In other words, our modelling of the dynamical evolution 
that follows electron attachment on a series of bent configurations for the title 
molecule indicates that the extra-electron of 6* symmetry (cr*-like wavefunction) 
couples very efficiently with that bending motion and can stabilize into a bound 
state of the anionic species. These findings therefore confirm the resonant 
complex to provide an efficient doorway state to the stabilization of the (NCCN)~ 
anion in its D2h deformed configuration with respect to the linear neutral. On the 
other hand, our calculations also suggest that metastable attachment of electrons 
to NCCN molecules does not lead to the formation of (CN)^ fragments via shape- 
resonance mechanisms along simple pseudo-ID molecular deformations, but may 
require instead a more comphcated set of nuclear motions. 



Concluding remarks 



The present thesis focuses on low-energy electron collisions with astrophysically 
relevant closed-shell molecules, therefore represents a theoretical/computational 

work which deals with an area placed at the boundary between (molecular) astro- 
physics, quantum collision thery, and of course theoretical chemistry. 

The findings presented in this work demonstrate that taking in consideration res- 
onant low-energy electron attachment processes associated not only to s-wave 
resonances, does provide many useful information. From the astrophysical point 
of view, such a treatment reveals infact to be exhaustive in the sense that it enables 
to locate many more channels to form different metastable anions associated with 
a given neutral molecule. 

In the framework of a rigorous ab initio approach, the ensuing characterization of 
the resonances by the 3D spatial maps describing the resonant wavefunctions en- 
ables to immediately visualize which are the likely reactive sites on the metastable 
anion, therefore allowing to also qualitatively surmise possible reaction paths in- 
volving that species. In line with that, keeping in mind that the stable anion forma- 
tion process by resonant electron collision is a unimolecular process, our pseudo- 
ID approach proves to be a powerful tool by which investigate the dynamical 
evolution of the resonant species; following the 'chemical intuition', when sup- 
ported by the necessary decreasing of both the real and the imaginary (longer 
lifetimes) component of the resonance energy as function of the choosen molecu- 
lar ID deformation, the 3D maps for different molecular geometries could be also 
considered as indicative of an IVR process in the sense that the localization of 
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the resonant electron on a specific fragment could suggest a non-ergodic IVR(| [), 
while, conversely, the 'broadening' of such spatial resonant electron density over 
a larger molecular region could suggest an actual ergodic multidimensional IVR. 
Although we are not able to estimate the time requested by the IVR occurrence, 
we can however take advantage of our resonance lifetime estimates to qualitatively 
infer if there is time enough for the IVR to be effective. 

Moreover, given a neutral molecule, which is astrophysically relevant, the qualita- 
tive comparison between the energies requested to the resonant species formation 
with the environmental conditions of the astrophysical region under investigation 
can then allow a deeper insigth into the feasibility of certain chemical networks; 
in this sense, one can further take advantage of the knowledge of the relevant ab 
initio integral, differential and momentum transfer cross sections to also estimate 
the metastable anion formation's rate coefficient, the latter being an important pa- 
rameter needed by astrophysicists involved in modeling the chemical networks in 
the ISM. 



The three molecular species that were computationally investigated about their 
behaviour under low-energy electron collisions, are the ortho-benzyne (0-C6H4), 
the coronene (C24H12), and the carbon nitride (NC2N), respectively. As already 
emphasized in the general introduction, due to their sizes, their peculiar structures, 
their chemical reactivity, their physical and chemical properties as well as accord- 
ing to astronomical observations, when available, the above three molecules be- 
long to different astrophysical contexts in the sense that each of them is linked to 
an important example of a specific astrophysical puzzle that appear to currently 
constitute an intriguing scientific challenge. 

The following important remarks can now de drawn from our present study: 
i) ortho-benzyne 

The computational investigation for this molecule focused on the possible role 
which could be played by the anionic ortho-benzyne as a very reactive intermedi- 

'^usually leading to fragmentation 
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ate in proto-planetary (PPN) atmospheres, where the physical conditions are suit- 
able to justify its actual existence. Due to the positive electron affinity that char- 
acterizes the parent neutral species, the 0-C6H4 anion could be thus also viewed 
as a reasonable aromatic precursor of larger species belonging to the polycyclic 
aromatic hydrocarbons' family. 

When considering the 0-C6H4 as frozen at its equilibrium geometry, the results 
indicate that this molecule indeed exhibits four fairly narrow resonances below 
10 eV, three of which can be classified as tt* resonances while one of them is a 
low-energy a* resonance exhibiting antibonding character across the triple bond. 
When compared with the benzene molecule, it shows lower energy resonances: 
since the excess energy carried by the extra electron has to come from the en- 
vironmental medium in such a stellar atmosphere, the attachment processes to 
benzene do require much higher electron temperatures than those needed for the 
formation of metastable anions of benzyne. Furthermore, the (o-C6H4~)* tem- 
porary negative ion is seen to have markedly longer lifetimes, thereby allowing 
more easily the dynamical coupling of the extra electron with the nuclear degrees 
of freedom which play a crucial role in the ensuing molecular break-up or also in 
possible radiationless non-dissociative stabilization paths. 

Keeping in mind that the benzene has a rather large negative electron affinity (- 
1.12±0.003 eV, BUl), it therefore appears as reasonable to argue about the possible 
surviving of the benzene as a neutral species against the more likely involvement 
of the longer lived negative ion 0-C6H4 in several subsequent chemical reactions. 
This provides a possible explanation about the presence of the former and, at the 
same time, the absence of the neutral counterpart [|44ll . 

The pseudo-ID quantum dynamics, which was employed to study the possible 
dynamical evolution once the metastable anion is formed, revealed the stabiliza- 
tion of the initial metastable anion by small ring deformations which cannot lead 
to ring-breaking effects, and, even more interestingly, the formation of a bound 
aromatic anion by small direct stretching of the triple bond, therefore suggesting 
the possibility of reaching the bound (0-C6H4) ground state by fast intramolec- 
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ular vibrational redistribution processes. The combination of the estimated rate 
coefficient value for the metastable ortho-benzyne anion with the possible fast ra- 
diationless stabilization of such a resonant anion, enabled us to suggest that the 
neutral 0-C6H4 might actually play a role in the circumstellar chemistry of CRL- 
618 and more in general in the PPN chemistry, if the neutral parent molecule is 
actually produced in that environment. In conclusion, in line with the above con- 
jectures, which are in turn proven to be realistic by the findings of the present 
quantum scattering calculations, it was then possible to also tentatively indicate 
the reactions with H/H+ and the following, albeit qualitative, chemical paths for 
an aromatic condensation reaction toward larger PAHs, 



H 




CH 



where the molecules chosen in the above figure indeed correspond to some of the 
most likely abundant carbonaceous species in the CRL-618 environment ( [|47l ). 
ii) coronene 

The present calculations for the coronene molecule, a larger partner when com- 
pared with the previous molecular system, are indeed the first example of a com- 
putational investigation of a medium-size member of the PAHs family, one which 
can be considered as likely to be present in several different astronomical envi- 



ronments [HOSilTQll where, together with larger PAHs, plays a role in contributing 
to both their chemical and physical properties HlOSi I102[ [82l iTQll . Tackling such 
a 'large' species with reliable ab initio methods has furthermore meant to provide 
a link with earlier studies on smaller members of that wide family of molecules, 
and therefore to clarify the relative importance of at least two emblematic sys- 
tems: ortho-benzyne and coronene. At the same time, one can also speculate that 
such an endeavour could be useful to better extrapolate the behaviour of polycyclic 
aromatic molecules interacting with low-energy electrons, moving from the small- 
size closed shell PAHs (benzene, benzyne) toward the large-size closed shell PAHs 
(with more than 25-20 C atoms): the coronene, in fact, as a medium-size PAH 
without permanent dipole, can be viewed as representing one possible ideal link 
between the previous two PAHs classes. In this connection, we should note that 
one important point that makes intriguing the comparison between the coronene 
and the o-benzyne, is provided by the electron affinity for these molecules; despite 
their marked difference in size (thus in the number of normal modes), as well as 
the presence of a rather large permanent dipole that characterizes the o-benzyne 
(1.68 D, [193), the EA value is similar: about -1-0.5 eV for C24H12 (-1-0.54 ± 0.1 
eV according to the kinetic measurements in II84I and -1-0.47 ± 0.09 eV as mea- 
sured by photodetachment experiments reported in [[89ll ) and -1-0.56 eV for 0-C6H4 
[I2l[l00l, respectively. 

The detailed search of the possible formation of threshold, metastable anions for 
the coronene molecule was carried out and discussed under the ligth of the most 
recent available experimental data: our findings point to the virtual state scattering 
(zero energy resonance), occurring for vanishing energies, as the likely doorway 
to the production of the parent metastable not dehydrogenated negative ion, conse- 
quently allowing us to also surmise that an important role might be played by the 
nonadiabatic couplings with the nuclear motions during the zero-energy electron- 
molecule encounter in order to access the stable bound state. Therefore we argue 
that the stable bound state, that experiments show to occur for an energy close to 
eV, could be accessed via the zero energy resonance evolution due to the nuclear 



182 



CHAPTER 3. CONCLUDING REMARKS 



motions that increase the induced dipole moments interaction. 

More in general, on the basis of the results obtained on its capability to form 
resonances in the 0-7 eV energy range at the equilibrium geometry, the coronene 
molecule could be viewed as a species that forms bound stable anions but almost 
only at very low-energy: according to experiments, in fact, at incident energies 
higher or equal to 2 eV, the present findings allow to surmise the occurrence of 
several resonant excitation processes instead of actual bound stable negative ions. 
From the astrophysical point of view, this means that the C24H12 molecule should 
not be able to strongly participate in soaking up a rather large amount of free 
electrons; on the other side we can speculate that it could strongly react with 
different atomic Rydberg species, if present in the same environment. 

In order to try to understand its possible behaviour, at the molecular level, as 
a deflector of low-energy electrons present in those astrophysical contexts, and 
therefore understand better its possible role as an energy absorber from the en- 
vironmental electrons or within a low-temperature molecular plasma produced 
in those astronomical surroundings, the elastic (rotationally summed) differential 
cross sections were investigated: the findings enable to look at this molecule as 
a relatively 'transparent' species for low-energy electrons, so that we can also 
conjecture that the C24H12 molecule is likely to have a limited role in disturbing 
radiowave propagation through interstellar plasmas since it should not participate 
in producing conductivity inhomogeneity. 

iii) dicyanogen 

This molecule is relevant for the chemical modeling of Titan's upper atmosphere, 
since according to the current reaction network [I119II , such a species is supposed 
to efficiently yield, by dissociative electron attachment, the most abundant an- 
ion observed in that environment: the CN^. In this connection, following the 
interaction with low-energy environmental electrons by means of a pseudo-ID 
modelling of excess energy redistributions, we have analysed in some detail the 
possible dynamical paths which could lead, under the conditions observed in Ti- 
tan's atmosphere, to the formation of both the above small anionic fragment and 



the stable (NCCN)^ . The latter, in fact, represents another possible product that 
could be efficiently produced, once the primary electron attachment have taken 
place so that the corresponding metastable anion, (NC2N^)*, is formed. 

The present scattering calculations revealed that a ^S^j resonance is formed at 
very low energy and with a very large lifetime, indicating its likely stabilization 
as a bound state by internal vibrational redistribution (IVR) processes of the small 
excess energy. Additionally, we found two other resonant anionic states, a 
state around 1.3 eV and a ^Ug one around 8 eV, respectively, which upon stretching 
the C-C bond follow different evolutions: the former, which at the equilibrium 
has a longer lifetime than the latter, moves up in energy with the C-C lengthening, 
and, as its lifetime decreases, it also appears to follow preferentially an electron 
detachment path. The other one, conversely, appears to stabilize into a very long- 
lived metastable state characterized by a strongly reduced residual energy, so that 
in this case we are able to argue that the C-C stretching should provide an efficient 
nuclear coupling by means of which the formation of the stable (NCCN)^ anion 
can occur. 

The CN^ fragment production is not however completely excluded: according to 
our pseudo-ID dynamical investigation, we can only argue that more complicated 
molecular deformations, following the C-C stretching, are likely to be involved in 
causing this fragmentation. 

We have further investigated the effects of bending deformations on the above 
resonances and carried out calculations in which the C=N groups are bent sym- 
metrically from the linear configuration and move the molecule to a geometry 
suggested by earher calculations ni25il to be the ground state of its bound nega- 
tive ion. 

Our calculations find now that the reduction of the symmetry removes the degen- 
eracy of the TT -resonances and that each component follows a different path (see 
our data in figures 4.31a, 4.31b). We see there that the cr*-like component (^By) 
of the TT* resonant electron is now dramatically lowering its energy location and 
lengthening its lifetimes. In comparison with the potential energy curve of the 
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neutral molecule CAg), our scattering calculations indeed suggest that the the real 
part of the metastable potential energy component for the ^-B„ anion now crosses 
the neutral potential around a bent structure of 155° (see figure 4.32), and there- 
fore stabilizes the bound anionic configuration suggested by the earlier quantum 
chemical calculations [I125L 

As a general conclusion from the astrophysical point of view, the above three 
in-depth studies thus provide important examples on the role of anionic species 
in different astronomical environment like a PPN atmosphere (where we surmise 
the involvement of o-benzyne anion in PAH synthesis), a dense interstellar cloud 
(where the coronene could represent a typical medium- size PAH) and a planetary 
ionosphere (the upper atmosphere of Titan where the anionic chemistry of nitro- 
genated polyynes as dicyanogen play a role) where, without any claim of com- 
pleteness, each of the above three molecular systems investigated in this work 
represents in fact an emblematic example. 



Appendix A 



The elastic ICS and DCS 



As briefly introduced in section 2.4.8, in order to obtain an expression for the dif- 
ferential cross section one needs to introduce then a coordinate system S' which 
is space-fixed (SF), the so-called LAB system, and another coordinate system 
S which is fixed to the molecular highest symmetry axis and that is identified 
as body-fixed (BF). The latter is obtained from S' by a rotation through the ap- 
propriate Euler angles (a, (3, 7). The previously discussed expression for the 
FN scattering amplitude (see eq. 2.66), through the symmetry adapted functions 
X^^(ko) and X^^(r), contains angular factors like Yim{^o) and Ygmi^') that are in 
turn defined with respect the BF system, S (kg and r are defined with respect the 
S reference). The corresponding expression in the SF frame is obviously given by 



where r is defined with respect to 5" and r with respect to 5* and the coefficient 
^mm' ("^^ 7) '^he Wigner functions [9]. So, to this purpose, one has to transfer 
the angular terms from the BF to the SF frame through the following relations 




(A.1) 




(A.2) 



m' 



X 




(A.3) 



m 
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where, since in the BF frame the incoming electron points toward the molecular 

main symmetry axis which is in turn taken as the z axis, in the last expression one 
should consider -Q-Q) and A'=0, therefore obtaining 

^.T(k;) = E 6?L^^L(«/37) (^) (A.4) 

By substituting eqs. A.2 and A.4 in eq. 2.66, the SF scattering amplitude is given 
by 



/™(kV . r'; = i: E Z E ^^^^i'-" (A.5) 
which, remembering the definition of T matrix as 1 — S becomes 



/..(kV ■ r'; a/37) = E E E E '^"^ ^'^ ^^"^ (A.6) 

ehme'h'm' p/i A ° 



The differential cross section (thus the angular distribution) for the scattering of 
an electron by randomly oriented molecules is then obtained by averaging the SF 



quantity 



over all the molecular orientations (a, /3, 7): 



da 



{k, . f') = ^ y dadpdj sin PlfsFiK ■ r';a(3j)\^. (A.7) 



It therefore follows that 



\fMK-r': a/jT)i^ = EE E E EEEp (^.s) 

ehm Ihfh t'h'm' eh'fh' P/* A,A ° 

^Am' Am' "^Om ^Om ^ih/'h'^ thfi h' 
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The last expression, when integrated as in eq. A.7, can be simplified by making 
use of the contraction theorem foe the spherical harmonics (see eq. 3.7.72 p. 216 
inHni) 



, . v-v- /(2f + 1) (2f + 1) (2L + 1) [ e F L \ e i L , 

LM V 47r \^ A -A M y \^ C 

as well as by using the following relations (see p. 216 eq.3.7.69 in [I128II and p.103 
eqs. 3.116-3.118 in [fT29l) 



L=\t 



l+£ L L 

M=-L M'=-L \ m m M J \ m' m' M' 



m' m' m' I \ m m m 



SO that the angular integration of the four rotation matrices in eq. A. 8 provides 



Sin /3 dadpd^ Vl^,Vi;^,V'orJ^L = (A.9) 



l+l +L 



L=\l-l\ M=~L 



\ m -m M 



l+l +L 

E E 

L=\i-i\ M=-L \ J \m -m M 



Y: E (2^ + 1) (-If I ^ ^ - I I ^ ; ^ I jdnvi^,vli,^ 



l+l +L 



L \ i i L 



Y E (2^+i)(-ir . 

L=\i-i\ M=-L \ I \m -m M 



SttM (-1) 

A -A / \ m' -m' -M 
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Noting now the orthogonality between 3j coefficients as expressed by writing (see 
eq. 2.32 p.51 lfT29H ) 



i' £' L \ e F L . , 
Vt \ a -a m y \^ a -a ' 

it is finally possible to write 



da 



+0O 



{k,-r') = J2AcVc{cos^), (A. 10) 

where Vcicosi!)) are Legendre polynomials and the coefficients Ac are therefore 
found to be 



^^ = i^(2^+i)'EEE E EE 7(2^+1) (2^>i) m'+i){2i') 

ehm ihm i'h'm' e'h'm' PfJ- Pf^ 



i i C \ i' i' C \ i i £ \ / f C 

00 J \o J \m -rh M J \m' -m' -M 

(_'\\fh+m' flVfl rrppfj, 1V^l \ ( 1pp.* rrppp* iVp* \ ^ A 1 n 

The total integral (cTj) and the momentum transfer ((t„J cross sections, (ICS 
and MTCS, respectively), follow from eqs. A. 10 and A.l 1: 



(Ti 



/d 
-^iil^-v)s\Yi'd d'dd^ = ATiAo (A.12) 



CT^ = y ^ (k'o ■r){l- cos ^) sin ^ dMip = 4n (^Aq - ^^i^ • (A- 13) 

In conclusion, the expressions for and Ai can be directly deduced from eq. 
A. 1 1 ; as an example, is given by 



^0 = E E E Kn.^z,i'H'^A'J" i-^r^"-' (A.14) 

ehm e'h'm' pfi 



Appendix B 

The rotational excitation cross 
sections 

In appendix A the way to obtain a rotationally summed cross section is pro- 
vided. In practice, however, it can be efficiently computed by sums over magnetic 
substates as briefly illustrated in the present appendix; see IfTOll for a full detailed 
treatment. 

So far the participation of the rotational (and vibrational) molecular degrees of 
freedom were ignored during the scattering process. In spite of the added com- 
plexity that such additional couplings introduce in the treatment of the collision, it 
is however still possible to extract some useful information by using the approx- 
imate treatment of the dynamics known as the adiabatic nuclei approximation 
(AN). 

Briefly, the physics involved requires to consider the effect of the kinetic energy 
operator of nuclear motion as negligible with respect to that of the corresponding 
operator of the electronic motion. As a consequence, the molecule is seen as hav- 
ing 'infinite' mass and to stand still during the scattering event, when the energy 
spacings of the molecular levels is treated as negligible with respect to the local 
kinetic energy of the impinging electron. In other words, the AN approximation 
is a special case of the FN approximation; in fact, if the nuclei are allowed to 
be fixed (not necessarily in their equilibrium position), it is possible to neglect 
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the coupling between the incoming particle's angular momenta and the rotational 
(and/or vibrational, in general) quantum number of the molecule. It therefore fol- 
lows that the electrons adiabatically respond to the istantaneous position of the 
nuclei. In the absence of long-range forces (non-polar molecular targets) this is 
indeed a very good approximation, but for polar molecules the coupling between 
the projectile and the nuclear motion must be taken into account: it can be re- 
laxed in the core-region (in the inner-region) but in the outer region the nuclear 
Hamiltonian should be incorporated into the scattering equation. 
Assuming the BO approximation is valid, the AN total scattering (electron -i- 
molecular target) wavefunction can be written as: 

^(re,r,R^) = Mre,r;RN)^N(RN) (B.l) 

where 'n' represents all the rotational (and vibrational) quantum numbers of the 
unperturbed nuclear wave function. The adiabatic nuclei rotation (ANR) approxi- 
mation for the scattering amplitude is obtained for the transition rij — > n^, follow- 
ing [2301, as 

^n) (B.2) 

where Ui and nj specify only the rotational molecular quantum numbers and 
fpNi^'o ■ r'; Kj\f\al3'y) is the FN elastic scattering amplitude in the SF frame (see 
sec. 2.4.8 and appendix A). In this connection, one should note that the collision 
time Tcoii may be comparable to the vibrational period; however, the rotational 
period Tj-ot ma still be such that Tcou ^ Trot, where this condition is satisfied even 
at thermal energies since for low-lying rotational states the rotational spacing is 
of the order of 10^^ ~10^'^ eV. Full details about the specific expressions for 
asymmetric top, symmetric top and spherical top molecules may be found in [|T0| . 
In the case of a symmetric top target molecule, the rotational wavefunction is 
given by the well-known formula 
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where M gives the projection of the total angular momentum J along the 'z' axis 
within the SF frame, while K gives its projection along the molecular figure axis 
that in turn corresponds with the 'z' axis in the BF frame; in this case, both M 
and K are good quantum numbers. 

One therefore writes the explicit scattering amplitude as calculated from eq. B.2, 
/{JiKiMi JfKfMf) and obtain from it the differential cross sections summed 
over all final degenerate states Mf and averaged over all initial degenerate mag- 
netic substates Mi 



where k''^ = k"^ + 2{Ej^Ki — Ej^Kj) and the formula above is appropriate for 
molecules like C24H12. 

In the case of asymmetric top (like for the 0-C6H4), K is no longer a good quan- 
tum number and the rotational wavefunction needs to be expressed as a linear 
combination of asymmetric top functions: 

K=-J 

where the coefficients a'^j^ are determined by diagonalizing the asymmetric top 
rotational Hamiltonian and can be chosen so that the asymmetric top wavefunc- 
tions are orthonormal and r is a pseudo quantum number which differentiate the 
different rotational levels with the same J . The scattering amplitude can there- 
fore be had once more from the B.2 formula and the corresponding cross section, 
summed over the final states Mf and averaged over initial states Mj is given by 



— ( T T - % 

dn^''''^'''^^-h{=ko)^^^^2j, + i 



E 



fiJiTiMi ^ JfTfMf) 



(B.6) 
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where in the present case k'"^ = fc^ + 2{Ej.^. — Ej^^-^). 

The corresponding rotational wavefunction for a spherical top molecule is the 
same as that for the symmetric top. There is a further degeneracy in that the 
(2 J + Vf states with the same J, but with different M and K, have the same 
energy Ej, where 



Ej^BJ{J + l) 



(B.7) 



where B is the molecular rotational constant (in unit of energy). The correspond- 
ing differential cross section is therefore given by 



da 



- - M^M L ^ l^<'-''''* " ''''' 



Mi,Mf 



(B.8) 
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